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FOREWORD 


To  ensure  that  the  U.S.  Army's  soldiers  acquire  the  skills 
and  knowledge  necessary  to  perform  their  jobs  successfully,  the 
U.S.  Army  Research  Institute  for  the  Behavioral  and  Social  Sci¬ 
ences  (ARI)  performs  behavioral  research  to  develop  methods  of 
training  that  can  improve  skill  acquisition.  Morse  code  copy 
training  is  one  skill  area  that  continues  to  challenge  the  re¬ 
search  community  to  identify  training  strategies  that  can  alter 
the  training  attrition  pattern,  particularly  during  the  speed 
building  phase. 

This  report  describes  a  quantitative  model  that  simulates 
the  perceptual-motor  skill  responsible  for  successful  Morse  code 
copy.  This  model  will  allow  detailed  characteristics  in  the 
underlying  human  information  processing  mechanics  to  be  simulated 
and  compared  to  actual  performance.  Further  application  of  this 
model  might  enable  the  "at-risk"  students  to  be  identified  early 
in  the  speed  building  phase  of  training  and  thus  designated  for 
specialized  training.  In  general,  this  research  adds  to  our 
understanding  of  the  acquisition  of  skilled  performance. 


EDdp^  M.  JOHNSON 
Acting  Director 
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MODELS  OF  MORSE  CODE  SKILL  ACQUISITION:  SIMULATION  AND  ANALYSIS 


EXECUTIVE  SUMMARY 


Requirement : 

The  U.S.  Army  Intelligence  School  at  Fort  Devens  (USAISD) 
has  been  experiencing  an  unacceptably  high  rate  of  attrition. 

(The  attrition  rate  was  recently  reduced  after  the  introduction 
of  a  new  training  device.)  The  Assistant  Deputy  Chief  of  Staff 
for  Training,  HQ  Training  and  Doctrine  Command  (TRADOC) ,  re¬ 
quested  that  the  U.S.  Army  Research  Institute  for  the  Behavioral 
and  Social  Sciences  (ARI)  examine  the  problem  and  pursue  a  course 
of  research  that  could  lead  to  a  reduction  in  the  attrition  rate. 
The  simulation  and  analysis  of  performance  and  models  to  differ¬ 
entiate  between  successful  students  and  attrites  were  important 
steps  in  the  approach. 


Procedure : 

In  a  briefing  to  HQ  TRADOC,  a  cognitive  process  model  of  the 
Morse  copy  skill  was  proposed  as  a  means  to  understand  the  intri¬ 
cacies  of  the  skill.  An  elaboration  of  this  model  was  briefed  to 
USAISD  staff.  Through  a  contractual  arrangement  with  the  U.S. 
Army  Research  Office,  researchers  developed  simulation  models 
that  applied  tools  in  queuing  networks  and  order-of-processing 
diagrams  to  the  task  of  a  student  learning  to  copy  Morse  code  at 
progressively  faster  speeds. 


Findings: 

A  simulation  model  was  derived  from  a  cognitive  analysis  of 
the  information  processing  demands  on  students  and  modeled  with 
order-of-processing  diagrams.  Various  assumptions  about  the 
••copy  behind"  phenomenon — not  responding  until  a  subsequent  char¬ 
acter  is  presented — were  made.  The  model  predicts  the  probabil¬ 
ity  of  a  correct  response,  an  incorrect  response,  a  period  (no 
guess)  response,  no  response,  and  a  correct  response  for  each  of 
the  five  serial  positions  in  a  group.  The  simulation  also  pre¬ 
dicts  the  time  it  takes  to  execute  both  a  correct  and  incorrect 
response  and  the  time  it  takes  to  execute  a  period  response.  A 
preliminary  test  was  conducted  with  response  data  from  students 
at  the  U.S.  Army  Intelligence  School,  Fort  Devens,  early  in  the 
speed  building  phase  of  training.  The  models  fit  the  observed 
data  quite  well,  with  only  the  predicted  response  times  for 
incorrect  responses  much  different  from  the  observation. 
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utilization  of  Findings: 

These  models  will  be  tested  with  data  from  students  who 
either  acquired  the  skill  rapidly  or  attrited  during  the  speed 
building  phase.  This  will  establish  whether  response  patterns 
early  in  the  speed  building  phase  can  distinguish  the  quick 
learners  from  the  likely  attrites.  These  findings,  in  turn,  will 
be  presented  to  USAISD  for  consideration  in  identifying  students 
for  specialized  training. 
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MODELS  OF  MORSE  CODE  SKILL  ACQUISITION: 
SIMULATION  AND  ANALYSIS 

INTRODUCTION 


Background 


The  U.S.  Army  Intelligence  School  at  Fort  Devens,  Massachusetts,  trains 
Morse  code  copying  to  about  1,000  individuals  per  year  from  the  four  armed 
services.  The  goal  is  to  train  military  personnel  to  copy  Morse  characters  at 
the  rate  of  20  groups  per  minute  (gpm;  1  group  equals  5  characters)  with  96 
percent  accuracy.  The  training  has  historically  had  a  high  rate  of  attrition, 
at  times  over  30  percent  (DOES,  1990) .  For  those  who  succeed  in  the  basic 
Morse  course,  the  training  time  has  a  wide  degree  of  individual  variability, 
ranging  from  200  to  1180  hours  to  achieve  proficiency.  By  far,  most  training 
time  is  spent  in  the  speed  building  phase  where  students  try  to  progress  in¬ 
crementally  from  6  gpm  to  20  gpm.  Since  attrition  is  most  likely  to  occur 
during  speed  building  and  well  into  the  training  program,  its  effect  can  be 
costly. 

Although  the  learning  of  Morse  code  has  long  been  a  topic  of  research  in 
psychology  (Bryan  and  Harter,  1897),  the  high  training  attrition  persists. 

The  learning  of  individual  Morse  characters  is  fairly  straightforward,  but  the 
ability  to  copy  groups  of  characters  accurately  at  increasing  rates  of  trans¬ 
mission  is  difficult  to  master  for  many  individuals,  even  after  hundreds  of 
hours  of  practice.  It  is  this  difficulty  in  being  unable  to  progress  to  a 
higher  speed,  not  the  lack  of  motivation  or  practice,  that  leads  to  the  costly 
attrition. 

The  copy  task  requires  listening  to  a  sequence  of  Morse  code  characters 
(each  character  is  itself  a  sequence  of  dahs  and  dits)  and  responding  with  the 
corresponding  character  on  a  standard  keyboard.  The  ability  to  copy  Morse 
code  characters  during  rapid  rates  of  transmission  (100  characters  per  minute) 
is  clearly  demanding.  To  copy  at  speeds  approaching  20  gpm  burdens  the  under¬ 
lying  cognitive  processes,  pushing  relevant  memories  to  their  limits.  It 
requires  the  simultaneous  processing  of  different  stimuli,  frequently  respond¬ 
ing  to  one  character  while  another  is  being  perceived.  Perhaps  the  memory 
capacity  of  some  operators  is  exceeded,  or  perhaps  the  concurrent  processing 
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abilities  of  other  operators  are  pushed  too  far.  Present-day  computer  technol¬ 
ogy  and  theoretical  frameworks  for  modeling  cognitive  processes  suggest  that  a 
closer  study  can  be  made  of  the  development  of  proficiency  during  speed  build¬ 
ing  than  was  previously  possible.  The  results  of  this  more  extensive  analysis 
can  then  lead  to  prescriptions  for  overcoming  the  persistent  training  barrier. 

A  closer  study  of  the  development  of  proficiency  during  speed  building 
requires  a  detailed  understanding  of  how  accuracy  and  response  time  change  as 
a  function  of  the  rate  of  presentation  (groups  per  minute) ;  the  serial  posi¬ 
tion  of  a  character  within  a  group;  the  serial  position  of  a  character  within 
a  block;  the  identity  of  a  character;  the  confusability  of  each  character  with 
all  other  characters;  the  times  to  perceive,  recognize,  and  program  a  response 
to  a  stimulus;  the  sizes  of  the  various  memories  involved  in  processing;  and 
the  "decay"  times  of  the  items  in  these  memories.  This  detailed  understanding 
can  only  come  from  a  precise  specification  of  the  cognitive  activity  underly¬ 
ing  Morse  intercept  behavior. 

Qualitative  Models.  Toward  this  end,  several  models  of  processing  in 
the  Morse  intercept  task  have  been  proposed  (Sabol,  Wisher  and  Medici,  1991; 
Townsend,  1991;  Wisher,  Kern  and  Sabol,  1990) .  These  models  all  derive  from 
the  two  descriptions  of  processing  first  presented  by  Wisher,  Kern  and  Sabol 
(1990).  Briefly,  in  their  early  state  cognitive  process  model,  it  was  assumed 
that  subjects  could  execute  only  one  process  at  a  time  and  that  the  buffers  in 
front  of  the  processors  were  limited  to  one  item  at  most.  In  the  advanced 
state  model.  Wisher  et  al.  assumed  that  subjects  could  execute  several 
processes  simultaneously.  In  addition,  they  assumed  that  the  buffer  sizes  were 
greater  than  one.  We  want  to  pay  particular  attention  to  this  advanced  state 
model. 


In  the  advanced  state  model  four  processes  were  arranged  in  series: 
auditory  perception,  character  recognition,  motor  organization,  and  response 
execution.  In  front  of  each  processor  was  a  buffer  that  held  information 
as  it  was  passed  along  through  the  system.  Each  buffer  had  a  limit  on  the 
number  of  items  that  it  could  hold  and  items  in  each  buffer  could  decay  over 
time.  Here  decay  refers  to  the  natural  loss  of  information  if  it  is  not  rehearsed. 
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Note  that  the  processing  time  for  each  server,  decay  times  for  items  in  each 
buffer  and  capacity  limits  of  each  buffer  were  left  unspecified. 

At  the  outset  of  processing,  the  stimulus  (a  sequence  of  square  waves  of 
sound  energy  or,  more  simply,  a  tone  sequence)  was  placed  in  the  auditory 
buffer  (echoic  memory) .  In  the  auditory  perception  stage  (which  accepted  in¬ 
formation  from  the  auditory  buffer) ,  it  was  assumed  that  the  sound  energy'  was 
transformed  into  a  string  of  elements  (i.e.,  that  sound  energy  was  transformed 
into  a  string  of  dits  and  dahs) .  These  element  strings  were  then  passed  to 
the  recognition  buffer.  In  the  recognition  stage  (which  accepted  items  from 
the  recognition  buffer) ,  it  was  assumed  that  an  element  string  was  recognized 
as  a  character.  The  character  code  was  passed  on  to  the  motor  organization 
buffer.  It  was  then  assumed  that  the  motor  program  needed  to  typ'=  the  par¬ 
ticular  character  was  set  u^  by  the  motor  organization  stage.  Finally,  the 
motor  program  was  passed  on  to  the  motor  execution  buffer  where,  eventually, 
the  response  was  initiated.  We  should  note  that  the  above  reflects  the 
details  of  just  one  of  the  late  stage  models  proposed  by  Wisher  et  al. 

Several  variations  on  the  above  were  suggested. 

Analytic  Models .  The  above  qualitative  model  suggests  a  number  of  points 
in  the  processing  of  stimuli  where  the  operator  may  have  difficulty.  However, 
it  ic  difficult  to  test  the  various  hypotheses  without  being  able  to  predict 
how  accuracy  and  response  time  change  a"  a  function  of  the  changes  in  che  in¬ 
dependent  variables  mentioned  above.  Towards  this  end  it  is  necessary  to 
quantify  the  model.  Unfortunately,  this  is  not  by  any  means  an  easy  tas)c. 

Briefly,  the  model  proposed  by  Wisher  et  a].  (1990)  is  quite  clearly  a 
queueing  model.  The  literature  on  queueing  models  is  extensive  (e.g..  Gross 
and  Harris,  1985,  discuss  the  basic  elements  of  queueing  theory;  Rouse,  1980, 
reviews  queueing  networks  that  have  been  used  to  model  person-machine 
systems) .  However,  the  majority  of  the  existing  analyses  assume  that  the 
queueing  system  has  reached  steady  state  or  equilibrium.  Such  will  not  be  the 
case  for  the  behavior  of  the  queueing  system  in  the  tasks  undertaken  by  the 
Morse  intercept  operator.  In  these  tasks,  the  system  is  unlikely  to  reach 
equilibrium  since  five  tone  sequences  are  presented  one  at  a  time,  followed  by 
an  interval  of  time  which  in  many  cases  is  long  enougr  for  the  system  to  clear 
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itself.  And  even  were  the  system  not  to  clear  itself,  the  fact  that  a  break 
exists  between  groups  makes  it  impossible  to  accept  the  assumption  of  steady 
state.  Thus,  one  needs  to  undertake  a  transient  analysis  of  the  queueing  sys¬ 
tem. 

Townsend  (1990)  building  on  results  reported  in  Fisher  and  Smith  (1987) , 
shows  the  form  that  the  analyses  must  take.  Basically,  it  is  necessary  to 
translate  the  queueing  network  into  an  Order-of-Processing  (OP)  diagram 
(Fisher  and  Goldstein,  1983;  Goldstein  and  Fisher,  1991) .  Once  so  translated, 
it  is  possible  to  obtain  closed  form  expressions  for  relatively  restricted 
queueing  networks.  The  analysis  for  two  networks  was  worked  out  in  detail  by 
Townsend. 

Briefly,  in  the  first  network,  there  were  two  nodes.  A  node  consists  of 
a  server  and  a  buffer  (possibly  of  size  zero) .  Townsend  assumed  that  the 
server  at  the  first  node  encodes  and  recognizes  the  stimulus  and  the  server  at 
the  second  node  organizes  and  executes  the  response.  The  buffer  at  the  first 
node  was  unlimited  in  capacity.  The  buffer  at  the  second  node  was  of  size 
zero.  If  an  item  was  encoded  and  the  response  server  was  executing,  that  item 
was  held  at  the  encoder  and  the  encoder  was  prohibited  from  accepting  further 
input  until  the  response  server  completed  its  execution  of  the  downstream 
item.  Townsend  assumed  that  n  items  were  present  in  the  first  buffer  at  the 
start  and  that  no  new  items  were  added  to  the  system.  Finally,  Townsend  as¬ 
sumed  that  the  service  times  were  independent,  exponentially  distributed 
random  variables.  For  this  system,  he  derived  analytic  expressions  for  t:  j 

time  on  average  it  took  to  respond  to  the  i^^  item  in  the  first  buffer. 

The  above  system  does  not  contain  a  mechanism  for  producing  errors. 

Thus,  by  itself,  the  system  cannot  be  the  one  that  subjects  are  using  since, 
among  other  things,  subjects  always  produce  a  significant  number  of  errors. 

One  reasonable  way  to  generate  errors  is  to  assume  that  the  subjects  can  ex¬ 
ecute  at  most  one  process  in  the  system  at  any  one  time.  If  an  item  arrives 
while  either  the  encoder  or  responder  is  busy,  that  item  gets  lost  from  the 
system.  In  the  previous  model,  the  item  was  simply  held  at  the  encoder.  For 
this  new  system,  Townsend  was  able  to  derive  both  the  probability  of  an  error 


and  the  time  on  average  it  takes  to  respond,  given  that  a  correct  response  is 
made  (he  assumed  that  if  an  item  was  lost  from  the  system  e''ther  no  response 
was  made  or  an  incorrect  response  was  made) . 

Goal  and  Objectives 

Goal.  There  are  many  reasons  operators  may  find  it  difficult  to  complete 
training.  One  hypothesis  (Wisher  et  al.,  1990)  is  that  some  operators  cannot 
simultaneously  listen  to  (encode)  an  incoming  signal  at  the  same  time  as  they 
attempt  to  recognize  and  respond  to  other  signals  that  have  already  arrived. 
This  ability  is  referred  to  as  the  ability  to  copy  behind.  The  overall  goal 
of  the  work  undertaken  as  part  of  this  effort  is  to  determine  the  extent  to 
which  the  copy  behind  hypothesis  is  supported  by  the  data  which  has  been  col¬ 
lected  in  the  Morse  code  copy  task. 

Objectives.  Attainment  of  the  above  goal  requires  the  meeting  of  three 
objectives.  First,  in  order  adequately  to  test  the  copy-behind  hypothesis,  it 
is  necessary  to  model  Morce  code  copy  behavior.  In  particular,  it  is  neces¬ 
sary  to  have  a  detailed  model  of  the  processes,  memories,  and  other  cognitive 
building  blocks  involved  in  the  Morse  intercept  process.  As  noted  above. 
Wisher  et  al.  (1990)  has  proposed  a  rather  extensive  qualitative  model.  What 
is  needed  at  this  point  is  a  more  precise  characterization  of  what  exactly  is 
meant  by  copy-behind,  when  exactly  a  period  response  will  be  produced,  when  no 
response  will  be  made,  and  so  on  (Objective  1) . 

Second,  once  the  generic  models  are  available,  it  is  necessary  to  quan¬ 
tify  the  behavior  of  the  models.  In  particular,  we  were  interested  in  having 
the  model  predict  the  probability  that  an  operator  types  a  character  cor¬ 
rectly,  types  a  character  incorrectly,  types  a  character  as  a  period  or  fails 
to  type  a  character.  (Periods  are  typed  when  an  operator  realizes  that  a 
character  has  been  transmitted,  but  is  unable  to  respond  quickly  enough.  The 
period  serves  to  maintain  the  format,  i.e.,  the  spacing.)  And  we  were  inter¬ 
ested  in  having  the  model  predict  the  time  that  it  takes  an  operator  to  type  a 
character  correctly,  the  time  that  it  takes  an  operator  to  type  a  character 
incorrectly,  and  the  time  that  it  takes  an  operator  to  type  a  period.  As 
above,  Townsend  (1990)  has  produced  analytic  expressions  for  response  times 
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and  errors  for  several  more  simple  models  than  the  one  developed  by  Wisher. 
What  is  needed  at  this  point  are  predictions  of  the  dependent  variables  for 
the  more  complex  model .  This  can  be  achieved  using  either  simulation 
(Objective  2a)  or,  where  possible,  deriving  the  relevant  analytic  expressions 
(Objective  2b) . 

Finally,  given  that  the  behavior  of  the  models  can  be  quantified,  it  is 
necessary  to  determine  just  how  well  the  models  fit  the  data.  None  of  the 
models,  either  the  qualitative  model  proposed  by  Wisher  et  al.  (1990)  or  the 
quantitative  models  described  by  Townsend  (1990)  has  been  fit  to  results. 

What  is  needed  now  is  an  evaluation  of  just  how  well  the  various  quantitative 
models  do  indeed  fit  the  data  (Objective  3) . 

A  MODEL  OF  MORSE  CODE  COPY 

The  models  of  Morse  intercept  that  we  developed  correspond  closely  to  the 
advanced  state  models  proposed  by  Wisher  et  al.  (1990) .  To  begin,  we  want  to 
tal)c  about  the  most  general  model  (GenCOPB)  that  we  developed,  one  which  al¬ 
lows  for  copy-behind  and  large-capacity  buffers  at  each  stage.  We  then  want 
to  describe  the  two  models  that  we  tested.  Both  assume  no  buffering  of  the 
input  at  any  stage.  The  first  allows  for  copy-behind  (LimCOPB) ;  the  second 
does  not  allow  for  copy-behind  (NoCOPB) . 

There  are  a  total  of  77  parameters  in  the  GenCOPB  model,  14  of  which  are 
free  and  63  of  which  are  fixed  by  the  data  and/or  condition.  There  are  a  to¬ 
tal  of  68  parameters  in  both  the  LimCOPB  and  NoCOPB  models,  5  of  which  are 
free  and,  again,  63  of  which  are  fixed  by  the  data  and/or  condition.  These 
parameters  are  described  in  detail  in  Appendix  A.  They  are  also  described  in 
context  in  the  material  below. 

We  assume  throughout  that  the  various  server  and  decay  times  are  inde¬ 
pendent  and  that  their  distribution  is  well  described  as  a  gamma.  Briefly, 
the  density  function  of  a  gamma  is  given  by  the  right  hand  side  of  the  equa¬ 
tion  below: 

f(x)  =  [rV‘^e"^''P]/r(o)  if  X  >  0, 
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0 


otherwise . 


The  parameter  a  is  referred  to  as  the  shape  parameter;  p  is  referred  to  as  the 
scale  parameter.  The  mean  and  variance  are  given  by,  respectively: 

E[X]  =  ap, 

VAR[X]  =  aP^. 

The  function,  r(a),  is  equal  to  a!  when  a  is  a  nonnegative  integer. 

General  Copy  Behind  Model  (GenCOPB) 

The  Stages.  We  begin  by  developing  a  very  general  model,  the  GenCOPB, 
which  allows  copy  behind.  Specifically,  we  assume  the  existence  of  three 
stages,  perception,  recognition  and  execution.  Each  stage  consists  of  a 
"server"  (some  process)  and  a  "buffer"  (some  memory) .  Specifically,  we  assume 
the  existence  of  three  servers.  We  assume  the  existence  of  three  buffers,  one 
in  front  of  each  of  the  servers.  And  we  assume  that  decay  can  occur  in  any 
one  of  the  three  buffers. 

We  assume  that  the  stages  are  in  series,  where  the  string  of  tones  input 
to  the  perception  stage  is  output  as  a  string  of  elements  (dits  and  dahs) . 

This  string  of  elements  is  input  to  the  recognition  stage  and  output  as  a 
character  code.  And  finally,  the  character  code  is  input  to  the  execution 
stage  and  output  as  an  actual  response. 

It  will  be  useful  to  have  a  more  formal  representation  for  the  various 
inputs.  To  begin,  note  that  there  are  31  different  possible  tone  inputs,  26 
letters  and  5  special  characters  (the  10  digits  were  not  included  when  fitting 
the  model  and  thus  are  not  discussed) .  Each  tone  input  t^  is  associated  with 

a  particular  element  string  e^.  Each  element  string  e^  is  associated  with  a 

particular  character  code  c^.  And  each  character  code  c^  is  associated  with  a 

particular  response  r^.  For  example,  the  tone  consisting  of  three  50  ms 

pulses  is  associated  with  the  element  string,  di-di-dit.  The  element  string 
di-di-dit  is  associated  with  the  code  for  the  character  s.  And  the  code  for 
the  character  s  is  associated  with  the  motor  program  for  pressing  the  letter  s 
on  the  keyboard. 
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Given  this  notation,  a  signal  is  output  correctly  if  tone  t^  is  the  input 
and  response  r^  is  the  output;  a  signal  is  output  as  an  incorrect  character  if 
tone  t^  is  input  and  response  r^  is  the  output  (j  =  1,,..,31,  i  /  j) .  A  sig¬ 
nal  is  output  as  a  period  if  tone  t^  is  the  input  and  a  period  is  the  output 

(a  period  is  indicated  by  the  integer  44) .  And  a  signal  is  output  as  no 
response  if  tone  t^^  is  the  input  and  no  response  is  made  to  this  tone. 

The  Servers .  Implicit  in  the  above  discussion  of  the  stages  is  the  fol¬ 
lowing  description  of  the  role  played  by  each  of  the  servers.  Specifically, 
as  in  the  model  proposed  by  Wisher  et  al.  (1990),  the  perception  server  maps 
sound  energy  into  a  string  of  Morse  elements,  i.e.,  dits  and  dahs.  The  recog¬ 
nition  server  maps  the  string  of  Morse  elements  into  a  character.  And  the 
execution  server  both  organizes  a  response  and  executes  it. 

We  made  several  simplifying  assumptions  about  the  distributions  of  the 
service  times.  To  begin,  we  assumed  that  the  duration  T^  of  the  recognition 

process  and  the  duration  T  of  the  execution  process  did  not  depend  on  the 

identity  of  the  input  to  the  process  or  on  the  output  of  the  process. 
Furthermore,  given  our  work  with  other,  similar  cognitive  tasks,  it  seemed 
reasonable  to  assimie  that  the  durations  of  the  recognition  and  execution 
processes  were  independent,  gamma  distributed  random  variables.  We  let  and 

represent  the  shape  and  scale  parameters  for  the  distribution  of  the  dura¬ 
tion  of  the  recognition  process;  and  we  let  and  p^  represent  the  shape  and 

scale  parameters  for  the  distribution  of  the  duration  of  the  execution 
process.  These  four  parameters  were  free. 

Although  we  did  not  assume  that  the  durations  of  the  recognition  and  ex¬ 
ecution  processes  depended  on  the  input  to  these  processes,  we  did  assiane  that 
the  duration  Tple^lt^)  of  the  perception  process  was  dependent  on  the  identity 

of  the  tone  sequence  t^  which  was  given  as  input.  At  this  point  we  do  not  as¬ 
sume  that  the  perception  time  depends  on  the  output  of  the  process.  However, 
we  believe  that  such  a  dependency  will  be  necessary  to  incorporate  in  future 


models  and  thus  we  keep  this  possibility  open  here  by  writing  Tpie^lt^)  in¬ 
stead  of  writing  more  simply,  Again,  we  assumed  that  the  durations  of 

the  perception  process  for  each  tone  sequence  were  independent,  gamma  dis¬ 
tributed  random  variables.  We  used  32  parameters  to  describe  these 
distributions,  31  shape  parameters  ap(i)/  i  =  1,...,31  and  one  common  scale 

parameter  p^.  These  parameters  were  fixed  by  the  data  in  the  latency  confu¬ 


sion  matrix  (see  Appendix  A  for  a  discussion  of  the  way  the  parameters  were 
set;  see  Appendix  B  for  a  discussion  of  the  latency  confusion  matrix) . 


The  Buffers .  We  made  several  critical  assumptions  about  the  buffers. 
Specifically,  we  assumed  that  when  a  stimulus  arrived  at  the  perception  buffer 
and  the  perception  buffer  was  empty,  it  went  immediately  to  the  perception 
server.  We  assumed  that  when  a  stimulus  arrived  at  the  perception  buffer  and 
the  buffer  had  space  available  for  it,  the  stimulus  queued  for  service  on  a 
first  come,  first  serve  basis.  And  we  assumed  that  when  a  stimulus  arrived  at 
the  perception  buffer  and  the  buffer  was  full,  the  stimulus  was  lost  from  the 
system.  This  generated  n^  response  from  the  model.  Similar  remarks  apply  to 
stimuli  arriving  at  the  recognition  and  execution  buffers. 


We  also  assumed  that  when  a  stimulus  decayed  from  any  buffer,  the 
stimulus  was  lost,  regardless  of  the  level  of  processing  which  had  been  com¬ 
pleted  immediately  prior  to  the  loss.  This  loss  also  generates  a  no  response 
in  the  model. 


A  total  of  six  parameters  were  needed  to  describe  the  decay  times  at  each 
of  the  three  buffers.  Specifically,  we  assumed  that  the  decay  times  were  in¬ 
dependent,  gamma  distributed  random  variables.  The  shape  and  scale  parameters 
of  the  distribution  of  the  decay  time  in  the  perception  buffer  are  noted  by, 
respectively,  a*  and  p*.  The  shape  and  scale  parameters  of  the  distribution 

of  the  decay  time  in  the  recognition  buffer  are  noted  by,  respectively,  o*  and 

P*.  And  the  shape  and  scale  parameters  of  the  distribution  of  the  decay  time 

in  the  execution  buffer  are  noted  by,  respectively,  a*  and  P|.  These  six 

parameters  are  free. 
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A  total  of  three  parameters  were  needed  to  describe  the  capacity  size  of 
the  three  buffers.  Specifically,  the  capacity  sizes  of  the  perception,  recog¬ 
nition  and  execution  buffers  are  noted  by,  respectively,  and  These 

three  parameters  are  free. 

Critical  Perception  Time.  We  made  one  critical  assumption  about  the  per¬ 
ception  server,  an  assumption  which  produced  a  period  response  from  the  model. 
Specifically,  we  assumed  that  if  the  time  to  perceive  the  stimulus  was  longer 
than  some  critical  value,  say  8,  then  the  model  coded  the  incoming  stimulus  as 
a  period.  The  rationale  for  doing  so  is  the  following.  As  the  time  it  takes 
an  operator  to  perceive  a  stimulus  increases,  the  data  suggest  that  accuracy 
decreases.  Thus,  rather  than  make  an  incorrect  response,  it  may  be  optimal 
for  the  operator  to  produce  a  period  as  a  response  (which,  during  training, 
earns  the  operator  a  smaller  point  loss  than  an  incorrect  response) . 

Note  that  we  assume  that  the  duration  of  the  perception  process  continues 
beyond  the  critical  value  8 .  It  is  only  after  the  perception  process  has  com¬ 
pleted  that  a  comparison  is  made  between  the  duration  of  this  process  and  the 
critical  value.  This  may  appear  counterintuitive  at  first  glance.  However, 
note  that  it  could  well  be  less  time  consuming  simply  to  perceive  a  stimulus 
and  then  determine  whether  the  critical  value  has  been  exceeded  than  concur¬ 
rently  both  to  perceive  the  incoming  stimulus  and  to  monitor  the  duration  of 
this  process.  The  exhaustive  scanning  mechanism  proposed  by  Sternberg  (1969) 
requires  for  its  justification  a  similar  argiraient.  The  critical  time  was  a 
free  parameter. 

Character  Errors .  Define  a  character  error  as  a  transformation  of  a 
character  (as  opposed  to  period)  input  at  one  stage  to  a  character  output  at 
the  same  stage  which  is  not  associated  with  the  input.  So,  a  character  error 
at  the  perception  stage  is  a  transformation  of  a  tone  sequence  t^  into  an  ele¬ 
ment  sequence  Cj,  where  i,j  =  1, ,31,  i  j.  A  character  error  at  the 

recognition  stage  is  a  transformation  of  an  element  string  e^  into  a  character 
code  Cj,  where  i,j  =  1,...,31,  i  ^  j.  And  a  character  error  at  the  execution 
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stage  is  a  transformation  of  a  character  code  c^  into  a  response  r^,  where 
i,  j  =  1, . .  .,31,  i  9^  j. 

We  need  to  consider  the  possibility  of  character  errors  at  all  three 
stages.  To  begin,  consider  character  errors  at  the  perception  stage.  When 
the  duration  of  the  perception  server  was  less  than  6 ,  on  some  trials  the  per¬ 
ception  stage  worked  correctly  (i.e.,  tone  sequence  t^  was  mapped  into  element 

string  e^^) ;  on  other  trials  it  produced  a  character  error.  Since  there  were 

31  tone  sequences,  we  need  to  specify  for  the  perception  stage  the  31  prob¬ 
abilities  p(e^lt^),  i  =  1,...,31.  These  31  probabilities  were  determined  by 

the  error  confusion  matrix  (see  Appendix  A  for  a  discussion  of  how  the 
parameters  were  estimated;  see  Appendix  B  for  a  discussion  of  the  error  confu¬ 
sion  matrix) .  To  keep  the  analysis  relatively  simple,  we  grouped  the  off- 
diagonal  confusions.  Thus,  the  probability  of  an  incorrect  perception 
transformation  was  set  equal  to  1  -  p(e^|t^). 

Next,  we  need  to  consider  the  possibility  of  character  errors  at  the 
recognition  and  execution  stages.  Implicit  in  the  above  is  the  assumption 
that  all  confusions  are  bundled  into  the  perception  stage.  That  is,  we  do  not 
also  allow  confusions  at  the  recognition  or  execution  stages.  This  is  not  be¬ 
cause  we  believe  that  such  confusions  do  not  occur.  Rather,  it  is  because  the 
preliminary  data  suggest  that  the  majority  of  confusions  occur  at  the  percep¬ 
tion  stage.  Of  course,  if  our  simplifying  assumption  is  wrong,  we  should  find 
the  GenCOPB  unable  to  account  for  significant  aspects  of  the  data. 

Limited  Copy  Behind  (LimCOPB)  and  No  Copy  Behind  (NoCOPB)  Models 

The  limited  copy  behind  (LimCOPB)  is  a  special  case  of  the  general  copy 
behind  (GenCOPB)  model.  Specifically,  in  the  LimCOPB  model,  we  assume  that 
the  buffers  are  of  size  zero  at  each  of  the  perception,  recognition,  and  ex¬ 
ecution  stages.  Note  that  this  means  we  do  not  need  parameters  to  describe 
the  duration  of  the  decay  in  these  three  buffers. 
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Strictly  speaking,  the  no  copy  behind  (NOCOPB)  model  is  not  a  special 
case  of  the  GenCOPB  model  since  in  the  NoCOPB  model  all  downstream  processors 
are  surveyed,  not  just  the  processor  at  the  next  stage.  Specifically,  if 
there  are  no  items  on  any  one  of  the  perception,  recognition  or  execution 
servers,  then  an  arriving  item  is  placed  on  the  perception  server.  Otherwise, 
the  arriving  item  is  lost  from  the  system  and  not  recoverable. 

SIMULATION 

We  now  want  to  describe  how  we  simulated  the  behavior  of  the  three 
models,  the  GenCOPB,  LimCOPB,  and  NoCOPB.  Recall  that  we  are  interested  in 
seven  dependent  variables,  four  related  to  accuracy  and  three  related  to  time. 
Specifically,  we  are  interested  in  simulating  the  probability  of  a  correct 
response,  the  probability  of  an  incorrect  character  response,  the  probability 
of  an  incorrect  period  response,  and  the  probability  of  no  response.  In  addi¬ 
tion,  we  are  interested  in  simulating  the  time  that  it  takes  to  make  a  correct 
response,  the  time  that  it  takes  to  make  a  incorrect  character  response,  and 
the  time  that  it  takes  to  make  an  incorrect  period  response. 

Because  the  general  copy-behind  model  was  complex,  we  decided  first  to 
simulate  its  behavior.  The  simulation  follows  directly  from  an  OP  repre¬ 
sentation  of  processing  in  such  a  task  (Fisher  and  Goldstein,  1983;  Goldstein 
and  Fisher,  1991) .  To  begin,  we  describe  the  state  variables  used  in  the 
simulation.  We  then  describe  the  actual  transition  rules.  Finally,  we 
describe  where  in  the  simulation  we  obtain  the  relevant  dependent  variables. 
(Note  that  we  will  wait  until  the  discussion  of  the  analytical  work  to 
describe  the  actual  OP  diagram.)  The  source  code  is  listed  in  Appendix  C. 

States 


Conceptually,  it  is  a  straightforward  matter  to  simulate  the  behavior  of 
the  system  we  have  described.  This  is  because  for  purposes  of  the  simulation 
we  can  classify  the  evolution  of  the  system  over  time  into  a  finite  number  of 
states.  Furthermore,  the  relevant  state  variables  themselves  are  constant 
throughout  the  duration  of  any  given  state.  Thus,  we  can  simulate  the  system 
as  a  series  of  discrete  events. 


12 


Each  state  is  defined  as  a  structure  consisting  of  a  (y^+  1)  X  5  ar¬ 
rival  matrix  M  (i) ,  a  (y  +  1)  X  5  perception  matrix  M  (i) ,  a  (y  +  1)  X  5 
a  p  p  r 

recognition  matrix  M^(i),  and  a  (y 1)  X  5  execution  matrix  M^(i) •  The  first 

row  of  the  arrival  (perception,  recognition,  execution)  matrix  contains  infor¬ 
mation  on  the  item  currently  on  the  arrival  (perception,  recognition, 
execution)  server.  The  second  and  subsequent  rows  contain  information  on  the 
items  in  the  arrival  (perception,  recognition,  execution)  buffer. 

As  noted  above,  the  31  characters  are  represented  by  the  integers  1-31 
(and  the  digits  by  the  numbers  32  -  41) .  The  intercharacter  interval  is  rep¬ 
resented  by  the  number  42.  The  intergroup  interval  is  represented  by  the 
number  43.  The  period  character  is  represented  by  the  number  44.  Note  that 
we  will  refer  to  the  generic  entity  which  passes  through  the  system  as  an 
item,  i.e.,  a  tone  sequence,  element  string,  character  code  and  response  are 
all  items. 

Arrival .  Consider  just  the  arrival  matrix  for  state  s^.  It  may  help 

some  to  display  a  particular  matrix  and  refer  to  the  entries  in  the  matrix  as 
we  discuss  them.  To  keep  things  simple,  we  will  display  only  two  rows  in  the 
arrival  matrix: 


(server  entries) 
(buffer  entries) 


To  begin,  consider  the  entries  in  the  first  row.  The  number  (23)  in  the  first 
column  and  first  row  indicates  the  identity  of  the  tone  sequence  in  state  s^ 

which  is  currently  on  the  arrival  server.  The  number  (23)  in  the  second 
column  and  first  row  indicates  the  identity  of  the  tone  sequence  which 
originally  was  sent.  In  the  arrival  matrix,  this  entry  and  the  preceding 
entry  are  identical  since  they  reference  the  same  item  before  any  processing 
has  begun.  The  entry  (1)  in  the  third  column  and  first  row  indicates  the 
serial  position  of  the  arriving  tone  sequence  in  the  group  in  which  it  was 
sent.  Since  there  are  five  tone  sequences  in  a  group,  the  serial  position 


23 

23 

1 

49 

0 

17 

17 

2 

ce 

0 
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will  always  lie  between  1  and  5,  inclusive.  The  number  (49)  in  the  fourth 
column  and  first  row  indicates  the  amount  of  time  whi::  remains  before  the 
tone  sequence  identified  in  columns  1  and  2  actually  arrives.  The  number  (0) 
in  the  fifth  column  and  first  row  indicates  the  amount  of  time  the  tone  se¬ 
quence  will  have  spent  in  the  system  when  it  finishes  processing  in  the 
arrival  stage.  This  time  is  zero  in  the  arrival  matrix  because  it  is  assumed 
that  timing  does  not  begin  until  the  arrival  process  is  complete. 

Consider  next  the  second  and  subsequent  rows  in  the  arrival  matrix.  The 
entries  (17  and  17)  in  the  first  and  second  columns  of  row  2  are  identical  and 
indicate,  respectively,  the  identity  of  the  next  tone  sequence  which  will  be 
sent  to  an  operator  and  the  identity  of  the  tone  sequence  which  was  originally 
sent  to  the  operator.  Note  that  the  next  tone  sequence  sent  to  the  operator 
is  the  first  tone  sequence  in  the  arrival  buffer.  The  number  (2)  in  the  third 
column  and  second  row  indicates  the  serial  position  in  the  group  in  which  it 
was  sent  of  the  tone  sequence  which  is  at  the  top  of  the  arrival  buffer.  If 
the  serial  position  of  the  arriving  tone  sequence  lies  between  1  and  4,  then 
the  serial  position  of  the  first  tone  sequence  in  the  buffer  is  increased  by 
1.  If  the  serial  position  of  the  arriving  tone  sequence  is  at  the  end  of  a 
group  (i.e.,  is  5),  then  the  serial  position  of  the  first  tone  sequence  in  the 
buffer  will  be  1.  The  serial  position  of  an  intercharacter  or  intergroup 
blank  is  set  at  0.  The  entry  («)  in  the  fourth  column  and  second  row  indic¬ 
ates  the  time  remaining  before  the  tone  sequence  which  is  at  the  top  of  the 
arrival  buffer  decays.  We  assume  that  this  time  is  infinite,  i.e.,  that  there 
is  no  decay  in  the  arrival  buffer.  Finally,  the  number  (0)  in  the  fifth 
column  and  second  row  indicates  the  total  time  that  the  tone  sequence  at  the 
top  of  the  arrival  buffer  will  have  spent  in  the  system  when  the  tone  sequence 
currently  being  processed  completes  service.  As  above,  this  is  zero  because 
it  is  assumed  that  timing  does  not  begin  until  the  arrival  is  complete. 

Perception.  The  entries  in  the  perception,  recognition  and  execution 
matrices  are  defined  in  much  the  same  fashion  as  the  above  entries.  To  begin, 
consider  the  entries  in  the  perception  matrix  and,  in  particular,  the  entries 
in  the  first  row.  Again,  it  may  help  to  display  an  actual  matrix,  one  where 
we  assume  that  the  buffer  can  hold  at  most  one  tone  sequence: 


(server  entries) 
(buffer  entries) 


16 

16 

2 

123 

650 

8 

8 

4 

500 

150 

The  number  (16)  in  the  first  column  and  first  row  indicates  the  identity  of 
the  tone  sequence  in  state  s^  which  is  currently  being  mapped  to  an  element 

string.  The  number  (16)  in  the  second  column  and  first  row  indicates  the 
identity  of  the  tone  sequence  which  originally  was  sent.  Since  the  arrival 
stage  always  outputs  the  same  tone  sequence  which  was  sent  as  input,  the  tone 
sequences  in  the  first  and  second  columns  are  identical  in  the  perception 
matrix  (as  well  as  the  arrival  matrix) ,  The  entry  (2)  in  the  third  column  and 
first  row  indicates  the  serial  position  of  the  tone  sequence  being  perceived 
in  the  group  in  which  it  was  sent.  The  nxmber  (123)  in  the  fourth  column  and 
first  row  indicates  the  amount  of  time  which  remains  before  the  tone  sequence 
in  the  process  of  being  perceived  actually  completes  this  process.  The  number 
(650)  in  the  fifth  column  and  first  row  indicates  the  amount  of  time  the  tone 
sequence  being  perceived  will  have  spent  in  the  system  when  it  finishes 
processing  in  the  perception  stage.  This  time  will  be  the  sum  of  the  times 
that  it  spent  in  the  perception  buffer  and  on  the  perception  server.  Thus, 
for  example,  if  tone  sequence  16  arrived  at  the  perception  stage  when  the 
server  was  occupied,  if  it  took  200  ms  to  complete  the  servicing  of  the  tone 
sequence  currently  on  the  server,  and  if  it  took  450  ms  to  service  tone  se¬ 
quence  16,  then  the  amount  of  time  tone  sequence  16  will  have  spent  in  the 
system  when  it  finishes  processing  in  the  perception  stage  will  be  650  ms. 

Consider  next  the  second  row  in  the  perception  matrix.  The  entry  (8)  in 
the  first  colinnn  of  row  2  indicates  the  identity  of  the  next  tone  sequence 
which  will  be  sent  to  the  perception  server  (i.e.,  the  identity  of  the  first 
item  in  the  perception  buffer) ,  The  entry  (8)  in  the  second  column  of  row  2 
indicates  the  identity  of  the  tone  sequence  in  the  first  row  in  the  perception 
buffer  which  originally  was  sent  to  the  arrival  stage.  Again,  this  is  the 
same  as  the  entry  in  the  first  column  and  second  row  since  we  assume  that  no 
errors  are  introduced  by  the  arrival  stage.  The  number  (4)  in  the  third 
column  and  second  row  indicates  the  serial  position  in  the  group  in  which  it 
was  sent  of  the  tone  sequence  which  is  at  the  top  of  the  perception  buffer. 
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Note  that  in  this  example  this  number  is  not  simply  one  more  than  the  cor¬ 
responding  number  in  the  first  row  and  third  column.  What  this  indicates  is 
that  the  third  tone  sequence  in  the  group  must  have  entered  the  perception 
stage  at  a  time  when  the  buffer  was  full.  Thus,  it  was  lost  from  the  system. 
The  number  (500)  in  the  fourth  column  and  second  row  indicates  the  time 
remaining  before  the  tone  sequence  which  is  at  the  top  of  the  perception  buff¬ 
er  decays.  Finally,  the  number  (150)  in  the  fifth  column  and  second  row 
indicates  the  total  time  that  the  tone  sequence  at  the  top  of  the  perception 
buffer  will  have  spent  in  the  system  when  the  tone  sequence  currently  on  the 
perception  processor  completes  its  service.  In  this  case,  tone  sequence  8 
could  have  arrived  when  there  was  150  ms  left  to  process  of  tone  sequence  16 
(of  course,  other  scenarios  could  have  led  to  the  same  time)  .  Similar  remarks 
would  apply  to  subsequent  rows  in  the  perception  matrix  were  we  to  assume  a 
larger  buffer. 

Recognition.  Next,  consider  the  entries  in  the  recognition  matrix  and, 
in  particular,  consider  the  entries  in  the  first  row  of  the  recognition 
matrix.  Assume  that  the  buffer  can  hold  only  one  element  string: 


Mp(i): 


10 

12 

5 

321 

1250 

44 

14 

1 

400 

621 

(server  entries) 
(buffer  entries) 


The  number  (10)  in  the  first  column  and  first  row  indicates  the  identity  of 
the  element  string  in  state  s^  which  is  currently  being  recognized.  The  num¬ 
ber  (12)  in  the  second  column  and  first  row  indicates  the  identity  of  the  tone 
sequence  which  originally  was  sent  and  which  has  been  transformed  into  the 
element  string  now  on  the  recognition  server.  Since  the  perception  stage  can 
introduce  errors,  these  two  numbers  need  not  be  identical  (and  are  not  in  this 
example) .  The  entry  (5)  in  the  third  column  and  first  row  indicates  the 
serial  position  of  the  element  string  which  is  on  the  recognition  server  in 
the  group  in  which  it  was  sent.  The  number  (321)  in  the  fourth  column  and 
first  row  indicates  the  amount  of  time  which  remains  before  the  element  string 
on  the  recognition  server  completes  processing  on  this  server.  The  number 
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(1250)  in  the  fifth  cclumn  and  first  row  indicates  the  amount  of  time  the  ele¬ 
ment  string  and  associated  tone  sequence  will  have  spent  in  the  system  when 
the  element  string  on  the  recognition  server  finishes  processing.  This  time 
will  be  the  sum  of  the  times  that  the  item  spent  in  the  perception  and  recog¬ 
nition  buffers  plus  the  sum  of  the  times  it  spent  on  the  perception  and 
recognition  servers. 

Consider  next  the  second  row  in  the  recognition  matrix.  The  entry  (44) 
in  the  first  column  of  row  2  indicates  the  identity  of  the  next  element  string 
which  will  be  sent  to  the  recognition  server  (i.e.,  the  identity  of  the  first 
element  string  in  the  recognition  buffer)  .  In  this  particular  exairiple,  the 
entry  represents  a  period  response,  an  indication  that  the  critical  time  on 
the  perception  server  was  exceeded  for  the  corresponding  tone  sequence  and 
thus  the  perception  server,  rather  than  generating  an  element  string  as¬ 
sociated  with  a  character  output,  generated  an  element  string  for  a  period  as 
an  output.  The  entry  (14)  in  the  second  column  of  row  2  indicates  the  iden¬ 
tity  of  the  tone  sequence  originally  sent  which  is  now  the  first  element 
string  in  the  recognition  buffer.  The  number  (1)  i:,  the  third  column  and 
second  row  indicates  the  serial  position  in  the  group  in  which  it  was  sent  of 
the  element  string  which  is  at  the  top  of  the  recognition  buffer.  The  number 
(400)  in  the  fourth  column  and  second  row  indicates  the  time  remaining  before 
the  element  string  which  is  at  the  top  of  the  recognition  buffer  decays. 
Finally,  the  number  (621)  in  the  fifth  column  and  second  row  indicates  the  to¬ 
tal  time  that  the  element  string  at  the  top  of  the  recognition  buffer  (and 
associated  tone  sequence)  will  have  spent  in  the  system  when  the  element 
string  current  on  the  recognition  processor  completes  its  service.  Similar 
remarJcs  apply  to  subsequent  rows  in  the  recognition  matrix  were  the  buffer  to 
hold  more  than  one  element  string. 

Execution.  Finally,  consider  the  entries  in  the  execution  matrix  and,  in 
particular,  consider  the  entries  in  the  first  row  of  the  execution  matrix. 
Again,  it  may  be  helpful  to  describe  these  entries  in  the  context  of  a  par¬ 
ticular  example.  Assume  that  the  execution  buffer  can  hold  only  one  character 
code: 
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20 

10 

3 

101 

1600 

18 

18 

4 

20 

2210 

(server  entries) 
(buffer  entries) 


The  number  (20)  in  the  first  column  and  first  row  indicates  the  identity  of 
the  character  code  in  state  s^^  for  which  a  response  is  currently  being  or¬ 
ganized.  The  number  (10)  in  the  second  column  and  first  row  indicates  the 
identity  of  the  tone  sequence  which  originally  was  sent  and  which  has  been 
transformed  into  element  string  20  at  the  perception  stage  and  character  code 
20  at  the  recognition  stage.  The  entry  (3)  in  the  third  column  and  first  row 
indicates  the  serial  position  of  the  character  code  which  is  on  the  execution 
server  in  the  group  in  which  it  was  sent.  The  number  (101)  in  the  fourth 
column  and  first  row  indicates  the  amount  of  time  which  remains  before  the 
character  code  on  the  execution  server  completes  processing  on  this  server. 

The  number  (1600)  in  the  fifth  column  and  first  row  indicates  the  amount  of 
time  the  character  code  and  its  associated  element  string  and  tone  sequence 
will  have  spent  in  the  system  when  character  code  finishes  processing  in  the 
execution  stage.  This  time  will  be  the  sum  of  the  times  that  the  item  spent 
in  the  perception,  recognition  and  execution  buffers  plus  the  sum  of  the  times 
it  spent  on  the  perception,  recognition  and  execution  servers. 

Consider  next  the  second  row  in  the  execution  matrix.  The  entry  (18)  in 
the  first  column  of  row  2  indicates  the  identity  of  the  next  character  code 
which  will  be  sent  to  the  execution  server  (i.e.,  the  identity  of  the  first 
character  code  in  the  execution  buffer) .  The  entry  (18)  in  the  second  column 
of  row  2  indicates  the  identity  of  the  tone  sequence  which  through  processing 
in  the  perception  and  recognition  stages  has  been  transformed  into  the  charac¬ 
ter  code  at  the  top  of  the  execution  buffer.  The  number  (4)  in  the  third 
column  and  second  row  indicates  the  serial  position  in  the  group  in  which  it 
was  sent  of  the  character  code  which  is  at  the  top  of  the  execution  buffer. 

The  number  (20)  in  the  fourth  column  and  second  row  indicates  the  time  remain¬ 
ing  before  the  character  code  which  is  at  the  top  of  the  execution  buffer 
decays.  Finally,  the  number  (2210)  in  the  fifth  column  and  second  row  indic¬ 
ates  the  total  time  that  the  character  code  at  the  top  of  the  execution  buffer 
(and  its  associated  element  string  and  tone  sequence)  will  have  spent  in  the 
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system  when  the  stimulus  current  on  the  execution  processor  completes  its 
service . 

Transition  Rules 

The  simulation  of  the  system  requires  knowledge  of  the  state  variables, 
as  explained  above.  It  also  requires  knowledge  both  of  the  events  which  lead 
to  a  transition  between  states  and  of  the  rules  used  to  relate  states  which 
follow  one  another.  We  now  want  to  describe  these  events  and  rules. 

Briefly,  a  transition  between  states  occurs  when  an  item  (i.e.,  a  tone 
sequence,  element  string  or  character  code)  completes  service  or  an  item 
decays.  The  determination  of  which  item  will  complete  or  decay  first  can  be 
made  in  a  straightforward  fashion  from  the  entries  in  the  fourth  column  of  the 
arrival,  perception,  recognition  and  execution  matrices.  Specifically,  the 
minimum  time  is  selected  from  the  set  of  times  consisting  both  of  the  arrival, 
perception,  recognition  and  execution  service  times  and  of  the  perception, 
recognition  and  execution  decay  times. 

Arrival .  We  now  need  to  specify  exactly  what  changes  are  made  when 
anyone  of  the  above  events  occur.  Although  quite  tedious,  the  details  vary 
enough  for  the  different  stages  to  require  a  relatively  complete  rendering  of 
the  transition  rules.  To  begin,  consider  the  arrival  stage.  The  item  in  the 
arrival  stage  which  completes  first  can  be  either  a  character  or  a  blank.  We 
need  to  consider  both.  And  the  item  which  completes  first  in  the  arrival 
stage  can  find  the  perception  buffer  full  or  not  full.  Again,  we  need  to  con¬ 
sider  both  possibilities. 

i)  Assume  that  the  item  which  completes  first  is  a  character  in  the  ar¬ 
rival  stage.  An  example  of  a  state,  say  s^,  where  the  first  item  to  complete 

is  a  tone  sequence  in  the  arrival  stage  is  given  below  on  the  left;  the 
entries  in  the  new  state,  say  s^,  which  follows  it  are  given  on  the  right. 

For  the  sake  of  simplicity,  we  assume  that  the  buffers  in  the  perception, 
recognition  and  execution  stages  can  hold  only  one  item.  And  we  represent 
only  the  top  most  item  in  the  arrival  buffer: 
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Arrival  Stage  (s^) 


23 

23 

1 

49 

0 

26 

26 

2 

00 

0 

Percept 

:ion  Stc 

ige  (s^l 

19 

19 

5 

305 

600 

- 

- 

- 

I 

lecognit 

;ion  Stc 

ige  (s^) 

14 

14 

3 

200 

953 

11 

11 

4 

780 

876 

Execut 

-ion  Stc 

ige  (s^) 

8 

8 

1 

510 

1451 

9 

9 

2 

452 

1008 

Arrival  Stage  (Sj) 


26* 

26* 

2* 

250* 

0 

22* 

22* 

3* 

00 

0 

Perc 

:eption 

stage 

19 

19 

5 

256* 

600 

23* 

23* 

1* 

1200* 

256* 

Recc 

>gnitioi 

i  Stage 

14 

14 

3 

151* 

953 

11 

11 

4 

731* 

876 

Exec 

:ution  i 

Stage  (: 

8 

6 

1 

461* 

1451 

9 

9 

2 

403* 

1008 

As  noted  above,  in  order  to  determine  which  event  completes  first,  we  need  to 
look  at  the  entries  in  column  four  of  the  matrices  defining  state  s^  and  find 

the  minimum.  Doing  such,  we  see  that  the  times  remaining  on  the  arrival,  per¬ 
ception,  recognition  and  execution  servers  are  respectively  49,  305,  200  and 
510  ms.  And  we  see  that  the  times  remaining  before  decay  occurs  in  the  ar¬ 
rival,  recognition  and  execution  buffers  are,  respectively,  «,  780  and  452  ms 
(the  perception  buffer  is  empty,  as  indicated  by  the  "-")  .  Thus,  the  minimum 
time  to  completion  or  decay  is  49  ms,  or  the  remaining  arrival  time  of  tone 
sequence  23. 

The  changes  in  the  arrival,  perception,  recognition  and  execution 
matrices  are  given  on  the  right  in  the  example  state  s^.  The  starred  (*) 

entries  indicate  where  changes  occur.  To  begin,  consider  the  change  in  the 
arrival  matrix.  Since  tone  sequence  23  completes,  we  move  it  off  the  arrival 
server.  We  move  the  first  tone  sequence  in  the  arrival  buffer,  tone  sequence 
26,  to  the  arrival  server  and  assign  it  the  appropriate  arrival  time  (say 
250) .  Note  that  we  have  displayed  only  the  first  two  rows  of  the  arrival 
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matrix.  However,  we  are  assuming  that  the  arrival  buffer  contains  additional 
items,  and  in  particular,  contains  tone  sequence  22  below  tone  sequence  26. 
Thus,  we  move  tone  sequence  22  up  to  the  first  position  in  the  arrival  buffer. 
By  implication  it  is  the  third  tone  sequence  in  the  current  group. 

We  also  need  to  change  the  perception  matrix.  Note  that  although  the 
server  of  the  perception  matrix  is  occupied,  the  buffer  of  the  perception 
matrix  is  empty.  Thus,  we  can  put  the  tone  sequence  which  just  arrived  (23) 
in  the  buffer  of  the  perception  matrix  and  assign  it  a  decay  time,  say  1200. 

We  need  to  make  two  additional  changes.  First,  the  time  remaining  to  process 
the  tone  sequence  on  the  perception  server  must  be  reduced  by  49  ms,  from  305 
to  256  ms.  Second,  the  amount  of  time  the  newly  arrived  tone  sequence  23  will 
have  spent  in  the  perception  buffer  after  the  tone  sequence  currently  on  the 
server  completes  its  execution  must  be  assigned.  In  this  case,  this  is  simply 
256  ms.  Note  that  the  amount  of  time  that  tone  sequence  19  will  have  spent  in 
the  system  after  it  completes  processing  is  unaffected  by  the  completion  of 
the  tone  sequence  in  the  arrival  stage.  Thus,  the  entry  in  the  perception 
matrix  in  column  5  and  row  1  in  both  states  s^  and  s^  is  600. 

Finally,  we  need  to  decrement  the  server  and  decay  times  in  the  recogni¬ 
tion  and  execution  matrices.  For  example,  the  time  remaining  until  element 
string  11  decays  in  the  recognition  buffer  is  780  ms  when  we  first  examine  the 
system.  After  tone  sequence  23  arrives,  49  ms  will  have  elapsed.  Thus  the 
time  remaining  until  element  string  11  decays  is  reduced  from  780  to  731  ms. 

ii)  Above,  a  character  in  the  arrival  stage  was  the  first  item  to  com¬ 
plete.  It  could  have  been  a  blank  (an  intercharacter  or  intergroup  blank) 
instead.  An  example  is  given  below  (recall  that  42  is  the  code  for  an  inter¬ 
character  blank) : 
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Perception  Stage  (s.) 


Perception  Stage  (s.) 


Recognition  Stage  (s. 


Recognition  Stage  (s.) 


14 

14 

11 

11 

Execution  Stage  (s.) 


1 

4 

14 

3 

J 

953 

11 

11 

4 

731* 

876 

Execution  Stage  (s.) 


OO 

OO 

9 

9 

8 

8 

1 

461* 

1451 

9 

9 

2 

403* 

1008 

In  this  example,  note  that  all  the  times  are  changed  as  above  and  that  a  new 
tone  sequence  is  added  to  the  arrival  server.  However,  the  blank  is  not  moved 
to  the  perception  buffer  since  blanks  are  not  processed  through  the  system. 

iii)  Finally,  suppose  that  the  item  which  completes  first  is  on  the  ar¬ 
rival  server  and  that  the  perception  buffer  is  full.  Then,  we  make  the 
changes  indicated  below: 


23 

23 

26 

26 

Arrival  Stage  (s.) 


49 


Perception  Stage  (s.) 


Arrival  Stage  (s^) 


26* 

26* 

2 

22* 

22* 

3 

250* 

0 

eo 

0 

Perception  Stage  (s.) 


19 

19 

4 

305 

600 

16 

16 

5 

838 

403 

Recognition  Stage  (s.) 


19 

19 

16 

16 

256* 

600 

787* 

403 

Recognition  Stage  (s.) 


14 

14 

2 

200 
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11 

11 

3 

780 

876 

14 

14 

2 

151* 

953 

11 

11 

3 

731* 

876 

Execution  Stage  (s^) 


Execution  Stage  (s^) 

\ _ I _ 1_ _ _ _ _ _ L _ -  -J _ I  I  I  »  I  1  I 


8 

8 

5 

510 

1451 

9 

9 

1 

452 

1008 

Note  that  the  only  difference  here  is  that  tone  sequence  23  is  not  added  to 
the  perception  buffer  since  this  buffer  is  already  full. 

Perception.  Next,  we  need  to  consider  the  perception  stage.  An  event 
(i.e.,  a  completion  or  decay)  in  this  stage  will  drive  the  transition  if 
either  the  item  on  the  perception  server  completes  or  an  item  in  the  percep¬ 
tion  buffer  decays  before  any  other  service  or  decay  occurs.  To  begin,  assume 
that  the  event  which  leads  to  the  transition  is  the  completion  of  the  process¬ 
ing  of  the  tone  sequence  on  the  perception  server.  Then,  if  the  duration  of 
the  service  is  greater  than  some  critical  time,  we  assume  that  the  tone  se¬ 
quence  on  the  server  is  replaced  by  a  period  response.  For  illustrative 
purposes,  let  us  assume  that  this  time  is  1600  ms.  If  the  duration  of  the 
perception  process  is  less  than  1600  ms,  then  the  tone  sequence  t^  on  the  per¬ 
ception  server  is  replaced  with  its  correct  element  string,  e^,  with 

probability  p(r.|t.)  and  with  an  incorrect  element  string  e.  (i  j)  with 
11  3 

probability  1  -  p(r^ltj^).  Furthermore,  it  will  be  the  case  that  sometimes  the 

tone  sequence  on  the  perception  server  completes  (regardless  of  its  duration) 
when  the  recognition  buffer  is  empty;  at  other  times  it  will  complete  when  the 
recognition  buffer  is  full. 

As  noted  above,  it  may  be  the  decay  of  a  tone  sequence  in  the  perception 
buffer  (rather  than  the  processing  of  a  tone  sequence  on  the  perception  serv¬ 
er)  that  drives  the  next  transition.  If  this  is  the  case,  then  we  need  to 
enforce  alternative  transition  rules. 

The  full  set  of  rules  is  discussed  below. 

i)  To  begin,  assume  that  it  is  the  perception  service  time  which  is  the 
minimum  of  the  times  remaining  until  an  item  either  completes  service  or 


decays  from  a  buffer.  And  assume  the  perception  service  time  is  greater  than 
the  critical  time  of  1600  ms.  Then,  the  matrices  are  changed  as  in  the  fol¬ 
lowing  example: 


There  are  several  things  to  note.  First,  the  tone  sequence  19  on  the  percep¬ 
tion  server  now  becomes  the  element  string  44  in  the  recognition  buffer  (the 
string  associated  with  a  period  response) .  Second,  the  actual  service  time  in 
the  perception  stage  is  greater  than  1600  ms  and,  in  this  case,  is  equal  to 
1888  ms,  assuming  that  no  time  was  spent  in  the  perception  buffer. 
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ii)  Assume  that  the  first  item  to  complete  or  decay  is  in  the  perception 
stage,  as  above.  However,  now  assume  that  the  perception  service  time  is  less 
than  the  critical  perception  time  of  1600  ms.  Then,  as  noted  above,  one  of 
two  things  can  happen.  The  tone  secjuence  19  can  be  transformed  correctly  into 
the  element  string  19  or  it  can  be  transformed  into  some  other  element  string. 
When  the  element  string  was  transformed  incorrectly,  we  mapped  tone  sequence  1 


to  tone  sequence  2  and  tone  sequences  2, _ ,31  to  tone  sequence  1.  For  ex¬ 

ample,  below  tone  sequence  19  is  mapped  to  tone  sequence  1: 
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Note  that  the  time  tone  sequence  19  will  have  been  in  the  system  when  it  com¬ 
pletes  service  on  the  perception  server  is  only  742  ms,  as  opposed  to  1888  ms 
in  the  previous  example.  Thus,  a  period  response  is  not  generated. 

iii)  In  both  of  the  above  examples,  there  was  room  in  the  recognition 
buffer  for  the  output  of  the  perception  server.  This  will  not  always  be  the 
case.  When  this  happens,  the  times  are  updated  appropriately  and  the  element 
string  output  from  the  perception  server  is  bumped  from  the  system,  increment¬ 
ing  by  one  the  number  of  no  responses  in  the  simulation.  We  have  already 
described  how  this  change  is  made  when  we  discussed  the  transition  rules  for 
the  arrival  stage. 
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iv)  Finally,  suppose  that  it  is  the  decay  of  the  tone  sequence  in  the 
perception  buffer  that  drives  the  transition.  Then,  the  state  entries  might 

appears  as  below; 
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Note  that  tone  sequence  16  in  the  perception  buffer  has  decayed  and  thus  ap¬ 
pears  in  state  s^,  but  not  in  state  s^. 

Recognition  and  Execution.  The  transition  rules  in  the  recognition  and 
execution  stage  are  similar  to  those  in  the  perception  stage  with  one  excep¬ 
tion.  In  particular,  the  transformation  of  an  element  string  into  its 
corresponding  character  code  is  considered  errorless  and  the  transformation  of 
each  character  code  into  its  corresponding  motor  program  is  considered  error¬ 
less. 
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Dependent  Variables 

Information  on  six  of  the  seven  dependent  variables  we  obtain  from  the 
execution  matrix  at  each  point  in  time  when  the  first  item  in  a  state  to  com¬ 
plete  is  on  the  execution  server.  An  example  execution  matrix  is  listed 
below: 


Execution  Stage 

(server  entries) 
(buffer  1  entries) 
(buffer  2  entries) 
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Now,  consider  how  we  obtain  six  of  the  seven  measures  from  the  above 
matrix.  Note  that  1451  is  the  final  measure  of  the  response  time  for  tone  se¬ 
quence  8.  This  will  be  a  sample  of  the  correct  response  time  (since  the 
entries  in  coliunns  1  and  2  of  row  1  are  identical) .  Thus  we  increment  by  one 
the  number  of  correct  responses,  giving  us  the  information  we  need  eventually 
to  compute  the  probability  of  a  correct  response.  When  we  return  at  some 
point  in  the  future  to  this  matrix,  as  long  as  the  first  character  code  in  the 
buffer  has  not  decayed,  a  period  (44)  will  be  output.  In  the  fifth  column  of 
the  first  row  will  now  appear  the  time  it  took  to  make  the  period  response 
when  tone  sequence  9  was  the  original  item.  Thus,  we  have  a  sample  of  a 
period  response  time.  Furthermore,  we  increment  by  one  the  number  of  period 
responses,  giving  us  the  information  we  need  eventually  to  compute  the  prob¬ 
ability  of  a  period  response.  Finally,  when  we  return  at  a  still  later  point 
to  this  matrix,  as  long  as  the  second  character  code  in  the  buffer  has  not 
decayed,  a  1  will  be  output.  In  the  fifth  column  of  the  first  row  will  appear 
the  incorrect  response  time.  And  we  increment  by  one  the  number  of  incorrect 
response,  giving  us  the  information  we  need  eventually  to  compute  the  prob¬ 
ability  of  ^  incorrect  response. 

We  obtain  information  on  the  seventh  dependent  variable  as  the  system 
evolves.  That  is,  rather  than  keeping  track  of  decayed  items  or  items  which 
have  been  bumped  from  the  system  because  they  were  output  at  a  time  when  the 
immediate  downstream  buffer  was  full,  we  simply  increment  by  one  the  number  of 
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items  to  which  an  operator  will  make  no  response  at  the  time  an  item  disap 
pears.  Thus,  we  obtain  the  information  we  need  in  order  to  compute  the 
probability  of  no  response . 


ANALYSIS 


Although  the  simulation  as  described  above  yields  the  desired  predic¬ 
tions,  it  has  one  weakness  (a  weakness  common  to  all  simulations) . 
Specifically,  it  is  never  truly  possible  to  identify  the  parameter  settings 
which  maximize  the  fit  of  a  given  model  to  the  data.  The  best  that  can  be 
done  with  a  simulation  is  to  make  the  search  through  the  parameter  space  a 
relatively  fine  one,  hoping  that  there  exist  no  unexplored  points  which  would 
radically  alter  the  fit  of  the  particular  model  being  simulated.  Exact  iden¬ 
tification  of  the  parameter  settings  which  maximize  the  fit  can  only  be  done 
when  closed  form  expressions  can  be  obtained  for  the  dependent  variables  of 
interest. 

Unfortunately,  the  GenCOPB,  LimCOPB  and  the  NoCOPB  are  themselves  ex¬ 
tremely  complex  models,  making  it  difficult  if  not  impossible  to  obtain  the 
desired  closed  form  expressions.  And  even  were  one  to  obtain  the  expressions, 
they  themselves  might  be  too  unwieldy  to  work  with.  The  question  at  this 
point  is  therefore  a  threefold  one.  First,  in  principle  can  closed  form  ex¬ 
pressions  be  obtained  for  the  GenCOPB,  LimCOPB  and  NoCOPB  models?  Second,  if 
in  principle  such  expressions  can  be  obtained,  are  there  some  experimental 
paradigms  which,  unlike  the  current  paradigm,  introduce  enough  constraints  to 
make  it  possible  actually  to  derive  these  expressions  in  a  reasonable  amount 
of  time?  And  third,  how  would  one  obtain  these  expressions  were  such 
paradigms  to  exist? 

OP  Diagrams 

We  can  answer  the  first  question  quite  easily.  Indeed,  it  is  possible  to 
obtain  closed  form  expressions  for  the  various  expressions  in  which  we  have  an 
interest.  In  order  to  prove  that  this  is  the  case,  we  first  need  to  describe 
the  general  structure  of  an  OP  diagram.  Then,  we  need  to  show  that  the  state 
space  of  the  models  we  have  developed  is  itself  an  OP  diagram.  Once  we  know 
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the  structure  of  our  models  is  an  OP  diagram,  we  can  use  the  work  in  Fisher 
and  Goldstein  (1983;  Goldstein  and  Fisher,  1991)  to  obtain  the  desired  expres¬ 
sions  given  certain  auxiliary  assumptions. 

General  Structure .  Formally,  an  OP  diagram  is  a  structure  <r,'F,Q>  con¬ 
sisting  of  a  set  r  of  processes,  set  4'  of  paths,  and  set  Q  of  states.  The  set 
r  consists  of  n(r)  processes,  labelled  x^^, . . . ,x^ .  Each  state  in  fl  cor¬ 
responds  to  a  partition  of  the  processes  in  F.  The  states  are  labelled 
®1' '  ‘  ■'®n(Q)  ■  that  not  all  partitions  of  F  need  be  in  li.  Each  path  a  in 

*F  consists  of  a  sequence  of  states.  We  identify  the  states  on  path  a,  in  or¬ 
der  ,ass  f  s  ,...,s  . 

1  2  n(a) 

Several  remarks  need  to  be  made  about  the  states.  First,  the  processes 
in  each  state  s^  must  be  resident  in  one  and  only  one  of  four  sets,  the  pre¬ 
active  set  A(s^),  the  current  set  C(s^),  the  completed  set  K(s^),  or  the  other 
set  O(s^).  Second,  the  first  state  along  any  path  must  have  empty  completed 

and  other  sets  and  at  least  one  process  in  the  current  set.  Third,  the  last 
state  along  any  path  must  have  empty  pre-active  and  current  sets  and  at  least 
one  process  in  the  completed  set. 

There  are  restrictions  not  only  on  the  character  of  the  sets  of  states  at 
the  beginning  and  end  of  each  path,  but  also  restrictions  on  the  character  of 
states  adjacent  to  one  another  on  a  path.  These  restrictions  are  displayed 
graphically  in  Figure  1  below: 


Figure  1 

Current  Completed 

A(s^)  — >  C(s.)  A(s.)  -/->  K(s.) 


P re -Active 
A(s^) - >  A(s.) 

C(s.)  -/->  A(s.) 

K(s^)  -/->  A(.?.) 

0(s.)  -/->  A(s.) 


C(s^) - >  C(s.) 

K(s^)  -/->  C(s.) 
0(s.) - >  C(s.) 


C(s^) - >  K(Sj) 

K(s^)  - >  K(Sj) 

0{s.)  -/->  K(s,) 


Other 

A(s.)  -/->  0(Sj) 

C(s.) - >  0(s.) 

K(s.)  -/->  0(Sj) 
0(s.) - >  0(s.) 
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Note  that  an  arrow,  — >,  joining  two  different  sets  in  adjacent  states,  say 
A(s^)  — >  A(Sj),  means  that  later  set  (in  this  case,  the  pre-active  set  A(Sj) 

can  contain  processes  in  the  earlier  set  A(s^) .  The  negated  arrow,  -r->,  con¬ 
nected  two  different  sets  in  adjacent  states,  say  C(s^)  -/->  A(Sj),  means  that 
the  later  set,  A(Sj),  does  not  contain  any  elements  of  the  earlier  set,  C(s^). 


In  words,  if  state  Sj  is  an  immediate  successor  of  state  s^  for  some  path 
in  r,  then  restrictions  apply  to  the  four  sets  which  define  state  Sj.  First, 
the  pre-active  set  A(Sj)  of  state  s^  can  contain  no  processes  which  are  not 
contained  in  the  pre-active  set  A(s^)  of  state  s^.  Second,  the  current  set 
C(Sj)  of  state  s^  can  contain  no  processes  in  the  completed  set  K(s^)  of  state 
s^.  Third,  the  completed  set  K(Sj)  of  state  s^  must  contain  every  process  in 
the  completed  set  K(s^)  of  state  s^  as  well  as  at  least  one  process  from  the 
current  set  C(S^)  of  state  s^.  Furthermore,  no  processes  from  the  pre-active 
set  A(s^)  or  other  set  O(s^)  of  state  s^^  can  appear  in  the  completed  set  K(Sj) 
of  state  Sj.  Finally,  the  other  set  0(Sj)  of  state  s^  can  contain  no 
processes  in  the  pre-active  set  A(s^)  or  completed  set  K(s^)  of  state  s^. 


If  a  network  is  an  OP  diagram,  then  using  the  techniques  of  Fisher  and 
Goldstein  (1983;  Goldstein  and  Fisher),  we  can  obtain  expressions  for  the  be¬ 
havior  of  the  network  if  two  additional  conditions  are  met.  First,  the  number 

of  states  along  any  path  in  the  OP  diagram  must  be  finite.  And  second,  any 

process  which  becomes  current  in  one  or  more  states  along  a  path  must  even¬ 
tually  enter  a  completed  set. 

We  can  easily  show  that  the  various  models  we  have  described  can  be  rep¬ 
resented  as  OP  diagrams.  Briefly,  in  our  particular  case,  we  associate  seven 

processes  with  each  tone  sequence,  four  service  processes  (one  each  for  the 

arrival,  perception,  recognition  and  execution  servers)  and  three  decay 
processes  (one  each  for  the  perception,  recognition  and  execution  servers) . 

At  the  outset,  the  pre-active  set  consists  of  each  of  the  seven  processes  for 
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all  of  the  different  tone  sequences  except  for  the  arrival  service  process  for 
the  first  tone  sequence.  This  will  appear  in  the  current  set  of  the  start 
state.  When  the  arrival  service  process  completes,  this  process  enters  the 
completed  set  of  the  second  state  and  both  the  perception  service  process  as¬ 
sociated  with  the  first  tone  sequence  and  the  arrival  service  process 
associated  with  the  second  tone  sequence  enter  the  current  set.  And  .‘jo  on. 

Analysis 

In  theory,  since  the  various  models  can  be  represented  as  OP  diagrams  we 
can  derive  the  moments  of  the  response  time  analytically.  However,  this  does 
not  mean  we  can  obtain  closed  form  expressions  for  these  models,  or  even  if  we 
could  obtain  such  expressions,  that  we  could  do  so  in  a  reasonable  length  of 
time.  In  fact,  although  we  cannot  obtain  closed  form  expressions  as  the 
models  are  currently  set  forth,  we  can  easily  modify  the  model  so  that  such 
expressions  can  be  obtained  in  theory.  The  problems  requiring  modification 
come  at  two  points. 

First,  recall  that  in  the  perception  stage  we  took  one  path  if  the  per¬ 
ception  service  time  were  greater  than  the  critical  perception  time,  another 
path  if  this  were  not  the  case.  Thus,  knowledge  of  the  last  finishing  process 
and  the  sequence  of  states  up  until  the  current  state  does  not  alone  determine 
the  path  taken  out  of  the  current  state.  Strictly  speaking  then,  the  various 
models  we  have  proposed  cannot  be  represented  as  OP  diagrams  (Goldstein  and 
Fisher,  1991) ,  a  seeming  contradiction  to  our  conclusion  in  the  previous  sec¬ 
tion. 


Fortunately,  we  can  easily  get  around  the  above  problem.  Specifically, 
we  can  represent  the  critical  time  itself  as  a  process.  The  critical  time 
process  is  then  in  a  race  with  the  perception  process.  If  the  critical  time 
process  finishes  before  the  perception  process,  then  the  tone  sequence  on  the 
perception  server  is  replaced  with  a  period  when  the  server  completes.  If  the 
perception  process  finishes  before  the  critical  time  process,  ^-hen  the  tone 
sequence  on  the  perception  server  is  replaced  with  its  associated  element 
string  when  the  server  completes. 
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Second,  it  is  not  possible  to  obtain  closed  form  expressions  for  gamma 
distributed  random  variables  unless  the  gamma  is  a  sum  of  independent  exponen¬ 
tials.  This  does  not  overly  restrict  the  shape  of  the  distributions  and  seems 
a  reasonable  enough  assumption.  Thus,  below  we  assume  that  the  distributions 
of  the  various  processes  are  independent,  gamma  distributed  random  variables 
subject  to  the  constraint  that  the  gamma  random  variables  are  themselves  sums 
of  independent  exponentials. 

Feasibility.  As  noted  above,  it  does  not  follow  that  because  we  can 
derive  closed  form  expressions  in  principle,  we  will  be  able  to  do  such  in 
practice.  In  fact,  it  is  much  too  time  consuming  to  derive  such  expressions. 
This  can  be  demonstrated  easily  enough.  Briefly,  if  250  stimuli  (tone 
sequences)  are  displayed  on  any  one  trial,  there  are  250!  different  orders  in 
which  any  one  combination  of  tone  sequences  could  be  presented.  Since  the 
analytic  expressions  for  the  response  time  depend  on  both  the  order  in  which 
the  tone  sequences  are  presented  and  their  identity,  it  simply  is  not  possible 
to  derive  all  250!  different  predictions  in  any  reasonable  time  frame. 

Fortunately,  there  is  a  way  around  the  above  problem.  Specifically,  we 
need  not  predict  response  times  for  all  250  stimuli.  We  could  predict 
response  times  for  just  the  first  five  stimuli  in  a  group.  Since  there  are  31 
characters,  there  are  31  X  30  X  29  X  28  X  27  =  20,389,320  different  predic¬ 
tions  of  response  time  we  would  need  to  make  (this  assumes  that  no  character 
appears  twice  in  the  same  group  of  five) .  Again,  this  is  obviously  way  to 
large.  However,  suppose  that  we  were  to  use  only  five  characters  throughout 
an  entire  experiment.  Each  trial  would  consist  of  a  random  permutation  of 
these  five  characters.  A  large  gap  would  appear  between  a  group,  say  3 
seconds,  enough  to  clear  the  system.  In  this  case  we  would  need  to  make  only 
5!  =  120  predictions  of  response  time.  This  may  be  small  enough  to  be  manage¬ 
able. 

Note  that  we  do  not  need  a  different  closed  form  expression  for  each  of 
the  120  different  predictions  of  response  time.  The  closed  form  expression 
remains  the  same  across  predictions.  What  changes  are  the  values  of  the 
parameters  we  substitute  into  the  closed  form  expression. 
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Computations.  We  need  now  to  represent  processing  as  an  OP  diagram  in 
order  to  show  how  we  would  compute  the  desired  closed  form  expression  for 
response  time.  Consider  just  the  NOCOPB  model  (it  is  simplest) .  And  consider 
just  one  path  in  the  evolution  of  the  system.  Suppose  that  five  stimuli  are 
presented,  tj^,...,t^.  Label  the  arrival,  perception,  critical  time,  recogni¬ 
tion  and  execution  processes  associated  with  the  i^^  tone  sequence 
respectively  as  a^,  p^,  c^,  r^,  and  e^.  Then,  an  example  path  through  the  OP 

diagram  would  appear  as  below: 
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^44 

13 

e  • . 

44 

Note  that  we 

have  not  represented  the  pre-active  set  simply  : 

would  be  too 

cumbersome. 

Briefly, 

.  the  figure  should  be  interpreted  as  follows. 

arrival  server  is  busy  processing  the  first  tone  sequence  t^  on  the  arrival 
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server.  In  state  S2  this  arrival  completes  (so  a^  is  placed  in  the  completed 
set)  and  the  first  tone  sequence  is  passed  on  to  the  perception  server  (p^^)  . 

A  process  (c^)  is  also  placed  on  the  perception  server  and  represents  the 
critical  time.  And  finally,  the  next  arrival  (a2)  is  now  placed  on  the  ar¬ 
rival  server.  Processing  continues  in  this  fashion  for  all  five  arrivals. 

Note  that  in  state  s^  we  assume  that  the  perception  process  (p^)  com¬ 
pletes  before  the  process  (c^^)  which  represents  the  critical  time.  Thus  c^  is 

placed  in  the  other  set  since  its  processing  is  interrupted.  Just  the  op¬ 
posite  happens  in  state  Sg.  Here  the  critical  perception  time  associated  with 

tone  sequence  t^  is  shorter  than  the  duration  of  the  perception  time  for  the 

same  tone  sequence.  Thus,  the  fifth  tone  sequence  is  sent  as  a  period  element 
string  (44)  to  the  recognition  server. 

It  is  now  a  straightforward  matter  to  derive  the  seven  dependent  vari¬ 
ables.  This  requires  )cnowing  how  to  compute  the  probability  that  a  particular 
path  is  ta)cen  and  the  expected  duration  of  the  path.  The  computations  are 
reasonably  simple  if  we  assume  that  the  durations  of  the  various  processes  are 
independent,  exponentially  distributed  random  variables.  In  this  case,  we 
need  only  one  parameter  to  describe  each  process  x^,  the  rate  parameter,  say 

1. . 

1 


The  probability  of  ta)cihg  a  path  from  the  start  state  up  to  and  including 

state  s  is  simply  the  product  of  i  conditional  probabilities,  the  first  of 
“i 

which,  P(s  |s  ) ,  is  the  ratio  of  the  rate  parameter  of  the  process  which 
“2  “1 

completes  when  a  transition  is  made  from  state  s  to  state  s  to  the  sum  of 

“1  “2 

the  rate  parameters  of  the  processes  current  in  state  s^  ,  the  second  of 


which,  P(s  Is  ) ,  is  the  ratio  of  the  rate  parameter  of  the  process  which 

“3  “2 
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completes  when  a  transition  is  made  from  state  s  to  state  s  to  the  sum  of 

“2  “3 

the  rate  parameters  of  the  processes  current  in  state  s^  ,  and  so  on.  The 

time  on  average  it  takes  to  traverse  a  path  up  to  and  including  state  s  is  a 

“i 

sum  of  i  conditional  expectations,  the  first  of  which,  E[T  |a],  is  the 

“l 

reciprocal  of  the  sum  of  the  rate  parameters  in  state  ,  the  second  of 
which,  E[T  |a],  is  the  reciprocal  of  the  sum  of  the  rate  parameters  of  the 

“2 

processes  current  in  state  s^  ,  and  so  on. 

Given  the  above  information,  we  can  go  on  to  compute  the  moments  of  the 
response  time  (Fisher  and  Goldstein,  1983;  Goldstein  and  Fisher,  1991) .  We 
will  not  develop  these  computations  here  since  to  do  so  would  require  a  fair 
amount  of  work,  work  which  may  well  change  as  we  become  more  familiar  with  the 
sort  of  model  we  eventually  want  to  fit. 

MODEL  EVALUATION 

We  were  most  interested  in  whether  the  LimCOPB  or  NoCOPB  models  better 
fit  the  results.  In  order  to  determine  this,  we  performed  some  preliminary 
analyses.  Below,  we  describe  these  preliminary  analyses  for  the  two  models. 

Briefly,  we  fit  the  results  from  19  subjects.  The  subjects  were  trained 
at  the  US  Army  Intelligence  School,  Fort  Devons,  Massachusetts.  The  subjects 
had  completed  the  character  learning  phase  and  were  beginning  the  speed  build¬ 
ing  phase  at  the  starting  rate  of  six  groups  per  minute  (5  characters  appeared 
in  each  group) .  For  these  subjects,  we  had  for  each  of  the  tone  sequences  t^, 

i  =  1  -  31,  the  probability  that  the  subjects  made  a  character  response  r^, 

j  =  1,...,31,  a  period  response,  or  no  response.  And  we  had  for  each  of  the 
tone  sequences  t^,  i  =  1  -  31,  the  time  on  average  that  it  took  the  subjects 

to  make  a  character  response  r^,  j  =  1,  ...,31,  or  a  period  response.  The 

relevant  data  are  given  in  Appendix  B. 
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LimCOPB 


Parameter  Search.  There  are  five  free  parameters  with  the  LimCOPB  model. 
We  did  not  engage  in  an  extensive  search  of  the  parameter  space  since  the  fit 
proved  quite  good  with  our  initial  estimates.  Specifically,  we  found  a  good 
fit  with  the  following  parameter  settings  for  a^,  a^,  and  8: 


Process 

Perception  Server 

Perception  Decay 
Recognition  Server 

Recognition  Decay 
Execution  Server 

Execution  Decay 


Shape  Scale  Mean  Standard  Deviation 
{parameters  determined  from  data;  depended  on 
individual  tone  sequence;  see  Appendix  A} 
{not  applicable  since  there  is  no  buffer} 

15  (a^)  10  (P^)  150  39 

{not  applicable  since  there  is  no  buffer} 

15  (Og}  lO(Pg)  150  39 

{not  applicable  since  there  is  no  buffer} 


The  critical  time  5  we  set  at  1500  ms  after  trying  several  alternatives.  The 
other  63  parameters  (32  parameters  to  describe  the  arrival  times  for  each 
character;  31  parameters  to  describe  the  errors)  were  fixed  by  the  confusion 
matrices  in  Appendix  2. 


With  the  above  parameter  settings,  we  ran  the  simulation.  The  predicted 
and  observed  probabilities  are  displayed  below.  The  simulated  results  are 
based  on  the  output  from  500  trials: 


P (CORRECT),  P(INC  CHAR),  P (PERIOD) ,  P (NO  RESP) 

Predicted  0.848000  0.064000  0.068000  0.020000 

Observed  0.849500  0.049745  0.069752  0.031003 


Superficially,  the  agreement  loo)cs  reasonable  enough.  We  also  obtained 
predictions  of  the  relevant  response  times  which  are  listed  below  together 
with  the  observations: 


RT (CORRECT),  RT(INC  CHAR),  RT (PERIOD),  RT (NO  RESP) 
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Predicted 

Observed 


1184.7144  1162.0707 

1230.8064  1528.5364 


1922,2770 

2080.2903 


0.0000 

0.0000 


Again,  the  match  looks  good  except  for  the  response  times  to  the  incorrect 
characters . 

Incorrect  Character  Response  Times.  The  fact  that  we  did  not  fit  the 
response  times  for  the  incorrect  characters  arises  directly  from  the  way  in 
which  we  assign  these  response  times.  In  particular,  note  that  the  time  to 
process  tone  sequence  t^  on  the  perception  server  is  based  on  the  time 

T(ri|ti)  that  it  takes  to  make  a  correct  character  response.  If  the  percep¬ 
tion  process  for  tone  sequence  t^  takes  longer  than  the  critical  perception 

time,  1500  ms  in  the  above  run,  then  it  is  assumed  to  be  transformed  into  the 
element  string  corresponding  to  the  period.  Since  only  times  greater  than 
1500  ms  go  into  the  average  of  the  period  response  time,  this  response  time 
will  be  at  least  as  great  as  1500  ms  even  though  the  underlying  distribution 
from  which  the  times  were  drawn  may  have  a  mean  much  smaller  than  1500  ms. 
However,  the  predicted  incorrect  character  response  times  will  equal  the  cor¬ 
rect  character  response  times  since  they  come  from  one  and  the  same  truncated 
distribution. 

Specifically,  if  the  perception  process  for  tone  sequence  t^  takes  less 
than  1500  ms,  then  the  tone  sequence  is  assumed  to  be  transformed  correctly 

A 

with  probability  p{r^|tj^)  and  incorrectly  with  one  minus  this  probability.  No 

adjustment  is  made  to  the  incorrect  response  times  however.  Thus,  we  find  a 
failure  to  predict  the  incorrect  response  times. 

We  need  in  future  work  to  extend  the  LimCOPB  model  so  that  it  can  predict 
the  incorrect  response  times.  Two  possibilities  suggest  themselves.  First, 
we  could  introduce  a  second  critical  time  with  a  value  less  than  the  first 
critical  time.  Thus,  say  we  have  5^^  and  62'  where  6^  <  {2-  Then,  a  tone  se¬ 
quence  would  be  transformed  into  a  period  element  string  if  the  perception 
service  time  were  greater  than  82.  A  tone  sequence  would  be  transformed  into 
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an  incorrect  element  string  if  the  perception  service  time  were  greater  than 
and  less  than  or  equal  to  62  •  And  a  tone  sequence  would  be  transformed 

into  the  correct,  associated  element  string  if  the  perception  service  time 
were  less  than  or  equal  to  6^^. 

Alternatively,  we  could  assign  one  distribution  to  the  correct  response 
times  and  a  second  distribution  to  the  incorrect  response  times.  Then,  as 
above  we  would  assume  that  any  response  time  greater  than  the  critical 
response  time  eventuated  in  a  period  response.  Currently,  we  are  exploring 
both  ways  of  handling  the  incorrect  character  response  times. 

Serial  Position  Information.  We  should  note  before  concluding  this  sub¬ 
section  that  the  simulation  predicts  not  only  the  seven  dependent  variables, 
but  also  serial  position  information.  Below,  we  see  that  there  there  appears 
to  be  a  primacy  effect,  but  not  a  recency  effect,  within  groups  of  five: 

SERIAL  POSITION  INFORMATION 
47  43  40  42  40 

This  is  what  we  would  expect  since  the  first  item  in  a  group  is  more  likely  to 
meet  an  unoccupied  perception  server  than  is  the  last  item  in  a  group. 

We  also  generated  predictions  of  response  time  across  the  five  serial 
positions  within  a  group: 

SERIAL  POSITION  INFORMATION 

1142.072  1130.202  1273.773  1170.171  1219.631 

Since  there  is  no  room  in  any  of  the  buffers,  any  item  that  made  it  through 
the  system  should  have  had  a  clear  path  all  the  way.  Thus,  we  would  not  ex¬ 
pect  to  see  serial  position  effects.  Without  further  statistical  analysis  it 
is  difficult  to  determine  whether  the  above  differences  are  significant. 
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NoCOPB 


Parameter  Search.  As  might  be  expected,  the  NoCOPB  model  fit  the  results 
every  bit  as  well  as  the  LimCOPB  model.  This  is  expected  since  at  the  slower 
rates  there  is  not  as  much  need  for  copying  behind.  Below,  we  list  the 
results  of  three  runs,  one  with  a  critical  time  of  1600  ms,  one  with  a  criti¬ 
cal  time  of  1400  ms,  and  one  with  a  critical  time  of  1300  ms.  The  first  run 
is  with  5  =  1600: 

CRITICAL  TIME  AT  PERCEPTION  SERVER  =  1600.000 


P (CORRECT) , 

P(INC  CHAR),  P (PERIOD), 

P(NO  RESP) 

Predicted 

0.892000 

0.056000 

0.004000 

0.048000 

Observed 

0.849500 

0.049745 

0.069752 

0.031003 

RT (CORREC) , 

RT(INC  CHR),  RT (PERIOD) 

,  RT(NO  RSP) 

Predicted 

1204.0682 

1220.6215 

2086.5015 

0.0000 

Observed 

1230.8064 

1528.5364 

2080.2903 

0.0000 

SERIAL  POSITION  INFORMATION:  NUMBER  CORRECT 

47  43  45  43  45 

SERIAL  POSITION  INFORMATION;  CORRECT  TIME 
1167.336  1214.751  1271.863  1187.347  1180.407 

Note  that  the  predicted  response  times  (except  for  the  incorrect  characters) 
match  well  the  observed  response  times.  However,  the  observed  probability  of 
a  period  response  (0.06975)  is  seriously  underestimated  (0.004). 

If  we  want  to  increase  this  probability  we  need  to  decrease  the  critical 
perception  time,  which  we  do  in  the  run  below; 

CRITICAL  TIME  AT  PERCEPTION  SERVER  =  1400.000 

P (CORRECT),  P(INC  CHAR),  P (PERIOD),  P (NO  RESP) 

Predicted  0.832000  0.040000  0.068000  0.060000 

Observed  0.849500  0.049745  0.069752  0.031003 
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RT(CORREC),  RT(INC  CHR) ,  RT (PERIOD),  RT (NO  RSP) 

Predicted  1152.5948  1265.1078  1830.5255  0.0000 

Observed  1230.8064  1528.5364  2080.2903  0.0000 

SERIAL  POSITION  INFORMATION:  NUMBER  CORRECT 

43  43  42  42  38 

SERIAL  POSITION  INFORMATION:  CORRECT  TIME 
1154.037  1142.396  1125.010  1159.207  1185.682 

The  predicted  probability  (0.068)  of  a  period  response  now  matches  closely  the 
observed  probability.  But  note  that  in  raising  the  predicted  probability  of  a 
period  response  we  have  lowered  the  average  period  response  time  to  the  point 
where  it  no  longer  is  closely  tied  to  the  observed  period  response  time. 

Thus,  the  model  is  clearly  sensitive  to  variations  in  the  parameters. 

Finally,  we  would  expect  even  worse  fits  were  we  to  decrease  the  critical 
response  time  still  further.  In  the  run  below,  we  lower  this  time  to  1300  ms: 

CRITICAL  TIME  AT  PERCEPTION  SERVER  »  1300.000 

P (CORRECT),  P(INC  CHAR),  P (PERIOD),  P (NO  RESP) 

Predicted  0.724000  0.044000  0.124000  0  108000 

Observed  0.849500  0.049745  0.069752  0.031003 

RT(CORREC),  RT(INC  CHR),  RT(PERIO),  RT(NO  RSP) 

Predicted  1169.4633  1225.6049  1854.8347  0.0000 

Observed  1230.8064  1528.5364  2080.2903  0.0000 


SERIAL  POSITION  INFORMATION:  NUMBER  CORRECT 

38  31  37  37  38 

SERIAL  POSITION  INFORMATION:  CORRECT  TIME 
1110.812  1167.096  1199.870  1182.710  1187.541 

Note  that  the  probability  of  a  correct  response,  the  probability' of  a  period 
response,  and  the  probability  of  no  response  are  now  seriously  off.  It  is 
clear  why  the  probability  of  a  period  response  should  increase  and  why  the 
probability  of  a  correct  response  should  decrease.  It  is  not  clear  to  us  why 
the  probability  of  no  response  should  also  increase.  Specifically,  the  dis¬ 
tribution  of  the  perception  service  times  does  not  change  as  a  function  of  the 
critical  perception  time,  nor  do  either  the  recognition  or  execution  service 
times  change.  Thus,  the  probability  of  a  period  response  ought  to  remain 
roughly  the  same.  We  need  to  determine  whether  the  above  variation  is  random 
or  systematic,  and  if  systematic  why  it  is  occurring. 

CONCLUSION 

The  simulation  developed  here  allows  the  user  to  predict  for  the  Morse 
code  copy  task  the  probability  of  a  correct  response,  the  probability  of  an 
incorrect  response,  the  probability  of  a  period  response,  the  probability  of 
no  response,  and  the  probability  of  a  correct  response  for  each  of  the  five 
serial  positions  in  a  group.  The  simulation  also  allows  the  user  to  predict 
the  time  that  it  takes  to  execute  a  correct  response,  the  time  that  it  takes 
to  execute  an  incorrect  response,  the  time  that  it  takes  to  execute  a  period 
response,  and  the  time  that  it  takes  to  execute  a  correct  response  for  charac¬ 
ters  occurring  in  each  of  the  five  serial  positions  of  a  group. 

In  the  most  general  case  the  simulation  required  the  estimation  of  77 
parameters,  63  of  which  were  fixed  by  the  data.  These  parameters  were  needed 
to  de.'^cribe  the  distribution  of  the  duration  of  the  service  and  decay  times  at 
each  of  the  perception,  recognition  and  execution  stages,  the  number  of 
stimuli  which  could  fit  into  the  three  buffers,  and  the  critical  perception 
time. 
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No  models  to  date  have  allowed  for  as  wide  a  specification  of  input 
parameters  or  for  as  rich  a  set  of  predictions.  These  models  should  make  it 
possible  to  explore  in  detail  a  number  of  hypotheses  about  the  origins  of  the 
training  problems  with  learning  to  copy  Morse  code  at  higher  input  rates. 

One  such  hypothesis  was  tested,  in  particular,  the  hypothesis  that  sub¬ 
jects  early  in  practice  could  not  process  two  stimuli  at  the  same  time  (the 
copy  behind  hypothesis) .  Several  variants  of  the  most  general  model  were 
created  to  test  this  hypothesis.  One,  the  LimCOPB,  allowed  copy  behind;  the 
other,  the  NoCOPB,  did  not.  The  models  fit  about  equally  well.  This  is  not 
all  that  surprising  given  that  copy  behind  is  not  particularly  useful  early  on 
since  the  intercharacter  interval  is  relatively  large.  It  still  remains  to  be 
determined  whether  at  the  faster  rates  subjects  who  make  it  through  the  train¬ 
ing  learn  to  copy  behind  whereas  those  that  don't  fail  to  learn  this  skill. 
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APPENDIX  A 

PARAMETERS  OF  GenCOPB,  LimCOPB  AND  NoCOPB  MODELS 

GenCOPB  Model .  There  are  a  total  of  77  parameters  in  the  GenCOPB  model. 
To  begin,  consider  the  determination  of  the  duration  of  the  perception,  recog¬ 
nition  and  execution  processes.  Thirty-six  parameters  are  needed  to  describe 
the  distributions  of  these  three  processes,  only  four  of  which  are  free.  We 
assume  that  the  time  T(r^|t^)  that  it  takes  to  generate  character  r^  as  a 

response,  given  tone  sequence  t^  was  presented,  is  the  sum  of  the  time 
Tp(eitti)  that  it  takes  to  perceive  a  tone  sequence  t^  as  element  string  e^, 
given  that  tone  sequence  was  presented,  plus  the  time  T^  that  it  takes  to 
recognize  element  string  e^  plus  the  time  T^  that  it  takes  to  organize  and  ex¬ 
ecute  a  response. 

T(rilti)  =  Tp(ejt.)  +  T^  +  T^. 

Note  that  we  assume  that  the  recognition  process  is  errorless  and  that  the 
duration  of  this  process  does  not  vary  from  one  element  string  to  the  next. 
Thus,  we  can  drop  both  the  indexed  element  string  e^  which  is  output  to  the 

recognition  process  and  the  character  code  c^  which  is  generated  from  this 

process.  Similarly,  we  assume  that  the  execution  process  is  errorless  and 
that  the  duration  of  this  process  does  not  vary  from  one  character  code  to  the 
next.  Thus,  we  can  drop  both  the  indexed  character  code  c^  and  the  response 

r^  which  is  generated  from  the  character  code. 

Since  we  assume  that  the  duration  T^  of  the  recognition  process  and  dura¬ 
tion  T^  of  the  execution  process  do  not  depend  on  the  identity  of  the  stimulus 

which  is  being  processed,  we  need  only  four  parameters  to  describe  the  dis¬ 
tribution  of  these  durations:  a  shape  parameter  and  scale  parameter  for 

the  recognition  stage  and  a  shape  parameter  and  scale  parameter  p^  for  the 

execution  stage.  These  four  parameters  are  free. 


A  -  1 


We  assume  that  the  distribution  of  a  gamma  with  individual 

shape  parameter  shape  parameters  Opli)  and  common  scale  parameter  Note 


that  we  could  have  individualized  both  the  shape  and  scale  parameters,  but 
decided  not  to  do  so  on  this  first  attempt.  The  common  scale  parameter  was 
selected  so  that  the  variance  of  the  predicted  observations  was  similar  to 
that  of  the  observed.  The  individual  shape  parameters  were  selected  so  that 


the  equation,  T(r^|t^)  -  E[T^]  -  E[T^]  =  ap(i)pp,  was  true,  where  T(r^|t^) 

is  the  average  correct  response  latency  (obtained  from  the  confusion  matrix) , 

E[T  ]  =  a  p  and  E[T  ]  =  a  p  ,  i.e., 

I  r-'  r'^r  '■  e  e'^e 


ap{i)  =  (T(r. |t^)  -  E[T^]  -  E[Tg]}/pp,  i  =  1,...,31  (Al.l) 

Thus,  although  there  are  32  parameters  [31  shape  parameters,  0^(1), 
i  =  1,...,31,  and  one  scale  parameter,  p^)  needed  to  describe  the  distribu¬ 
tions  of  the  durations  of  the  31  correct  perception  times,  these  parameter.' 
are  all  fixed  by  the  confusion  matrix  data  (for  response  times) . 


The  next  10  parameters  are  free.  To  begin,  consider  the  parameters 
needed  to  described  the  decay  time  in  the  perception,  recognition,  and  execu¬ 
tion  buffers.  We  assume  that  the  distribution  of  the  decay  time  does  not 
depend  on  the  identity  of  the  stimulus  being  processed  or  the  fact  that  the 
stimulus  is  a  letter  as  opposed  to  a  period.  Thus,  we  need  a  total  of  6 
parameters,  2  for  the  perception  decay  time  (a*  and  p*) ,  2  for  the  recognition 

decay  time  (a*  and  p^  and  two  for  the  execution  decay  time  (a*  and  p*) . 

Next,  we  need  three  parameters  for  the  capacity  of  the  perception,  recognition 
and  execution  buffers.  We  will  label  these  7^,  and  y^.  Finally,  we  will 

also  need  a  parameter  for  the  critical  time  which  we  will  label  5 .  Again, 
these  last  ten  parameters  are  free. 


To  complete  the  GenCOPB  model,  we  need  31  parameters  to  determine  when  a 
character  is  perceived  correctly  and  when  it  is  perceived  incorrectly  as 
another  character,  given  that  the  duration  of  the  perception  process  is  less 
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than  the  critical  time  8 .  These  parameters  are  fixed  by  the  data  in  the  con¬ 
fusion  matrix  (errors).  Specifically,  the  probability  p(e^|t^)  that  element 

string  e^  is  generated  at  the  perception  stage,  given  that  tone  sequence  t^ 

was  presented  and  that  the  duration  of  the  perception  process  was  less  than  8, 

A 

is  set  equal  to  the  observed  probability  p(r^|t^)  that  response  r^  is  made, 

given  that  tone  sequence  t^  was  presented.  The  probability  that  element 

string  e^  is  generated  at  the  perception  stage,  given  that  tone  sequence  t^ 

was  presented  (i  /  j)  and  that  the  duration  of  the  perception  process  was  less 
than  the  critical  time  8,  is  simply  the  complement  of  the  above,  i.e.,  simply 

A 

1  -  p(r^lt^).  In  the  simulation,  the  determination  on  each  trial  of  whether  a 

correct  or  incorrect  element  string  was  generated  was  made  by  sampling  from  a 
uniform  [0,1]  distribution.  If  the  value  sampled  lay  in  the  interval 

A 

[0,p (r^ |t^) ] ,  then  the  tone  sequence  t^was  transformed  into  the  correct  ele¬ 
ment  string;  otherwise,  it  was  assumed  that  the  element  string  was  incorrect. 
No  attempt  was  made  to  determine  which  of  the  incorrect  element  strings  was 
produced. 

Note  that  we  do  not  require  additional  parameters  to  describe  the  output 
of  the  recognition  and  execution  stages.  This  is  because  we  assume  that  the 
probability  p(c^|e^)  that  character  code  c^^  is  generated  by  the  recognition 

stage,  given  that  element  string  e^  is  input  to  the  recognition  stage,  is 

equal  to  1,  as  is  the  probability  p(r^lc^)  that  the  execution  stage  generates 

response  r^,  given  that  character  code  c^  was  generated. 

LimCOPB  and  NoCOPB  Models.  In  order  to  reduce  the  number  of  free 
parameters,  we  made  several  simplifying  assumptions.  Specifically,  we  assumed 
that  there  was  no  room  in  the  perception,  recognition  or  execution  buffers. 
And,  therefore,  the  three  free  buffer  size  parameters  and  six  free  decay 
parameters  were  no  longer  necessary.  Thus,  the  number  of  free  parameters  was 
reduced  by  nine,  from  14  to  5.  The  number  of  fixed  parameters  remained  at  63. 
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APPENDIX  B: 

ERROR  AND  LATENCY  CONFUSION  MATRICES 


Below,  we  report  the  results  from  the  19  subjects  we  used  to  determine 
the  relative  fit  of  the  LimCOPB  and  NoCOPB  models.  First,  we  list  the  data  in 
the  error  confusion  matrix.  The  first  entry  in  row  1  (0.9095)  is  the  prob¬ 
ability  that  tone  sequence  t^  was  recalled  as  character  r^;  the  second  entry 

(0)  in  row  1  is  the  probability  that  tone  sequence  t^  was  recalled  as  charac¬ 
ter  Z2‘,  and  so  on  through  the  41^^  entry.  The  42^^  entry  in  the  i^^  row  is 

the  probability  that  the  i^^  character  is  recalled  correctly.  For  the  first 
row,  it  is  the  probability  (0.9095)  that  the  first  entry  is  recalled  cor- 

rectly.  The  43^^  entry  in  the  i  row  is  the  probability  that  a  tone  sequence 
is  recalled  as  a  character  other  than  the  one  which  is  associated  with  it.  In 
the  first  row,  this  probability,  0.0181,  is  equal  to  the  product  (7) (.0016). 

The  44^^  entry  in  the  i^^  row  is  the  probability  of  a  period  response  to  the 

i  character.  And  the  45  entry  in  the  row  is  the  probability  of  that  no 

response  is  made  to  the  i^^  character.  Note  that  the  first  "row"  includes 
the  four  lines  indented  underneath  the  first  line;  space  did  not  permit  the 
printing  of  all  entries  in  a  row  of  the  confusion  matrix  as  a  single  line: 
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Next,  we  present  the  results  from  the  latency  confusion  matrix, 
that  the  only  difference  between  these  results  and  the  results  above 
there  is  not  an  entry  for  the  latency  of  no  response. 
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0 

0 

0 

1177 

1622.2775 

1771 

0 

896 

0 

0  0  1546 

0  3808 

0 

0 

0 

0 

0 

983  0 

2906 

0  0 

0 

0 

0 

1399 

1203 

1312  0 

0 

0  0 

0 

0 

0 

1404 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1312 

1616.3990 

2166 

0 

0  1126 

1447 

867  870 

0 

0  2203 

0 

0 

1416 

2125 

0  0 

0 

1171  0 

0 

0 

0 

0 

2279 

0  1358 

3190 

1606  943 

1339 

0 

0 

3041 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1358 

1602.5048 

2055 

0 

0  957 

1355 

0  0 

0 

0  0 

0 

1706 

865 

2372 

0  0 

2023 

1146  1734 

0 

0 

1522 

1730 

3811 

0  1992 

1157 

0  0 

0 

0 

0 

0 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1157 

1628.7236 

2140 

0 

1558  1771 

0  727 

0 

0  1776 

0 

0 

0 

0 

0 

0  0 

1901 

1308  1699 

0 

0 

0 

1986 

0 

0  1205 

0 

1338  0 

2597 

0 

1892 

1975 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1338 

1703.9824 

2095 

0 

1274 

0  1132  0 

0 

0  0 

0 

0 

0 

1619 

1516 

0  0 

2035 

0  0 

0 

0 

0 

0 

0 

0  0 

0 

1256  1303 

0 

0 

1005 

0 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1303 

1446.3634 

1962 

0 

0  0 

0 

648  1293 

1395 

438  0 

0 

0 

0 

1286 

0  0 

0 

1901  0 

0 

765 

0 

0 

989 

2412  2186 

0 

0  1948 

1321 

0 

935 

1419 

0 

0  0 

0 

0  0 

0 

0 

0 

0 

1321 

1361.5017 

1740 

0 
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0  1131 


0  1506  0  0  0  0 

0000  1436 
0  0  6043  0  0 

0  0  0  0  0 

0  1111  1548.1082  2201 
1161  3424  0  1375  0 

0000  1283 
0  0  0  0  0 

5232  0  0  0  0 

0  1373  179.7826  2146 
00000  1788 
0  1278  000 

1081  0  1384  902  1136 

1216  0  0  0  0 

0  1216  1398.9321  1943 


0  1270 
000  2570  0 

000  1111  1514 
0  0  0  0  0 

0 

0  2036  1930  0  1591 

1940  0  978  0  0 

000  1449  1373 
0  0  0  0  0 

0 

000  2000 
0000  1855 
1043  0  0  0  0 

0  0  0  0  0 

0 
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APPENDIX  C: 

THE  FORTRAN  SIMULATION  CODE 


MOST  RECENT  UPDATE:  11/21/91  (APPROXIMATE) 

THIS  PROGRAM,  MORSE. FOR,  WRITTEN  8/19/91  ON  VAX  AT  UNIVERSITY  OF 
MASSACHUSETTS  COLLEGE  OF  ENGINEERING  FOR  DR,  ROBERT  WISHER 
BY  DONALD  L.  FISHER.  THE  PROGRAM  SIMULATES  THE  BEHAVIOR  OF 
AN  OPERATOR  TRYING  TO  COPY  MORSE  CODE, 

THE  UNDERLYING  MODEL  IS  ONE  IN  WHICH  THE  FOLLOWING  ASSUMPTIONS  ARE  MADE: 

1)  THERE  ARE  FOUR  PROCESSES:  ARRIVAL,  AUDITORY  PERCEPTION,  CHARACTER 
RECOGNITION  AND  MOTOR  ORGANIZATION  AND  EXECUTION.  THE  DURATION 

OF  THE  ARRIVAL  PROCESS  IS  A  CONSTANT  EQUAL  TO  THE  DURATION 
OF  THE  PARTICULAR  MORSE  CHARACTER  THAT  IS  SENT.  THE  DURATION  OF  THE 
PERCEPTION,  RECOGNITION  AND  EXECUTION  PROCESSES  ARE  EACH  GENERAL 
GAMMA  RANDOM  VARIABLES. 

2)  THERE  ARE  FOUR  BUFFERS:  THE  ARRIVAL  BUFFER  IS  STIMBUF  AND  IS 
UNLIMITED.  IT  CONTAINS  BOTH  THE  CHARACTER  AND  INTERCHARACTER 
STIMULI,  THE  PERCEPTION,  RECOGNITION  AND  EXECUTION  BUFFERS  ARE 
EACH  SPECIFIED  BY  THE  USER  IN,  RESPECTIVELY,  INFO(2,2),  INFO(2,3) 

AND  INFO (2, 4) .  NOTHING  DECAYS  OUT  OF  THE  ARRIVAL  BUFFER.  THE 
DURATION  OF  THE  DECAY  PROCESSES  IN  THE  THREE  REMAINING  BUFFERS 

IS  MODELED  AS  A  GENERAL  GAMMA  RANDOM  VARIABLE. 

3)  THE  EXECUTION  BUFFER  HOLDS  PERIODS,  NOT  CHARACTERS.  SPECIFICALLY, 
IF  A  STIMULUS  EXITS  THE  RECOGNITION  STAGE  AND  THE  EXECTUTION  SERVER 
IS  FULL,  THEN  THE  STIMULUS  ENTERS  THE  EXECUTION  BUFFER  AS  A  PERIOD. 

THE  SIZE  OF  THE  EXECTUTION  BUFFER  SHOULD  BE  KEPT  RELATIVELY  LARGE. 
SIMILARLY,  BECAUSE  WE  DO  NOT  WANT  RAPID  DECAY  OF  PERIODS,  THE  DURATION 
OF  A  PERIOD  IN  THE  EXECUTION  BUFFER  SHOULD  BE  SET  RELATIVELY  LARGE. 


DEFINE  VARIABLES: 

NCHR:  NUMBER  OF  DIFFERENT  CHARACTERS  IN  INPUT  STREAM:  26  (JUST  LETTERS) , 
31  (LETTERS  +  SPECIAL  CHARACTERS),  41  (LETTERS,  SPECIAL 
CHARACTERS  AND  DIGITS) ;  INITIALIZED  IN  SUBROUTINE  SINV 

NCHRBUF:  DIMENSION  OF  CHRBUF;  MUST  BE  GREATER  THAN  NGROUP*NREP 

NCNT:  SET  EQUAL  TO  INFO(2,l)  AT  START;  INFO(2,l)  IS  INITIALIZED 
IN  THE  FILE  MORSE. IN 

NELEM:  NEXT  ELEMENT  IN  STIMULUS  BUFFER  (STIMBUF) ;  READ  IN  FROM 
SUBROUTINE  SINV;  ALWAYS  START  AT  ELEMENT  0, 
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NGPM:  GROUPS  PER  MINUTE;  READ  IN  FROM  SUBROUTINE  SINV 

NGROUP:  NUMBER  OF  CHARACTERS  IN  A  GROUP 

NOUT:  SET  EQUAL  TO  1  IF  FULL  OUTPUT  (INFO, STATEC, STATED, STATEX, 
STATEP,STATEY)  DESIRED;  OTHERWISE  SET  TO  0;  READ  IN  FROM 
SUBROUTINE  SINV 

NPER:  THE  NUMBER  OF  PERIOD  RECODINGS  PRODUCED  WHEN  A  STIMULUS 
LEAVES  THE  PERCEPTION  SERVER. 

NREP:  NUMBER  OF  REPETITIONS  OF  A  GROUP 

NSTIM:  NUMBER  OF  "STIMULI"  =  2*NGROUP*NREP  (NOTE  THAT  EVERY  CHARACTER 
IS  FOLLOWED  BY  EITHER  AN  INTERCHARACTER  BLANK  OR  AN  INTER¬ 
GROUP  BLANK) ;  NOTE  THAT  NELEM  CAN  REACH  NSTIM  AS  AN  UPPER 
LIMIT 

RCRIT:  IF  IT  TAKES  LONGER  THAN  RCRIT  MILLISECONDS  TO  SERVICE  A 
STIMULUS  IN  THE  PERCEPTION  NODE,  THEN  THE  STIMULUS  IS 
RECODED  AS  A  PERIOD. 


DEFINE  MATRICES: 

ALPHA (2, 4),  BETA (2, 4);  IT  IS  ASSUMED  THAT  EACH  SERVICE  TIME  AND  EACH 
DECAY  TIME  (EXCEPT  ARRIVAL)  IS  A  GAMMA  RANDOM  VARIABLE 
WITH  PARAMETERS  ALPHA  (a)  AND  BETA  (b) ,  I.E., 


-a  a-1  -x/b 

f(x)  =  [b  X  e  ] /gaitima  (a) . 

IN  THIS  CONTEXT,  THE  PARAMETERS  FOR  THE  FOUR  SERVICE  AND 
DECAY  TIMES  ARE: 


ARRIVAL  TIME  :  ALPHA (1,1), 
HOWEVER,  COULD  CHANGE 
ARRIVAL  DECAY:  ALPHA (2,1), 
AGAIN,  COULD  CHANGE 


PERCEP 

PERCEP 

RECOG 

RECOG 

EXEC 

EXEC 


TIME 

DECAY 

TIME 

DECAY 

TIME 

DECAY 


ALPHA(1,2) , 
ALPHA(2,2) , 
ALPHA (1,3), 
ALPHA(2,3), 
ALPHA(1,4) , 
ALPHA(2,4) , 


BETA (1,1) 
TO  GAMMA 
BETA (2,1) 

BETA (1,2) 
BETA (2, 2) 
BETA (1,3) 
BETA (2, 3) 
BETA (1,4) 
BETA (2, 4) 


ASSUMED  CONSTANT; 
ASSUMED  UNLIMITED; 


CHRBUF(410) :  PROGRAM  WRITTEN  TO  HANDLE  50  GROUPS  OF  5  CHARACTERS 
EACH.  NOTE  THAT  STIMBUF  HAS  BOTH  INTERACHARACTER  AND 
INTERGROUP  BLANKS  AND  THUS  SHOULD  BE  TWICE  AS  LARGE  AS 
CHARBUF.  ARRAY  OF  SIZE  410  BECAUSE  EVERY  41  CHARACTERS 
IS  GUARANTEED  TO  CONTAIN  DIGITS  1-41.  IF  31  CHARACTERS 
IS  USED,  THEN  CHRBUF  IS  FILLED  TO  POSITION  31*10  =  310. 
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IF  26  CHARACTERS  ARE  USED  THEN  ARRAY  IS  FILLED  TO  26*10  =  260. 
SINCE  260  IS  GREATER  THAN  STANDARD  50*5  =  250  CHARACTERS, 

I  WENT  WITH  THE  RATHER  ODD  410. 

CONMAT(41,45)  : 

A)  INITIALLY,  IN  ROW  I,  COLUMN  J  IS  PROBABILITY  THAT  CHARACTER 
I  WAS  SENT  AND  CHARACTER  J  WAS  TYPED.  NOTE  THAT  PROBABILITIES 
IN  A  ROW  SHOULD  ADD  TO  ONE.  HOWEVER,  THIS  DOES  NOT  MEAN  THAT 
AN  ELEMENT  IN  ROW  I  IS  NEVER  MISSED  OR  TYPED  AS  A  PERIOD. 

THIS  PROBLEM  OCCURS  LATER  IN  PROCESSES  (CONMAT  IS  USED  IN 

THE  PERCPETUAL  BUFFER) .  NOTE  THAT  IN  ROW  I  AND  COLUMN  44  I 
CONDITIONAL  PROBABILITY  OF  TYPING  A  PERIOD  GIVEN  CHARACTER  j. 

AND  IN  ROW  I  AND  COLUMN  45  IS  CONDITIONAL  PROBABILITY  OF 
TYPING  NOTHING  GIVEN  CHARACTER  I.  THE  OUTPUT  INFORMATION 
ON  PERIODS  AND  NO  RESPONSES  IS  CONTAINED  IN  KILL:  THE  NUMBER 
OF  ENTRIES  IN  THE  PERCEPTION  AND  RECOGNITION  NODES  IS  THE 
NUMBER  OF  NO  RESPONSES;  THE  NUMBER  OF  ENTRIES  IN  THE 
EXECUTION  NODE  IS  THE  NUMBER  OF  PERIODS. 

B)  AFTER  READING  IN  CONMAT,  IF  P(J|I)  IS  THE  PROBABILITY  THAT 
J  IS  TYPED  WHEN  I  IS  PRESENTED,  THEN  IN  ROW  I,  COLUMN  1  WE 
PLACE  P(1|I),  IN  ROW  I,  COLUMN  2  WE  PLACE  P(1|I)  +  P(2|I), 

AND  SO  ON. 

CONFUS(41,41) :  IN  ROW  I  AND  COLUMN  J  IS  THE  NUMBER  OF  TIMES  IN  THE 
SIMULATION  THAT  CHARACTER  J  REPLACES 

,CTER  I  AT  THE 

END  OF  SERVICE  ON  THE  PERCEPTION  NODE.  THIS  INFORMATION  IS 
OBTAINED  FROM  CONMAT. 

DATERR(41,45) :  ERROR  DATA  COMES  FROM  SUBJECTS.  IN  ROW  I  AND  COLUMN  J 
IS  PROBABILITY  THAT  A  SUBJECT  RESPONDS  WITH  STIMULUS  J 
GIVEN  SUBJECT  HAS  BEEN  PRESENTED  STIMULUS  I.  I  =  1,...,41 
SINCE  NUMBER  OF  STIMULI  CAN  GO  UP  TO  41.  J  =  1,...,41, 

44,45.  J  =  44  IS  PERIOD  RESPONSE.  J  =  45  IS  NO  RESPONSE. 

J  =  42,43  ARE  USED  AS  STIMULUS  VALUES  (INTERCHARACTER  AND 
INTERGROUP  BLANKS),  NOT  RESPONSES.  HOWEVER,  IN  COLUMN  42 
I'VE  PLACED  THE  PERCENTAGE  CORRECT  RESPONSES,  I.E.,  IN 
DATERR(IA,42)  I'VE  PLACED  DATERR ( lA, lA) .  AND  IN  COLUMN 
43  I'VE  PLACED  THE  INCORRECT  CHARACTER  ONLY  ERRORS. 

THUS,  THE  PROBABILITIES  SUMMED  WITHIN  A  ROW  ACROSS 
COLUMNS  42  -  45  OUT  TO  EQUAL  1. 

INFO (2, 4):  CONTAINS  INFORMATION  ON  PARAMETERS  OF  ARRIVAL,  PERCEPTION, 
RECOGNITION  AND  EXECUTION  PROCESSES  IN,  RESPECTIVELY, 

COLUMNS  1,  2,  3  AND  4: 

INFO (1,1)  =  NUMBER  OF  CHARACTERS  ON  ARRIVAL  SERVER  (1  OR  0) 

INFO (2,1)  =  TOTAL  NUMBER  OF  ARRIVALS  IN  SESSION  LEFT;  NOTE 
THAT  ARRIVALS  INCLUDE  INTERCHARACTER  AND  INTER¬ 
GROUP  INTERVALS 

INFO (1,2)  =  NUMBER  OF  CHARACTERS  ON  PERCEPTION  SERVER  AND 
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AND  IN  PERCEPTION  BUFFER  (E.G.,  IF  THERE  IS  ONE 
CHARACTER  CURRENTLY  BEING  PERCEIVED  AND  TWO 
CHARACTERS  IN  THE  BUFFER,  THEN  INFO (1,2)  =  3) 

INFO (2, 2)  =  MAXIMUM  NUMBER  OF  CHARACTERS  WHICH  CAN  BE 
SERVICED  AND  BUFFERED  FOR  PERCEPTION 

INFO (1,3)  =  NUMBER  OF  CHARACTERS  ON  EXECUTION  SERVER  AND  IN 
RECOGNITION  BUFFER. 

INFO (2, 3)  =  MAXIMUM  NUMBER  OF  CHARACTERS  WHICH  CAN  BE 
SERVICED  AND  BUFFERED  FOR  RECOGNITION. 

INFO (1,4)  =  NUMBER  OF  CHARACTERS  ON  EXECUTION  SERVER  AND  IN 
EXECUTION  BUFFER. 

INFO (2, 4)  =  MAXIMUM  NUMBER  OF  CHARACTERS  WHICH  CAN  BE 
SERVICED  AND  BUFFERED  FOR  EXECUTION. 

INTERCdS)  ,INTERG(15)  ;  CONTAINS  IN  COLUMN  I  THE  INTERCHARACTER 
(INTERGROUP)  INTERVAL  FOR  PATES  OF,  RESPECTIVELY,  6,8,10,12 
14,16,18,20  GROUPS  PER  MINUTE. 

KEEP (41, 46):  CONTAINS  NUMBER  OF  TIMES  IN  KEEP (I, J)  CHARACTER  J  IS 
TYPED  WHEN  CHARACTER  I  WAS  PRESENTED.  NOTE  THAT  ENTRIES 
IN  KEEP(I,J)  CAN  BE  NO  GREATER  THAN  ENTRIES  IN  CONFUS(I,J) 
BECAUSE  CHARACTERS  CAN  DECAY  AT  A  NODE  OR  NOT  GAIN  ENTRY 
TO  A  NODE  AFTER  THEY  ARE  CONFUSED  IN  THE  PERCPETION  SERVER. 
NOTE  THAT  BEFORE  SUBROUTINE  SOUT,  KEEP (lA, 42) , KEEP (lA,  43) 

AND  KEEP (lA, 45)  ARE  BLANK;  KEEP (lA, 44)  CONTAINS  THE  PERIOD 
COUNT.  SUBROUTINE  SOUT  PUTS  INTO  KEEP (lA, 42)  THE  ENTRY 
IN  KEEP(IA,IA),  I.E.,  THE  NUMBER  OF  CORRECT  RESPONSE.  IT 
PUTS  INTO  KEEP (lA, 43)  THE  NUMBER  OF  INCORRECT  CHARACTER 
CHARACTER  RESPONSES.  AND  IT  PUTS  INTO  KEEP (lA, 45)  THE  NUMBER 
OF  NO  RESPONSES  (BY  SUMMING  OVER  KILL (lA, J) , J=l,  8) . 

KEEPSP(5):  KEEPSP(I)  CONTAINS  NUMBER  OF  CHARACTERS  IN  ITH  POSITION 
IN  A  GROUP  WHICH  MADE  IT  THROUGH  EXECUTION  STATE.  YIELDS 
STANDARD  SERIAL  POSITION  CURVES.  ONLY  CORRECT  CHARACTERS 
ARE  COUNTED. 

KILL (41, 8) :  CONTAINS  NUMBER  OF  TIMES  CHARACTER  IN  ROW  IS  LOST 

BECAUSE  OF  BUMPOUT,  DECAY  OR  PERIOD.  COLUMNS  1-4  CONTAIN 
ARRIVAL,  PERCEPTION,  RECOGNITION  AND  EXECUTION  BUMPOUT; 

NEXT  4  COLUMNS  CONTAIN  DECAY.  FOR  EXAMPLE,  IF  KILL (11, 3)  =  4, 
THIS  MEANS  THAT  CHARACTER  11  ON  LEAVING  THE  PERCEPTION  SERVER 
HAS  BEEN  BUMPED  4  TIMES  BECAUSE  THE  RECOGNITION  BUFFER  IS 
FULL.  IF  KILL (21, 6) =5,  THIS  MEANS  THAT  CHARACTER  21  ON  THE 
PERCEPTION  BUFFER  HAS  DECAYED  5  TIMES  BECAUSE  THE  PERCEPTION 
SERVER  IS  FULL.  NOTHING  IS  BUMPED  OUT  FROM  THE  ARRIVAL 
QUEUE  NOR  DOES  ANY  DECAY  OCCUR.  FURTHERMORE,  NOTHING  IS 
BUMPED  OUT  OF  THE  RECOGNITION  BUFFER  NOR  DOES  DECAY  WITH 
COMPLETE  LOSS  OCCUR.  RATHER,  PERIOD  IS  TYPED.  THUS, 

IF  KILL (10, 4)  =  6,  THIS  MEANS  THAT  SIX  PERIODS  WERE 
TYPED  BECAUSE  THE  CHARACTER  BUFFER  SPACE  IN  THE  EXECUTION 
BUFFER  WAS  EXCEEDED.  IF  KILL (10, 8)  =  2,  THIS  MEANS  THAT 
TWO  PERIODS  WERE  TYPED  BECAUSE  THE  DECAY  TIME  IN  THE 
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EXECUTION  BUFFER  WAS  EXCEEDED. 

DATLAT(41,45) :  CONTAINS  RESPONSE  TIMES  FROM  SUBJECTS.  LIKE  MATRIX 
ERROR  OTHERWISE  EXCEPT  FOR  NO  ENTRY  IN  COLUMN  45  (SINCE 
NO  TIMES  FOR  NO  RESPONSES) . 

STATEC (30, 4) , STATED (30,4) :  COLUMN  1  CONTAINS  INFORMATION  ON  THE 
ARRIVAL  PROCESS;  COLUMN  2  CONTAINS  INFORMATION  ON  THE 
PERCEPTION  PROCESS;  COLUMN  3  CONTAINS  INFORMATION  ON 
THE  RECOGNITION  PROCESS;  AND  COLUMN  4  CONTAINS  INFORMATION 
ON  EXECUTION  PROCESS  IN  BOTH  STATEC  AND  STATED. 

STATEC (1,1)  =  MORSE  CHARACTER  CURRENTLY  ARRIVING 
STATED (1,1)  =  REMAINING  DURATION  OF  MORSE  CHARACTER 

STATEC (1,2)  =  MORSE  CHARACTER  CURRENTLY  BEING  PERCEIVED;  NOTE 
THAT  THIS  IS  NOT  NECESSARILY  THE  CHARACTER  WHICH 
IS  SENT  IF  THERE  IS  A  CONFUSION 
STATED (1,2)  =  REMAINING  DURATION  OF  PERCEPTION  PROCESS 
STATEC (2, 2)  =  MORSE  CHARACTER  IN  PERCEPTION  BUFFER  (IF  ANY) 
STATED (2, 2)  =  REMAINING  DECAY  TIME  OF  ABOVE  CHARACTER 
STATEC (3, 2)  =  ADDITIONAL  MORSE  CHARACTER  IN  PERCEPTION 
BUFFER  (IF  ANY) 

STATED (3,2)  =  REMAINING  DECAY  TIME  OF  THIS  MORSE  CHARACTER 
(AND  SO  ON  FOR  EACH  ADDITIONAL  MORSE  CHARACTER  IN  PERCEPTION 
BUFFER) 

STATEC (1,3)  =  MORSE  CHARACTER  CURRENTLY  BEING  RECOGNIZED 
STATED (1,3)  =  REMAINING  DURATION  OF  RECOGNITION  PROCESS 
STATEC (2, 3)  =  MORSE  CHARACTER  IN  RECOGNITION  BUFFER  (IF  ANY) 
STATED (2, 3)  =  REMAINING  DECAY  TIME  OF  ABOVE  CHARACTER 
STATEC (3, 3)  =  ADDITIONAL  MORSE  CHARACTER  IN  RECOGNITION 
BUFFER  (IF  ANY) 

STATED (3, 3)  =  REMAINING  DECAY  TIME  OF  THIS  MORSE  CHARACTER 
(AND  SO  ON  FOR  EACH  ADDITIONAL  MORSE  CHARACTER  IN  RECOGNITION 
BUFFER) 

STATEC (1,4)  =  MORSE  CHARACTER  CURRENTLY  BEING  RESPONDED  TO 
STATED (1,4)  =  REMAINING  DURATION  OF  ABOVE  CHARACTER 
STATEC (2, 4)  =  MORSE  CHARACTER  IN  RESPONSE  BUFFER  (IF  ANY) 
STATED (2, 4)  =  REMAINING  DECAY  TIME  OF  ABOVE  CHARACTER 
(AND  SO  ON) 

STATEP (30,4) :  CONTAINS  IN  ROW  I  AND  COLUMN  J  POSITION  OF  CHARACTER 
IN  STREAM.  FOR  EXAMPLE,  IF  THE  GROUPS  ARE  OF  SIZE  5 
AND  THE  CHARACTER  CURRENTLY  ON  THE  RECOGNITION  SERVER  IS 
THE  SECOND  CHARACTER  IN  A  STREAM,  THEN  STATEP (1,2)  =  2. 

STATEX(30,4) ,STATEY(30,4) :  STATEX  IS  INDENTICAL  TO  STATEC  EXCEPT 
FOR  THE  FACT  THAT  THERE  IS  NO  CHANGE  TO  THE  IDENTITY  OF 
A  CHARACTER  MISPERCIEVED .  THIS  MAKES  IT  POSSIBLE  TO  TRACK 
THE  RESPONSE  TIME  TO  THE  CHARACTER  WHICH  WAS  SENT.  STATEY 
KEEPS  TRACK  OF  THE  TOTAL  TIME  THAT  A  CHARACTER  SPENDS  IN 
THE  SYSTEM.  MEASUREMENT  STARTS  FROM  THE  POINT  WHERE  THE 
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CHARACTER  HAS  ARRIVED  (I.E.,  MEASUREMENT  STARTS  WHEN  THE 
CHARACTER  ENTERS  THE  PERCEPTION  QUEUE  OR  SERVER)  AND 
CONTINUES  UNTIL  THE  RESPONSE  HAS  BEEN  EXECUTED. 

STIMBUF(820) ;  STIMBUF  CONTAINS  THE  CHARACTER,  INTERCHARACTER  CODE  OR 
INTERGROUP  CODE.  THE  SPACE  BETWEEN  CHARACTERS  IS  IDENTIFIED 
BY  THE  INTEGER  42;  THE  SPACE  BETWEEN 

’S 

IS  IDENTIFIED  BY  THE  INTEGER  43.  THE  FIRST  CHARACTER  APPEARS 
IN  STIMBUF (1),  THE  SECOND  CHARACTER  IN  STIMBUF (2),  AND  SO  ON. 

NOTE  THAT  STIMBUF  MUST  BE  TWICE  THE  SIZE  OF  CHRBUF. 

STIMDUR(43):  CONTAINS  DURATIONS  OF  1  -  26  LETTERS,  5  SPECIAL 
CHARACTERS,  10  DIGITS  AND  BOTH  INTERCHARACTER  (42)  AND 
INTERGROUP  (33)  INTERVALS.  NOTE  THAT  DURATION  OF 
INTERCHARACTER  AND  ITNERGROUP  INTERVALS  WILL  VARY  WITH  GPM. 

THIS  WILL  SUPERCEDE  WHATEVER  IS  READ  INTIALLY  INTO 
STIMDUR. 

TIME (41, 45):  CONTAINS  IN  TIME (I, J)  FOR  EVERY  TIME  CHARACTER  I 
IS  SENT  AND  CHARACTER  J  IS  TYPED  A  RUNNING  SUM  OF  THE 
RESPONSE  TIME.  AT  THE  END,  IN  SUBROUTINE  SOUT, 

THE  AVERAGE  IS  COMPUTED  IN  TIME  BY  DIVIDING  TIME (I, J) 

BY  KEEP (I, J) . 

TIMESP(5):  TIMESP(I)  CONTAINS  RESPONSE  TIMES  OF  ALL  CORRECTLY 
RECOGNIZED  CHARACTERS  FROM  POSITION  I  IN  A  GROUP.  AT  END 
TIMESP  YIELDS  AVERAGE.  TOTAL  OBSERVATIONS  FOR  ITH  CHARACTER 
IN  STREAM  IS  KEPT  IN  KEEPSP(I). 


DEFINE  FOTTCTIONS 

SGAMMA (...):  GENERATES  PSEUDO  RANDOM 

GAMMA  RANDOM  VARIABLE  USING  ROUTINE  IN  AVERILL  M.  LAW  AND 
W.  DAVID  KELTON,  SIMULATION  MODELING  AND  ANALYSIS, 

NEW  YORK,  MCGRAW  HILL,  1982,  PAGE  257. 


DEFINE  SUBROUTINES: 

SARR (...):  THIS  CREATES  NEXT  STATE  GIVEN  MINIMUM  HAS  OCCURRED  ON 
ARRIVAL  SERVER.  ACTIONS  MUST  BE  TAKEN  TO  THE  INFO  AND 
STATE  MATRICES  FOR  BOTH  THE  ARRIVAL  AND  PERCEPTION  PROCESSES: 
A-S)  IF  THERE  ARE  MORE  ARRIVALS  CHANGE  STATEC(1,1)  TO 

STIMBUF (NELEM) ;  CHANGE  STATED (1,1)  TO  RANDOM  DURATION. 
OTHERWISE  STATECd,!)  =  STATED  (1,1)  =  0. 

A-I)  IF  THERE  ARE  MORE  ARRIVALS  DECREMENT  INFO (2,1); 

OTHERWISE  SET  INFO (1,1)  TO  ZERO. 

P-S)  IF  ARRIn/AL  is  a  character  (NOT  A  BLANK  TIME)  AND 
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IF  THERE  IS  ROOM  IN  THE  PERCEPTION  BUFFER  THEN 
SET  STATEC (INFO (1,2) +1,2)  =  OLD  STATEC(1,1);  SET 
STATED (INFO (1,2) +1,2)  =  RANDOM  DURATION.  IF 
ARRIVAL  IS  A  CHARACTER  AND  NO  ROOM  IN  RECOGNITION 
BUFFER,  PUT  DISPLACED  CHARACTER  IN  KILL  DEPOSITORY 
P-I)  IF  ARRIVAL  IS  A  CHARACTER  AND  IF  THERE  IS  ROOM  IN  THE 
PERCEPTION  BUFFER,  THEN  INFO (1,2)  =  INFO (1,2)  +  1 


SCHRBUF ( . . . ) :  PUTS  INTO  CHRBUF(410)  RANDOM  PERMUTATIONS 

OF  DIGITS  1-41  (CORRESPONDING  TO  IDENTIFIERS  FOR  MORSE 
CHARACTERS,  INCLUDING  SPECIAL  CHARACTERS) 

SEXEBUF (...):  PERFORMS  SAME  OPERATIONS  AS 
SPERBUF  IN  EXECUTION  BUFFER 

SGPM (,..):  PUTS  INTO  STIMDUR(42)  THE  INTERCHARACTER 

DURATION  FROM  INTERC(I)  WHERE  INTERC(l)  IS  6  GPM  AND 
INTERC(15)  IS  20  GPM.  SIMILARLY  FOR  STIMDUR(43)  AND 
INTERG(J) . 

SMINI (...):  IDENTIFIES  CHARACTER  WITH  MINIMUM 

REMAINING  PROCESSING  TIME  AND  SETS  KROW  AND  KCOL  TO  THE 
ROW  AND  COLUMN  ENTRY  IN  STATE  OF  THE  DURATION  OF  THIS 
CHARACTER.  FOR  EXAMPLE,  SUPPOSE: 

STATED (1,1)  =  50,  STATED (1,2)  =  100,  STATED (2, 2)  =  110, 

STATED (1,3)  =  200,  STATE (1,3)  =  120  AND  STATE (1,4)  =  80 

THEN  KROW  =  1  AND  KCOL  =  1  SINCE  50  IS  THE  MINIMUM  DURATION. 

SPER(...) :  IN  STATES  WHERE  THE  MINIMUM  DURATION  IS  IN  THE 

PERCEPTUAL  SERVER,  SPER  MAKES  THE  NECESSARY  CHANGES.  NOTE 
THIS  IS  WHERE  A  PERIOD  IS  INTRODUCED.  SPECIFICALLY, 

IF  THE  TIME  IT  TAKES  TO  PERCEIVE  A  STIMULUS  IS  LONGER 
THAN  SOME  TIME  RCRIT,  THEN  THE  STIMULUS  IS  CODED  AS 
A  PERIOD.  OTHERWISE  (I.E.,  IF  THE  SERVICE  TIME  IS 
LESS  THAN  THE  CRITICAL  TIME)  THEN  WE  USE  THE  CONFUSION 
MATRIX  TO  DETERMINE  WHETHER  THE  CHARACTER  IS  CODED 
CORRECTLY  OR  INCORRECTLY. 

SPERBUF (...):  IN  STATES  WHERE  THE  MINIMUM  DURATION 

IS  IN  THE  PERCEPTUAL  BUFFER,  SPERBUF  REMOVES  THE 
MINIMUM  DURATION  PROCESS  FROM  STATEC  AND  STATED  AND  MOVES 
EVERYTHING  ABOVE  IT  UP  ONE  IN  THE  QUEUE.  FOR  EXAMPLE,  SUPPOSE 
WE  HAVE  FIVE  ITEMS  IN  THE  PERCEPTUAL  BUFFER  AND  ONE  ON  THE 
SERVER.  WE  APPLY  SMINI  AND  OBTAIN: 

STATED (1,2)  =  50,  STATED (2, 2)  =  150,  STATED (3, 2)  =  0, 

STATED (4, 2)  =  170,  STATE (5,2)  =  120  AND  STATE (6, 2)  =  130. 

SPERBUF  MOVES  THE  THIRD  (STATED (4, 2) ) ,  FOURTH  (STATED (5, 2) ) 

AND  FIFTH  (STATED (6, 2) )  UP  ONE  AND  REPLACES  STATE(6,2)  WITH  0, 
SO  THAT  WE  HAVE: 
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STATED  (1,2)  =  50,  STATED  (2, 2)  =  150,  STATED  (3, 2)  =  170, 
STATED (4, 2)  =  120,  STATE (5, 2)  =  130  AND  STATE (6,2)  =  0. 

NOTE  THAT  THE  SAME  CHANGES  ARE  MADE  TO  STATEC 


SPERERR(JTEMPC,CONMAT,CONFUS,SEED)  :  LOOKS  AT  CONFUSION  MATRIX  AND  DECIDES 
WHETHER  TO  REPLACE  CURRENT  CHARACTER  WHICH  HAS  JUST  BEEN 
PERCEIVED  WITH  ANOTHER  CHARACTER.  FOR  EXAMPLE,  SUPPOSE  THAT 
THERE  WERE  FOUR  CHARACTERS.  SET  P(J|I)  EQUAL  TO  THE 
PROBABILITY  OF  PERCEIVING  J  GIVEN  I.  SUPPOSE  P(l|l)  =  .4, 

P(2|l)  =  .3,  P(3|l)  =  .2  AND  P  (4  ID  =  .1.  THEN,  IN  CONMAT 
IN  ROW  1  WE  CODE  .4,  .7,  .9  AND  1.0.  WE  THEN  GENERATE 
A  VALUE  FOR  A  UNIFORM [0,1]  RANDOM  VARIABLE  AND  DETERMINE 
WHETHER  THIS  VALUE  IS  IN  THE  INTERVAL  [0,.4).  IF  SO, 

THEN  THE  CHARACTER  1  IS  CORRECTLY  CODED  AS  A  1.  IF  NOT, 

THEN  WE  DETERMINE  WHETHER  THE  VALUE  IS  IN  THE  INTERVAL 
[.4,.7) .  IF  SO,  THEN  A  2  IS  INCORRECTLY  CODED  AS  A  1. 

AND  SO  ON. 

SRECBUF (STATEC, STATED, INFO, KROW,KCOL, KILL) :  PERFORMS  SAME  OPERATIONS  AS 
SPERBUF  IN  RECOGNITION  BUFFER 

SSTIM(CHRBUF,STIMBUF,NGROUP,NREP,NSTIM) :  PUTS  FIRST  CHARACTER  IN  CHRBUF 
INTO  STIMBUF,  THEN  INTERCHARACTER  IDENTIFIER,  THEN  NEXT 
CHARACTER  IN  CHRBUF,  AND  SO  ON  UNTIL  NGROUP  CHARACTERS  HAVE 
BEEN  PLACED  IN  STIMBUF,  AT  WHICH  POINT  AN  INTERGROUP 
IDENTIFIER  IS  PLACED  IN  STIMBUF.  THIS  IS  REPEATED  FOR 
EACH  SUCCEEDING  GROUP. 

SSUB  (STATEC, STATED, INFO, KROW, KCOL) :  SUBTRACTS  THE  MINIMUM  DURATION  OF  THE 
PROCESS  FOUND  IN  SMIN  FRCW  THE  DURATION  OF  EACH  OF  THE 
PROCESSES  IN  STATED,  LEAVING  THE  REMAINING  DURATION  FOR  EACH 
PROCESS.  FOR  EXAMPLE,  AFTER  SSUB  OPERATION  ON  THE  OUTPUT 
OF  SMIN  ABOVE,  WE  WOULD  HAVE: 

STATED (1,1)  =  0,  STATED (1,2)  =  50,  STATED (2, 2)  =  60, 

STATED (1,3)  =  150,  STATE (1,3)  =  120  AND  STATE (1,4)  =  30. 


TYPE  VARIABLES 
INTEGER  NCHR 
INTEGER  NCHRBUF 
INTEGER  NGPM 
INTEGER  NCNT 

INTEGER  KROW, KCOL  !ROW  AND  COLUMN  MARKERS  OF  MIN  DUR  IN  STATE 

INTEGER  NELEM 

INTEGER  NGROUP, NREP,NSTIM 

INTEGER  NOUT 

INTEGER  NPER 

INTEGER* 4  SEED 

REAL  RCRIT 
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TYPE  MATRICES 

INTEGER  CHRBUF(410) 

INTEGER  CONFUS(41,41) 

INTEGER  INFO (2, 4) 

INTEGER  KEEP (41, 46) 

INTEGER  KEEPSP(5) 

INTEGER  KILL (41, 8) 

INTEGER  STATEC(30,4) 

INTEGER  STATEP(30,4) 

INTEGER  STATEX(30,4) 

INTEGER  STIMBUF(820) 

REAL  ALPHA(2,4) ,BETA(2,4) 

REAL  CONMAT(41,45) 

REAL  DATERR(41,45) ,DATLAT(41,45) 
REAL  INTERC(15) ,INTERG(15) 

REAL  STATED (30, 4) 

REAL  STATEY(30,4) 

REAL  STIMDUR(43) 

REAL  TIME (41, 45) 

REAL  TIMESP(5) 


DATA  INITIALIZATION  STATEMENTS 
DATA  SEED/98753/ 

DATA  CONFUS/1681*0/ 

DATA  INFO/ 8*0/ 

DATA  INTERC/1260., 1030., 860., 729., 624., 537., 466., 405., 353., 
C  308. ,268. ,233. ,203. ,175. ,150./ 

DATA  INTERG/2940 . , 2402 . , 2008 . , 1701 . , 1455 . , 1254 . , 1087 . , 

C  945., 824., 718., 626., 545., 473., 408., 350./ 

DATA  KEEP/1886*0/ 

DATA  KEEPSP/5*0/ 

DATA  KILL/328*0/ 

DATA  STATEC/120*0/ 

DATA  STATED/120*0./ 

DATA  STATEP/120*0/ 

DATA  TIME/1845*0./ 

DATA  TIMESP/5*0./ 

C 

C  OPEN  FILES 

OPEN (UNIT=1, FILE=' MORSE . IN' , STATUS=' OLD' , READONLY) 

OPEN (UNIT=2,FILE=' MORSE. OUT' ,STATUS=' NEW' ) 

OPEN (UNIT=3, FILE=' MORSE .DAT' , STATUS^' OLD' , READONLY) 

C 

C  READ  IN  MATRIX  INFORMATION:  SIMULATION  AND  DATA 

CALL  SINP (STATEC, STATED, INFO, STIMDUR, ALPHA, BETA, 

C  STATEX, STATEY) 

CALL  SDAT (DATERR, DATLAT, CONMAT) 
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BOUNCE  BACK  HERE  FOR  SIMULATION  USING  NEW  GPM 
'  CONTINUE 

INITIALIZE  VARIABLES 

CALL  SINV ( INFO,  NCHR, NCHRBUF , NGPM, NGROUP , NREP , NSTIM, NELEM, NCNT , 

C  NOUT,NPER,RCRIT,KCOPB,SEED) 

RANDOMIZE  INPUT  STREAM  OF  CHARACTERS 
CALL  SCHRBUF (SEED, CHRBUF, NCHR, NOUT) 

INSERT  INTERCHARACTER  AND  INTERGROUP  IDENTIFIERS  INTO  STREAM 
CALL  SSTIM (CHRBUF , STIMBUF , NGROUP , NREP , NSTIM, NOUT) 

SET  INTERCHARACTER  AND  INTERGROUP  DURATION  FOR  A  PARTICULAR  RATE  (GROUPS 
PER  MINUTE) 

CALL  SGPM(STIMDUR, INTERC, INTERG,NGPM) 

BEGIN  CREATING  STATES 

0  CONTINUE  (BOUNCE  BACK  HERE  IF  MORE  STATES 

NO  MORE  STATES  IF  NOTHING  ACTIVE  ON  SERVERS 
IF( (INFO(l,l)+INFO(l,2)+INFO(l,3)+INFO(l,4) ) .EQ.O)  GOTO  999 

OBTAIN  MINIMUM  PROCESSING  DURATION 

CALL  SMINI (STATEC, STATED, INFO, KROW, KCOL, NOUT) 

SUBTRACT  MINIMUM  PROCESSING  DURATION 

CALL  SSUB (STATEC, STATED, INFO, KROW, KCOL) 

IF  ARRIVAL  COMPLETES  FIRST  THEN: 

IF(KCOL.EQ.l)  THEN 

NELEM  =  NELEM  +  1 

CALL  SARR (STATEC, STATED, INFO, KROW, KCOL, STIMBUF , NELEM, 

C  KILL,  SEED,  STIMDUR, ALPHA, BETA, STATEX, STATEY, STATEP, 

C  NOUT , D ATLAT , KCOP  B ) 

GO  TO  100 
END  IF 

IF  CHARACTER  ON  PERCEPTION  SERVER  COMPLETES  FIRST  THEN: 

IF( (KCOL.EQ,2) .AND. (KROW.EQ.l) )  THEN 

CALL  SPER (STATEC, STATED, INFO, KROW, KCOL, STIMBUF, NELEM, 

C  KILL, CONMAT, CONFUS, SEED, ALPHA, BETA, STATEX, STATEY, 

C  STATEP, NOUT, NPER,RCRIT,DATLAT) 

GO  TO  100 
END  IF 
C 

C  IF  CHARACTER  IN  PERCEPTION  BUFFER  COMPLETES  FIRST,  THEN 
C  SQUASH  PERCEPTION  BUFFER 

IF( (KCOL.EQ.2) .AND. (KROW.GT.l) )  THEN 

CALL  SPERBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, STATEX, 

C  STATEY, STATEP , NOUT) 

GO  TO  100 
END  IF 
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C  IF  CHARACTER  ON  RECOGNITION  SERVER  COMPLETES  FIRST  THEN: 
IF((KC0L.EQ.3) .AND. (KROW.EQ.l) )  THEN 

CALL  SREC (STATEC, STATED, INFO, KROW, KCOL, STIMBUF, NELEM, 

C  KILL, SEED, ALPHA, BETA, STATEX, STATEY, STATEP, NOUT, 

C  DATLAT) 

GO  TO  100 
END  IF 
C 

C  IF  CHARACTER  IN  RECOGNITION  BUFFER  COMPLETES  FIRST,  THEN 
C  SQUASH  RECOGNITION  BUFFER 

IF((KCOL.EQ.3) .AND. (KROW.GT.l))  THEN 

CALL  SRECBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, STATEX, 

C  STATEY, STATEP , NOUT) 

GO  TO  100 
END  IF 

IF  CHARACTER  ON  EXECUTION  SERVER  COMPLETES  FIRST  THEN: 

IF( (KCOL.EQ.4) .AND. (KROW.EQ.l) )  THEN 

CALL  SEXE (STATEC, STATED, INFO, KROW, KCOL, STIMBUF, NELEM, 

C  KEEP , KILL, ALPHA, BETA, STATEX, STATEY, TIME , STATEP , KEEPSP , 

C  TIMESP, NOUT, DATLAT) 

GO  TO  100 

END  IF 

IF  CHARACTER  IN  EXECUTION  BUFFER  COMPLETES  FIRST,  THEN 
SQUASH  EXECUTION  BUFFER 
IF ( (KCOL.EQ.4) .AND. (KROW.GT.l) )  THEN 

CALL  SEXEBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, 

C  STATEX, STATEY, STATEP , NOUT) 

GO  TO  100 

END  IF 

GET  HERE  WHEN  NO  MORE  STATES 
'9  CONTINUE 


WRITE  SELECTED  OUTPUT 

CALL  SOUT ( INFO, KEEP , KEEPSP , KILL, CONFUS, TIME, TIMESP , NCHR, NCNT, 
C  NOUT, NGPM, NPER, DATERR, DATLAT,  RCRIT,  KCOPB) 

ZERO  INPUT  TO  START  SIMULATION  AT  NEW  RATE 

CALL  SZERO (CONFUS, NCHR, INFO, NCNT, KEEP, KEEPSP, KILL, STATEC, 

C  STATEP, STATEX, STATED, STATEY, TIME, TIMESP, NPER) 

GO  TO  99 


END  MAIN 
END 
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SUBROUTINE  SINV 


8 

10 

12 

14 

16 

18 

C 


C 

C 

C 


C 


SUBROUTINE  SINV ( INFO, NCHR, NCHRBUF, NGPM, NGROUP , NREP , NSTIM, NELEM, 
C  NCNT, NOUT, NPER, RCRIT, KCOPB,  SEED) 

INTEGER  KCOPB 

INTEGER  NCHR 

INTEGER  NCHRBUF 

INTEGER  NCNT 

INTEGER  NGPM 

INTEGER  NELEM 

INTEGER  NGROUP , NREP , NSTIM 

INTEGER  NOUT 

INTEGER  NPER 

INTEGER* 4  SEED 

INTEGER  INFO (2, 4) 

REAL  RCRIT 


FORMAT ('  ENTER  1  TO  CONTINUE,  0  TO  STOP  =  ',$) 
FORMAT {'  ENTER  ELEMENT  IN  STREAM:  0  =  ',$) 
FORMAT ('  ENTER  GROUPS  PER  MINUTE:  6-20  =  ',$) 
FORMAT ('  DO  YOU  WANT  FULL  OUTPUT:  1/0  =  ',$) 
FORMAT ('  ENTER  CRITICAL  VALUE  IN  MS:  =  ',$) 
FORMAT ('  ENTER  1  FOR  COPY  BEH;  0  NO  COPB  =  ',$) 


BEGIN  MAIN  BODY  OF  PROGRAM 
WRITE (5,*) 

WRITE (5, 8) 

READ (5,*)  NSTOP 
IF(NSTOP.EQ.O)  THEN 
WRITE (5,*) 

WRITE (5,*)  '  OUTPUT  IN  FILE  MORSE. OUT  ' 

WRITE (5,*) 

STOP 
END  IF 

JNCNT  =  (INFO (2,1)  -l)/2  +  1 
WRITE (5,*) 

WRITE (5,*)  '  NUMBER  OF  STIMULI  IN  INPUT  STREAM  =  ', INFO (2,1) 

WRITE (5,*)  '  NUMBER  OF  CHARACTERS  IN  INPUT  STREAM  =  ', JNCNT 

—  NCNT 

NCNT  =  INFO (2,1) 

NELEM 
NELEM  =  0 
NGPM 
WRITE (5,*) 

WRITE (5, 12) 

READ (5,*)  NGPM 

IF( (NGPM.LT.6) .OR. (NGPM.GT.20) )  THEN 

WRITE (5,*)  '  GROUP  SIZE  INCORRECT;  STOP  PROGRAM  ' 

STOP 
END  IF 

—  RCRIT 
WRITE (5,*) 
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WRITE (5, 16) 

READ (5,*)  RCRIT 

C - KCOPB 

WRITE (5,*) 

WRITE (5, 18) 

READ (5,*)  KCOPB 

C - NOUT 

WRITE (5,*) 

WRITE(5,14) 

READ (5,*)  NOUT 

C  -  NGROUP 

NGROUP  =  5 

C - NREP 

NREP  =  50 

C - NSTIM 

NSTIM  =  2*NGROUP*NREP 

C  CHECK  THAT  THERE  ARE  NOT  TOO  MANY  STIMULI  FOR  PROGRAM  TO  HANDLE 


IF(INFO(2,l) .GT. NSTIM)  THEN 

WRITE (5,*)  '  NUMBER  OF  ARRIVING  STIMULI  TOO  LARGE:  ' 
WRITE (5,*)  '  INFO(2,l)  GT  NSTIM  =  NREP*NGROUP  ' 
STOP 
END  IF 


C - NCHR 

NCHR  =31 

C  -  NCHRBUF 

NCHRBUF  =410 

C  CHECK  THAT  THERE  ARE  NOT  TOO  MANY  STIMULI  THAT  MUST  BE  CONSTRUCTED 


IF (NGROUP*NREP . GT . NCHRBUF)  THEN 

WRITE (5,*)  '  NEED  TO  ENLARGE  MATRICES  CHRBUF  AND  STIMBUF  ' 
WRITE (5,*)  '  SO  THAT  NGROUP*NREP  GT  NCHRBUF  WHERE  CHRBUF  ' 
WRITE (5,*)  '  HAS  DIMENSION  NCHRBUF  AND  STIMBUF  HAS  DIMEN-  ' 
WRITE (5,*)  '  SION  2*NCHRBUF  ' 

STOP 
END  IF 

C - NPER 

NPER  =  0 

RETURN 

END 
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SUBROUTINE  SINP 
C 
C 
C 
C 

c 

SUBROUTINE  SINP (STATES, STATED, INFO, STIMDUR, ALPHA, BETA, 

C  STATEX, STATEY) 

INTEGER  INFO (2, 4) 

INTEGER  STATES ( 30, 4) , STATEX (3 0,4) 

REAL  ALPHA (2, 4) , BETA (2, 4) 

REAL  STATED (30, 4) , STATEY (30, 4) 

REAL  STIMDUR (43) 

10  FORMAT(1X,4I4,10F10.4) 

12  FORMAT(20(1X,12F5.0/)  ) 

14  FORMAT (8 (IX, F7. 2)) 

S 

S  BEGIN  INPUTTING  INFORMATION 
DO  100  lA  =  1,2 
READd,*)  (INFOdA,  J)  ,  J=l,4) 

100  SONTINUE 

DO  240  lA  =  1,2 

READd,*)  (ALPHA(IA,J),J=1,4),  (BETA(IA,  J) ,  J=l,  4) 

240  SONTINUE 

DO  110  lA  =  1,30 
DO  110  IB  =  1,4 
STATES (lA, IB)  =  0. 

STATED (lA, IB)  =  0. 

110  SONTINUE 

STATES (1,1)  =  42  UNITIALIZE  OTH  STIMULUS  TO  42 

STATED (1,1)  =  50.  IWITH  A  DURATION  OF  50  MS 

DO  112  lA  =  1,30 

DO  112  IB  =  1,4 

STATEX (lA, IB)  =  STATES (lA, IB) 

STATEY (lA, IB)  =  STATED (lA, IB) 

112  SONTINUE 

DO  120  lA  =  1,42,3 

READd,*)  (STIMDUR(J)  ,  J=IA,IA+2) 

120  SONTINUE 

READd,*)  STIMDUR (43) 

S 

S  OUTPUT  INFORMATION  TO  FILE  ON  WHAT  WAS  READ  IN 
WRITE (5,*) 

WRITE (5,*)  '  WARNING  -  WARNING  :  STIM  DURATION  OF  SPES  SHAR  WRONG  ' 
WRITE (5,*) 

WRITE (2,*)  '  STIMULUS  DURATIONS  OF  1  -  31  SHARASTERS  AND  BLANKS  ' 
WRITE (2, 12)  (STIMDUR (J) ,J=1, 33) 

WRITE (2,*) 

WRITE (2,*)  '  DATA  IN  ALPHA  AND  BETA  AT  START  ' 

DO  260  lA  =  1,2 

WRITE (2, 14)  (ALPHA (IA,J) ,J=1,4) , (BETA (lA, J) , J=l, 4) 

260  SONTINUE 
WRITE (2,*) 

WRITE(2,*)  '  PERSEPTION  (SERV/DES) :  ' ,ALPHA(1,2) *BETA(1,2) , 
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C  ALPHA(2,2) *BETA(2,2) 

WRITE (2,*)  '  PERCEPTION  (SERV/DEC) :  ' , ALPHA (1, 3) *BETA (1, 3) , 

C  ALPHA(2,3)*BETA{2,3) 

WRITE{2,*)  '  PERCEPTION  (SERV/DEC):  ' , ALPHA (1, 4) *BETA (1, 4) , 

C  ALPHA(2,4) *BETA(2,4) 

WRITE (2,*) 

WRITE (2,*)  '  DATA  IN  INFO  AT  START  ' 

DO  220  lA  =  1,2 

WRITE  (2,*)  (INFOdA,  J)  ,  J=l,4) 

220  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  INPUT  IN  STATEC  AND  STATED  AT  START  ' 

DO  200  lA  =  1,5 

WRITE (2, 10)  (STATEC (lA,  J)  ,J=1, 4) , (STATED (lA, J) , J=l, 4) 

200  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  INPUT  IN  STATEX  AND  STATEY  AT  START  ' 

DO  210  lA  =  1,5 

WRITE (2, 10)  (STATEX (lA, J) ,J=1,  4) ,  (STATEY (lA, J) , J=l, 4) 

210  CONTINUE 

CHECK  ON  VALIDITY  OF  INPUT  INFORMATION 
DO  400  lA  =  1,2 
DO  400  IB  =  1,4 
IF(ALPHA(IA,IB) .LE.l)  THEN 

WRITE (5,*)  '  ALPHA  LESS  THAN  1;  ALTERNATIVE  GAMMA  GENERATOR  ' 
WRITE (5,*)  '  IS  NEEDED;  SEE  LAW  AND  KELTON,  PAGE  255  ' 

STOP 
END  IF 
400  CONTINUE 
RETURN 
END 
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SUBROUTINE  SDAT 
C 
C 
C 

c 

c 

SUBROUTINE  SDAT (DATERR, DATLAT, CONMAT) 

REAL  CONMAT (41, 45) 

REAL  DATERR (41, 45) , DATLAT (41, 45) 

10  FORMAT ( IX, lOFl 0.4) 

DO  100  lA  =  1,41 

READ (3, 10)  (DATERR (lA, J) ,J=1, 45) 

100  CONTINUE 

DO  200  lA  =  1,41 

READ (3, 10)  (DATLAT (lA, J) ,J=1, 45) 

200  CONTINUE 

DO  300  lA  =  1,41 

RDEN  =  DATERR (lA,  42)  +  DATERR (lA, 43) 
IF(RDEN.GT.O)  THEN 

IF(IA.EQ.l)  THEN 

CONMAT (IA,1)  =  DATERR (lA, 42) /RDEN 
CONMAT (lA, 2)  =1,0 

ELSE 

CONMAT (IA,1)  =  DATERR (lA, 43) /RDEN 
CONMAT (IA,IA)  =1.0 
END  IF 
END  IF 

300  CONTINUE 
RETURN 
END 
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SUBROUTINE  SCHRBUF 
C 
C 
C 
C 
C 

SUBROUTINE  SCHRBUF (SEED,CHRBUF,NCHR,NOUT) 

INTEGER  NOUT 

INTEGER  NCHR 

INTEGER* 4  SEED 

INTEGER  CHRBUF(410) 

INTEGER  ICHR(41) 

REAL  RCHR(41) 

10  FORMAT (2 (IX, 1614, /) ) 

12  F0RMAT(3(1X,12F6.4,/) ) 

DO  900  IX  =  1,10 
DO  100  lA  =  1,NCKR 
RCHRdA)  =  RAN  (SEED) 
iOO  CONTINUE 

DO  200  lA  =  1,NCHR 
RMIN  =  2. 

DO  210  IB  =  1,NCHR 
IF(RCHR(IB) .LT.RMIN)  THEN 
JPOS  =  IB 
RMIN  =  RCHR(IB) 

END  IF 

210  CONTINUE 

ICHR(JPOS)  =  lA 
RCHR(Jt'OS)  =  2. 

200  CONTINUE 

DO  300  lA  =  1,NCHR 
JA  =  NCHR* ( IX- 1)  +  I A 
CHRBUF(JA)  =  ICHR(IA) 

300  CONTINUE 
900  CONTINUE 
RETURN 
END 
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SUBROUTINE  SOUT 


C 

C 

C 

C 

C 

SUBROUTINE  SOUT ( INFO, KEEP , KEEPSP, KILL, TONFUS, TIME, TIMESP , NCHR, 
C  NCNT, NOUT, NGPM, NPER, DATERR, DATLAT, 

C  RCRIT,KCOPB) 

INTEGER  KCOPB 

INTEGER  NCHR,  NCNT, NOUT, NGPM, NPER 
INTEGER  INFO (2, 4) 

INTEGER  CONFUS(41,41) 

INTEGER  KEEP (41, 46) 

INTEGER  KEEPSP (5) 

INTEGER  KILL (41, 8) 

REAL  RCRIT 

REAL  DATERR (41, 45) , DATLAT (41, 45) 

REAL  DATOUT(2,8) 

REAL  PRDERR(41,45) 

REAL  PRDOUT(2,8) 

REAL  TIME  (41, 45) 

REAL  TIMESP  (5) 

10  FORMAT (IX, 1014) 

12  FORMAT (15 (IX, 14)) 

14  FORMAT (IX, 13, IX, F7, 4, IX, 10013) 

16  FORMAT (IX, 13, IX, 15F7 . 1/5X, 15F7 . 1/5X, 15F7 .1) 

18  FORMAT(1X,10F12.6) 

20  FORMAT(1X,10F12.4) 

C  INITIALIZE  MATRICES,  VARIABLES 
DO  40  lA  =  1,2 
DO  40  IB  =  1,8 
PRDOUT(IA,IB)  =  0. 

DATOUT(IA,IB)  =  0. 

40  CONTINUE 

DO  42  lA  =  1,41 
DO  42  IB  =  1,45 
PRDERR(IA,IB)  =  0. 

42  CONTINUE 
JKEEP  =  0 
JKILL  =  0 
JOFF  =  0 
JPER  =  0 
KPD  =  0 
KPB  =  0 

DO  80  lA  =  1,NCHR 
DO  80  IB  =  1,45 
IF(IB.EQ.44)  THEN 

JPER  =  KEEP (lA, 44)  +  JPER 
END  IF 

IF(IB.NE.IA)  THEN 

JOFF  =  KEEP  (lA,  IB)  +  JOFF 
END  IF 

JKEEP  =  KEEP  (lA,  IB)  +  JKEEP 
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80  CONTINUE 

DO  85  lA  =  1,NCHR 
DO  85  IB  =  1,8 

IF(IB.EQ.2)  KPB  =  KPB  +  KILL{IA,IB) 

IF(IB.EQ.6)  KPD  =  KPD  +  KILL(IA, IB) 

JKILL  =  KILL (lA, IB)  +  JKILL 

85  CONTINUE 

JNCNT  =  (NCNT  -  l)/2  +  1 
RNCNT  =  JNCNT 
JKEEPS  =  0 
DO  86  lA  =  1,5 

JKEEPS  =  KEEPSP(IA)  +  JKEEPS 

86  CONTINUE 

C  FIND  PREDICTED  ERROR  PERCENTAGES:  NEED  TO  ASSIGN  COUNTS  TO 
C  KEEP (lA, 42)  =  KEEP(IA,IA),  KEEP (lA, 43)  =  CHARACTERS 

C  TYPED  AS  JA.NE.IE,  KEEP (lA, 44)  =  PERIODS, 

C  KEEP (lA, 45)  =  NO  RESPONSES,  AND  KEEP (lA, 46)  =  KEEP(IA,1) 

C  +  ...  +  KEEP (lA, 45).  USED  AS  DENOMINATOR  FOR  ROW 

C  ERROR  PROBABILITIES  IN  PRDERR 

C  TO  BEGIN,  CONSTRUCT  KEEP (lA, 42) 

DO  400  lA  =  1,41 

KEEP (lA, 42)  =  KEEP(IA,IA) 

C  NEXT  CONSTRUCT  KEEP (lA, 43) 

JA  =  0 

DO  403  IB  =  1,41 
IF(IA.NE.IB)  THEN 

JA  =  KEEP (lA, IB)  +  JA 
END  IF 

403  CONTINUE 

KEEP  (lA,  43)  =  JA 

C  KEEP (lA, 44)  ALREADY  CONSTRUCTED  SINCE  RUNNING  COUNT  KEPT  OF  PERIODS 
C  CONSTRUCT  KEEP (lA, 45) 

DO  404  IB  =  1,8 

KEEP (lA,  45)  =  KILL (lA, IB)  +  KEEP (lA,  45) 

404  CONTINUE 

C  KEEP (lA, 46)  SET  TO  SUM  OF  ALL  KEEP  +  KILL  RESPONSES  TO  STIMULUS  lA. 
C  AT  THE  SAME  TIME  COMPUTE  KEEP  TO  ERROR  PERCENTAGES. 

DO  406  IB  =  1,45 
IF(KEEP(IA,46) .NE.O)  THEN 
RNUM  =  KEEP (lA,  IB) 

RDEN  =  KEEP (lA, 46) 

PRDERR (lA, IB)  =  RNUM/RDEN 
END  IF 

406  CONTINUE 
400  CONTINUE 

C  COMPUTE  AVERAGE  ERROR  PERCENTAGES  FOR  CHARACTER  CORRECT, 

C  CHARACTER  INCORRECT,  PERIOD  AND  NO  RESPONSE 

DO  410  lA  =  1,41 
DO  412  IB  =  42,45 

PRDOUTd,  IB-41)  =  KEEP  (lA,  IB)  +  PRDOUT  (1,  IB-41) 

412  CONTINUE 
410  CONTINUE 

DO  414  lA  =  1,4 
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PRD0UT(2,IA)  =  PRDOUTd,  lA) /RNCNT 
414  CONTINUE 

C  ALSO  NEED  TO  FIND  OBSERVED  SUMMARY  ERROR  PERCENTAGES 
DO  420  lA  =  1,41 
RSUM  =  0. 

DO  422  IB  =  42,45 
RSUM  =  DATERR(IA, IB)  +  RSUM 
422  CONTINUE 

DO  424  IB  =  42,45 
IF(RSUM.GT.O)  THEN 

DATOUTd,  IB-41)  =  DATOUT  (1,  IB-41)  +  1. 

DATOUT{2,IB-41)  =  DATOUT (2, IB-41)  +  DATERR (lA, IB) 

END  IF 

424  CONTINUE 
420  CONTINUE 

DO  426  lA  =  1,4 

IF (DATOUT (1, lA) .NE.O)  THEN 

DATOUT (2, lA)  =  DATOUT (2, lA) /DATOUT (1, lA) 

END  IF 

426  CONTINUE 
C 

C  ALSO  NEED  TO  COMPUTE  RAW  SUMS  FOR  PREDICTED  LATENCIES;  RAW  SUMS 
C  ALREADY  RECORDED  FOR  TIMEdA,J),  J  =  1,...,41.  TIME  (lA,  42) 

C  SET  TO  TIMEdA,IA);  TIME  (lA,  43)  SET  TO  SUM  OF  CHARACTER 

C  ONLY  RESPONSE  TIMES;  AND  TIME (lA, 44)  ALREADY  SET.  NOTE 

C  THERE  IS  NO  ENTRY  IN  TIME (lA, 45)  SINCE  THERE  IS  NO  TIME 

C  FOR  NO  RESPONSES. 

DO  440  lA  =  1,41 

TIME (lA, 42)  =  TIME(IA,IA) 

DO  440  IB  =  1,41 
IF(IA.NE.IB)  THEN 

TIME (lA, 43)  =  TIME (lA, IB)  +  TIME (lA, 43) 

END  IF 

440  CONTINUE 

C  FIND  WEIGHTED  AVERAGE  SIMULATED  (PREDICTED)  LATENCIES. 

DO  450  lA  =  1,41 
DO  450  IB  =  42,45 

PRDOUTd,  IB-37)  =  PRDOUT  (1,  IB-37)  +  TIME  (lA,  IB) 

450  CONTINUE 

DO  452  lA  =  5,8 

IF (PRDOUT (l,IA-4) .NE.O)  THEN 

PRDOUT (2, lA)  =  PRDOUT (1,IA) /PRDOUT (l,IA-4) 

END  IF 

452  CONTINUE 

C  COMPUTE  TIMES  IN  EACH  CELL;  CANNOT  DO  BEFORE  ABOVE  BECAUSE 
C  NEED  TO  KEEP  SUMS  IN  TIME  IN  ORDER  TO  GET  CORRECT 

C  WEIGHTED  AVERAGES. 

DO  454  lA  =  1,NCHR 
DO  454  IB  =  1,NCHR 
IF (KEEP (lA, IB) .NE.O)  THEN 
RA  =  KEEP (lA, IB) 

TIME (lA, IB)  =  TIME(IA,IB)/RA 
END  IF 

454  CONTINUE 
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C  FIND  AVERAGE  OBSERVED  LATENCIES 
DO  460  lA  =  1,41 
DO  464  IB  =  42,45 
IF(DATLAT(IA, IB) .GT.O)  THEN 

DATOUTd,  IB-37)  =  DATOUT  (1,  IB-37)  +  1. 

DATOUT(2,IB-37)  =  DATOUT (2, IB-37)  +  DATLAT (lA, IB) 

END  IF 

IF  (IB. EQ. 42)  THEN 
END  IF 

464  CONTINUE 
460  CONTINUE 

DO  466  IB  =  5,8 
IF  (DATOUTd,  IB)  .NE.O)  THEN 

DATOUT (2, IB)  =  DATOUT (2, IB) /DATOUT (1, IB) 

END  IF 

466  CONTINUE 
C 

C  WRITE  OUT  SUMMARY  INFORMATION 
RKEEPS  =  JKEEPS 
RPC  =  RKEEPS/RNCNT 
WRITE (2,*) 

WRITE (2,*)  '  COPY  BEHIND  =  1,  NO  COPY  BEHIND  =  0  ',KCOPB 
WRITE (2,*) 

WRITE (2,*)  '  GROUPS  PER  MINUTE  (GPM)  =  ',NGPM 

WRITE (2,*) 

WRITE (2,*)  '  CRITICAL  TIME  AT  PERC  SERVER  =  ',RCRIT 
WRITE (2,*) 

WRITE (2,*)  '  TOT  CHAR,  KEEP,  KILL  =  ' , JNCNT, JKEEP, JKILL 

WRITE (2,*) 

WRITE (2,*)  '  TOT  CHAR,  TOT  COR,  %  COR  =  ', JNCNT, JKEEPS, RPC 

WRITE (2,*) 

WRITE (2,*)  '  P (CORRECT),  P (BAD  CHAR),  P (PERIOD) ,  P (NO  RESP)  ' 

WRITE (2, 18)  (PRDOUT(2, J) , J=l,4) 

WRITE (2, 18)  (DATOUT(2, J) , J=l,4) 

WRITE (2,*) 

WRITE (2,*)  '  RT(CORREC),  RT (BAD  CHR) ,  RT(PERIO),  RT(NORSP)  ' 

WRITE (2,20)  (PRDOUT(2,J) ,J=5,8) 

WRITE (2, 20)  (DATOUT(2,J),J=5,8) 

WRITE (2,*) 

WRITE (2,*)  '  SERIAL  POSITION  INFORMATION:  CORRECT  KEEP  ' 

WRITE (2,*)  (KEEPSP (J) , J=l,5) 

DO  88  lA  =  1,5 
RA  =  KEEPSP (lA) 

IF (KEEPSP (lA) .NE.O)  TIMESP(IA)  =  TIMESP (lA) /RA 
88  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  SERIAL  POSITION  INFORMATION:  CORRECT  TIME  ' 

WRITE (2,*)  (TIMESP (J),J=1, 5) 

C 

C  WRITE  OUT  DETAILED  INFORMATION 
IF(NOUT,EQ.O)  GO  TO  999 
WRITE  (2,*) 

WRITE (2,*)  '  KEEP  INFORMATION  ' 

DO  90  lA  =  1,NCHR 
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WRITE  (2, 12)  (KEEP  (lA,  J) ,  J=l,45) 

90  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  KILL  INFORMATION  ' 

DO  100  lA  =  1,NCHR 

WRITE  (2, 10)  lA,  (KILL  (lA,  J)  ,  J=l,  8) 

100  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  CONFUSION  INFORMATION  FROM  SIMULATION  ' 

DO  105  lA  =  1,NCHR 
RSUM  =  0. 

DO  110  IB  =  1,NCHR 
RSUM  =  CONFUS(IA,IB)  +  RSUM 
110  CONTINUE 

RNUM  =  CONFUS(IA,IA) 

IF(RSUM.NE.O)  THEN 

RPER  =  RNUM/RSUM 

WRITE (2, 14)  IA,RPER, (CONFUS (lA, J) , J=l, NCHR) 

ELSE 

WRITE (2,*)  lA, '  RSUM  =  0:  NO  CHARACTERS  SENT  OR  ALL  BUMPED  ' 
END  IF 

105  CONTINUE 
WRITE (2,*) 

WRITE (2,*)  '  RESPONSE  TIME  INFORMATION  ' 

DO  210  lA  =  1,NCHR 

WRITE (2, 16)  lA, (TIME (IA,J) ,J=1, NCHR) 

210  CONTINUE 
999  CONTINUE 
C 

C  PERFORM  CHECKS  ON  COMPUTATIONS 
C  CHECK  1 

NTOT  =  0 

DO  300  lA  =  1,NCHR 
DO  300  IB  =  1,NCHR 
NTOT  =  CONFUS (I A, IB)  +  NTOT 
300  CONTINUE 

IF(NTOT.NE. (JKEEP  +  JKILL  -  KPD  -  KPB  -  NPER) )  THEN 
WRITE (2,*) 

WRITE (2,*)  '  CHECK  1  PROBLEMS:  CONFUS  SUM  NE  KEEP  +  KILL  -  ' 

WRITE (2,*)  '  KILLPERBUFDECAY  -  KILLPERBUFBUMP  -  #  REC  PER' 

WRITE(2,*)  '  NTOT, JKEEP, JKILL, KPD, KPB, NPER 
WRITE (5,*) 

WRITE (5,*)  '  CHECK  1  PROBLEMS:  CONFUS  SUM  NE  KEEP  +  KILL  -  ' 

WRITE (5,*)  '  KILLPERBUFDECAY  -  KILLPERBUFBUMP  -  #  REC  PER' 

WRITE(5,*)  '  ',NTOT, JKEEP, JKILL, KPD, KPB, NPER 

ELSE 

WRITE (5,*) 

WRITE (5,*)  '  CHECK  1  OK:  ' 

WRITE (5,*)  '  CONFUS  =  KEEP  +  KILL  -  KPD  -  KPB  -  NPER' 
WRITE(5,*)  '  ', NTOT,  JKEEP, JKILL, KPD, KPB, NPER 

END  IF 

C  PERFORM  CHECKS  ON  COMPUTATIONS 
C  CHECK  2 

IF (JNCNT.NE. (JKEEP+ JKILL) )  THEN 
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WRITE (2,*) 

WRITE (2,*)  '  CHECK  2  PROBLEMS:  NUMBER  OF  STIMULI  IN  STREAM 

WRITE (2,*)  '  NOT  EQUAL  TO  KEEP  +  KILL  COUNT  ' 

WRITE (2,*)  '  ',JNCNT, JKEEP,JKILL 

WRITE (5,*) 

WRITE (5,*)  '  CHECK  2  PROBLEMS:  NUMBER  OF  STIMULI  IN  STREAM 

WRITE (5,*)  '  NOT  EQUAL  TO  KEEP  +  KILL  COUNT  ' 

WRITE (5,*)  '  JNCNT, JKEEP, JKILL 

ELSE 

WRITE (5,*) 

WRITE (5,*)  '  CHECK  2  OK:  KEEP  +  KILL  =  #  CHAR  IN  INPUT  ' 
WRITE (5,*)  '  JNCNT, JKEEP, JKILL 

END  IF 

PERFORM  CHECK  ON  COMPUTATIONS 
CHECK  3 

IF (JNCNT. NE. ( JKEEPS+JKILL+JOFF) )  THEN 
WRITE (2,*) 

WRITE (2,*)  '  CHECK  3  PROBLEMS:  NUMBER  OF  STIMULI  IN  STREAM 

WRITE  (2,*)  '  NOT  EQUAL  TO  CORRECT  KILL  +  OFF  DIAG  KEEP  ' 

WRITE (2,*)  '  ', JNCNT, JKEEPS, JKILL, JOFF 

WRITE (5,*) 

WRITE (5,*)  '  CHECK  3  PROBLEMS:  NUMBER  OF  STIMULI  IN  STREAM 

WRITE (5,*)  '  NOT  EQUAL  TO  CORRECT  +  KILL  +  OFF  DIAG  KEEP  ' 

WRITE (5,*)  '  JNCNT, JKEEPS, JKILL, JOFF 

ELSE 

WRITE (5,*) 

WRITE (5,*)  '  CHECK  3  OK;  NUMBER  OF  STIMULI  IN  STREAM  ' 

WRITE (5,*)  '  IS  EQUAL  TO  CORRECT  +  KILL  +  OFF  DIAG  KEEP  ' 

WRITE (5,*)  '  JNCNT, JKEEPS, JKILL, JOFF 

END  IF 
RETURN 
END 
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SUBROUTINE  SSTIM 


SUBROUTINE  SSTIM (CHRBUF, STIMBUF, NGROUP,NREP, NSTIM, NOUT) 
INTEGER  CHRBUF (410) , STIMBUF (820) 

INTEGER  NGROUP,NREP, NSTIM 
10  FORMAT (IX, 2013) 

JCHR  =  0 

DO  100  lA  =  l,NSTIM,2*NGROUP 
DO  120  IB  =  lA, IA+2*NGROUP-l,2 
JCHR  =  JCHR  +  1 
STIMBUF (IB)  =  CHRBUF (JCHR) 

STIMBUF (IB+1)  =  42 
120  CONTINUE 

STIMBUF (lA  +  2*NGROUP  -1)  =  43 
100  CONTINUE 

IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  DATA  IN  CHRBUF  AT  START:  ONLY  FIRST  10  ' 
WRITE (2, 10)  (CHRBUF (J) ,J=1, 10) 

WRITE (2,*) 

WRITE (2,*)  '  DATA  IN  STIMBUF  AT  START:  ONLY  FIRST  20  ' 
WRITE (2, 10)  (STIMBUF (J) ,J=1, 20) 

END  IF 
RETURN 
END 
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SUBROUTINE  SMINI 


SUBROUTINE  SMINI (STATEC, STATED, INFO, KROW, KCOL, NOUT) 
INTEGER  NOUT 
INTEGER  INFO (2, 4) 

INTEGER  STATEC (30, 4) 

REAL  STATED (30, 4) 

INTEGER  KROW, KCOL 
RMIN  =  10000000. 

DO  100  lA  =  1,4 
IF(INFO(l,IA) .GT.O)  THEN 

DO  120  IB  =  l,INFO(l,IA) 

IF (RMIN. GT. STATED (IB, lA) )  THEN 
KROW  =  IB 
KCOL  =  lA 

RMIN  =  STATED (KROW, KCOL) 

END  IF 

120  CONTINUE 

END  IF 

100  CONTINUE 

IF(NOUT.EQ.l)  THEN 

WRITE (2,*)  '  MIN  DURATION  IN  STATE  KROW, KCOL, RMIN 
END  IF 
RETURN 
END 
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SUBROUTINE  SSUB 


SUBROUTINE  SSUB (STATEC, STATED, INFO, KROW, KCOL) 

INTEGER  INFO (2, 4) 

INTEGER  STATEC (30, 4) 

REAL  STATED (30, 4) 

INTEGER  KROW, KCOL 
10  FORMAT(1X,4I4,10F10.4) 

RMIN  =  STATED (KROW, KCOL)  !THE  MINIMUM  FOUND  IN  SMIN 

DO  100  lA  =  1,4 
IF(INFO(l,IA) .EQ.O)  GOTO  100 
DO  100  IB  =  l,INFO(l,IA) 

STATED  (IB,  lA)  =  STATED  (IB,  lA)  -  RMIN 
100  CONTINUE 
RETURN 
END 
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SUBROUTINE  SARR 


SUBROUTINE  SARR (STATEC, STATED,  INFO, KROW,  KCOL, STIMBUF, NELEM, 

C  KILL, SEED, STIMDUR, ALPHA,  BETA,  STATEX, STATEY, STATEP , 

C  NOUT, DATLAT, KCOPB) 

INTEGER  KROW, KCOL 
INTEGER  NELEM 
INTEGER  NOUT 
INTEGER* 4  SEED 
INTEGER  INFO (2, 4) 

INTEGER  STATEC (30,4), STATEX (30,4) 

INTEGER  STATEP (30, 4) 

INTEGER  STIMBUF(820) 

INTEGER  KILL (41, 8) 

REAL  ALPHA (2, 4) , BETA (2, 4) 

REAL  DATLAT (41, 45) 

REAL  STATED (30, 4) , STATEY (30, 4) 

REAL  STIMDUR (43) 

0  FORMAT(1X,4I4,10F10.4) 

MAKE  SURE  THAT  SARR  IS  CALLED  WHEN  IT  WAS  POSSIBLE  FOR  ARRIVAL  TO 
HAVE  HAD  MINIMUM  DURATION 
IF(INFO(l,l) .EQ.O)  THEN 

WRITE (2,*)  '  NOT  POSSIBLE  FOR  MINIMUM  TO  HAVE  BEEN  ARRIVAL  ' 
END  IF 

C  CHANGE  ARRIVAL  STATES 
JTEMPC  =  STATEC (1,1) 

JTEMPP  =  STATEP (1,1) 

IF(INFO(2,l) .EQ.O)  THEN  ! NO  MORE  ARRIVALS 

STATEC (1,1)  =  0 
STATED  (1,1)  =  0. 

STATEP (1,1)  =  0 
STATEX (1,1)  =  0 
STATEY (1,1)  =  0. 

END  IF 

IF(INFO(2,l) .GT.O)  THEN  (MORE  ARRIVALS 

STATEC (1,1)  =  STIMBUF (NELEM) 

STATED (1,1)  =  STIMDUR (STIMBUF (NELEM) ) 

IF (STATEC (1, 1) .EQ. 43)  THEN 
STATEP (1,1)  =  0 
END  IF 

IF  (STATEC  (1, 1)  .LT.  42)  THEN 

STATEP (1,1)  =  STATEP (1,1)  +  1 
END  IF 

STATEX (1,1)  =  STATEC (1,1) 

STATEY (1,1)  =  0. 

END  I'' 

C  CHANGE  ARRIVAL  INFORMATION 
IF (INFO (2,1) .EQ.O)  THEN 
INFO (1,1)  =  0 
END  IF 
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IF(INF0(2,1) .GT.O)  THEN 

INFO(2,l)  =  INFO(2,l)  -  1 
END  IF 

C  CHANGE  PERCEPTION  STATES 
IF(JTEMPC.GE.42)  THEN 
C 

SERVER  IS  BLANK 

GOTO  999 
END  IF 


!NO  CHANGE  NEEDED  IF  CHAR  ON 

ARRIVAL 


CODE  BELOW  USED  TO  MODEL  BEHAVIOR  OF  OPERATORS  WHO  CAN  WORK  ON 
AT  MOST  ONE  STIMULUS  AT  A  TIME,  I.E.,  OF  OEPRATORS  WHO 
CANNOT  COPY  BEHIND 
IF(KCOPB,EQ.O)  THEN 

RSTAT  =  INFO (1,2)  +  INFO (1,3)  +  INFO (1,4) 

IF(RSTAT.GT.O)  THEN  !AT  LEAST  1  SERVER  OCCUPIED 

KILL(JTEMPC,2)  =  KILL ( JTEMPC, 2)  +  1 
IF(NOUT.EQ.l)  THEN 


CHARACTER:  PER  BUFF  FULL 


WRITE (2,*) 

WRITE (2,*)  '  KILLED 


JTEMPC , KILL ( JTEMPC ,  2 ) 

END  IF 
GOTO  999 
END  IF 
END  IF 

IF(INFO(l,2) .EQ.O)  THEN  ! ROOM  ON  PER  SERVER 

STATEC(1,2)  =  JTEMPC 
STATEX(1,2)  =  JTEMPC 
STATEP(1,2)  =  JTEMPP 

STATED (1,2)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,2,1, JTEMPC) 

STATEY(1,2)  =  STATED (1,2) 

END  IF 

IF( (INFO(l,2) .GT.O) .AND. (INFO(l,2) .LT. INFO (2, 2) ) )  THEN 

!MORE  ROOM  IN  PER  BUFFER 

STATEC (INFO (1,2) +1,2)  =  JTEMPC 
STATEX (INFO (1,2) +1,2)  =  JTEMPC 
STATEP (INFO (1,2) +1,2)  =  JTEMPP 

STATED (INFO (1,2) +1,2)  =  SGAMMA (ALPHA,  BETA, DATLAT, 

C  SEED, 2, 2,1, JTEMPC) 

STATEY (INFO (1,2) +1,2)  =  STATED (1,2) 

END  IF 

IF (INFO (1, 2) .GE. INFO (2,2) )  THEN  !NO  MORE  ROOM  IN  BUFFER 

KILL (JTEMPC, 2)  =  KILL (JTEMPC, 2)  +  1 
IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILLED  CHARACTER:  PER  BUFF  FULL  ', 

C  JTEMPC, KILL (JTEMPC, 2) 

END  IF 
END  IF 

CHANGE  PERCEPTION  INFORMATION 
IF(INFO(l,2) .LT.INFO(2,2) )  THEN 
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INF0(1,2)  =  INF0(1,2)  +  1 
END  IF 

999  CONTINUE 

IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  INPUT  AFTER  SARR  CHANGES  ' 

DO  910  lA  =  1,2 

WRITE (2, 10)  (INFO(IA, J) , J=l,4) 

910  CONTINUE 

DO  900  lA  =  1,2 

WRITE  (2, 10)  (STATECdA,  J) ,  J=l,4)  ,  (STATED  (lA,  J) ,  J=l,  4) 
900  CONTINUE 

DO  902  lA  =  1,2 

WRITE (2, 10)  {STATEX(IA,J),J=1,4), (STATEY (lA, J) , J=l, 4) 
902  CONTINUE 

DO  904  lA  =  1,2 

WRITE  (2, 10)  (STATEPdA,  J),J=1,4) 

904  CONTINUE 

END  IF 
RETURN 
END 
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SUBROUTINE  SPER 


SUBROUTINE  SPER (STATEC, STATED, INFO, KROW, KCOL, STIMBUF, NELEM, 

C  KILL, CONMAT, CONFUS, SEED, ALPHA, BETA, STATEX, STATEY, STATEP, 

C  NOUT,NPER,RCRIT,DATLAT) 

INTEGER  KROW, KCOL 
INTEGER  NELEM 
INTEGER  NOUT 
INTEGER  NPER 
INTEGER* 4  SEED 
INTEGER  CONFUS (41,41) 

INTEGER  INFO (2, 4) 

INTEGER  STATEC (30,4), STATEX (30, 4) 

INTEGER  STATEP (30, 4) 

INTEGER  STIMBUF (820) 

INTEGER  KILL (41, 8) 

REAL  RCRIT 

REAL  ALPHA (2, 4) , BETA (2, 4) 

REAL  DATLAT(41,45) 

REAL  STATED (30, 4) , STATEY (30, 4) 

REAL  CONMAT (41, 45) 

FORMAT(1X,4I4,10F10.4) 

MAKE  SURE  THAT  SPER  IS  CALLED  WHEN  IT  WAS  POSSIBLE  FOR  PERCEPTION  TO 
HAVE  HAD  MINIMUM  DURATION 
IF(INFO(l,2) .EQ.O)  THEN 

WRITE (2,*)  '  NOT  POSSIBLE  FOR  MINIMUM  TO  HAVE  BEEN  PERCEPTION  ' 
END  IF 

INITIALIZE  VARIABLES 
JTEMPC  =  STATEC (1,2) 

JTEMPX  =  STATEX (1,2) 

JTEMPP  =  STATEP (1,2) 

RTEMPY  =  STATEY (1,2) 

IF  THE  PERCEPTION  TIME  IS  LONG,  THEN  A  PERIOD  WILL  BE  RECORDED 
IF (RTEMPY.GT. RCRIT)  THEN 
JTEMPC  =  44 
NPER  =  NPER  +  1 
END  IF 

IF  PERCEPTION  TIME  IS  LESS  THJ^  RCRIT,  THEN  GOTO  CONFUSION  MATRIX 
TO  CREATE  PERCEPTUAL  ERRORS,  E.G.,  IF  THE  CODE 
IS  ON  THE  PERCEPTION  SERVER,  THEN  WTIH  SOME  PROBABILITY 
THE  CODE  MIGHT  BE  PERCEIVED  AS,  SAY, 

IF (RTEMPY. LE. RCRIT)  THEN 

CALL  SPERERR (JTEMPC, CONMAT, CONFUS, SEED) 

END  IF 

C  CHANGE  PERCEPTION  STATES 

IF(INFO(l,2) .EQ.l)  THEN  ! NOTHING  IN  PERCEPTION  BUFFER 

STATEC (1,2)  =  0 
STATED (1,2)  =  0. 

STATEX (1,2)  =  0 
STATEP (1,2)  =  0 


C  -  30 


STATEY(1,2)  =  0. 

END  IF 

IF(INFO(l,2) .GT.l)  THEN  !MORE  IN 

PERCEPTION  BUFFER 

STATEC(1,2)  =  STATEC(2,2)  ! ASSIGN  DURATION  PERC  SERVER 

STATED (1,2)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,2,2, JTEMPX) 

STATEX(1,2)  =  STATEX(2,2) 

STATEP(1,2)  =  STATEP(2,2) 

STATEY(1,2)  =  STATEY(2,2)  +  STATED (1,2) 

DO  100  lA  =  2,INFO(l,2) 

STATEC(IA,2)  =  STATEC (IA+1, 2) 

STATED (lA, 2)  =  STATED (IA+1, 2) 

STATEX(IA,2)  =  STATEX (IA+1, 2) 

STATE?  (lA,  2)  =  STATE?  (IA-^1,2) 

IF (IA.LT.1NF0(1,2) )  THEN 

STATEY(IA,2)  =  STATEY (IA+1, 2)  +  STATED (1,2) 

ELSE 

STATEY (lA, 2)  =  0. 

END  IF 

100  CONTINUE 

END  IF 

C  CHANGE  PERCEPTION  INFORMATION 

INFO(l,2)  =  INFO(l,2)  -1 
C  CHANGE  RECOGNITION  STATES 

IF (INFO (1, 3) .LT. INFO (2,3) )  THEN  IMORE  ROOM  IN  RECOG  BUFFER 

STATEC (INFO (1,3) +1,3)  =  JTEMPC  ! ASSIGN  DURATIONS 
STATEX(INFO(l,3)+l,3)  =  JTEMPX 
STATEP (INFO (1,3) +1,3)  =  JTEMPP 
IF(INFO(l,3) .EQ.C)  THEN 

STATED (1,3)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,3,2, JTEMPX) 

ELSE 

STATED (INFO (1,3) +1,3)  =  SGAf^lA (ALPHA, BETA, DATLAT, 

C  SEED, 2, 3, 2, JTEMPX) 

END  IF 

STATEY (INFO (1,3) +1,3)  =  RiCMPY  +  STATED (1,3) 

END  IF 

IF (INFO (1,3) .GE. INFO (2,3) )  THEN  !NO  MORE  ROOM  IN  RECOG  BUFFER 

KILL (JTEMPX, 3)  =  KILL (JTEMPX. 3)  +  1 
IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILLED  CHARACTER:  REC  BUFF  FULL 
C  JTEMPX, KILL (JTEMPX, 3) 

END  IF 
END  IF 

C  CHANGE  RECOGNITION  INFORMATION 
IF(INFO(l,3) .LT.INFO(2,3) )  THEN 
INFO(l,3)  =  INFO(l,3)  +  1 
END  IF 

999  CONTINUE 

IF(NOUT.EQ.l)  THEN 
WRITE  (2,*) 

WRITE (2,*)  '  INPUT  AFTER  SPER  CHANGES  ' 
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DO  910  lA  = 

1,2 

WRITE (2, 10) 

(INFOdA,  J) ,  J=l, 

910 

CONTINUE 

DO  900  lA  = 

1,2 

WRITE (2, 10) 

(STATECdA,  J) ,  J= 

900 

CONTINUE 

DO  902  lA  = 

1,2 

WRITE (2, 10) 

(STATEX (lA, J) , J= 

902 

CONTINUE 

DO  904  lA  = 

1,2 

WRITE (2, 10) 

(STATE? (lA, J) ,J= 

904 

CONTINUE 

END  IF 
RETURN 
END 


4) 

a, 4) , (STATED (lA, J) ,J=1, 4) 
a, 4), (STATEy(IA,J),J=l,4) 
1,4) 
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SUBROUTINE  SREC 


SUBROUTINE  SREC (STATES, STATED, INFO, KROW, KCOL, STIMBUF , NELEM, 

C  KILL, SEED, ALPHA, BETA, STATEX, STATEY, STATED, NOUT, 

C  DATLAT) 

INTEGER  KROW, KCOL 
INTEGER  NELEM 
INTEGER  NOUT 
INTEGER* 4  SEED 
INTEGER  INFO (2, 4) 

INTEGER  STATEC(30,4) ,STATEX(30,4) 

INTEGER  STATEP(30,4) 

INTEGER  STIMBUF (820) 

INTEGER  KILL (41, 8) 

REAL  ALPHA(2,4) ,BETA(2,4) 

REAL  DATLAT (41, 45) 

REAL  STATED (30, 4) , STATEY (30, 4) 

10  FORMAT(1X,414,10F10.4) 

C  MAKE  SURE  THAT  SREC  IS  CALLED  WHEN  IT  WAS  POSSIBLE  FOR  PERCEPTION  TO 
C  HAVE  HAD  MINIMUM  DURATION 

IF(INFO(l,3) .EQ.O)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  NOT  POSSIBLE  FOR  MINIMUM  TO  HAVE  BEEN  RECOGNITION  ' 
END  IF 

C  CHANGE  RECOGNITION  STATES 
JTEMPC  =  STATEC(1,3) 

JTEMPX  =  STATEX (1,3) 

JTEMPP  =  STATEP(1,3) 

RTEMPY  =  STATEY (1,3) 

IF(INFO(l,3) ,EQ.l)  THEN  (NOTHING  IN  RECOGNITION  BUFFER 

STATEC(1,3)  =  0 
STATED (1,3)  =  0. 

STATEX (1,3)  =  0 
STATEP(1,3)  =  0 
STATEY (1,3)  =  0. 

END  IF 

IF(INFO(l,3) .GT.l)  THEN  (MORE  IN 

RECOGNITION  BUFFER 

STATEC(1,3)  =  STATEC(2,3)  (ASSIGN  DURATION  TO  RECOG  SERVER 

STATED (1,3)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,3,3, JTEMPX) 

STATEX (1,3)  =  STATEX (2, 3) 

STATED (1,3)  =  STATED (2, 3) 

STATEY (1,3)  =  STATEY (2, 3)  +  STATED (1,3) 

DO  100  lA  =  2, INFO (1,3) 

STATEC(IA,3)  =  STATEC (IA+1, 3) 

STATED (lA, 3)  =  STATED (IA+1, 3) 

STATEX (lA, 3)  =  STATEX (IA+1, 3) 

STATED (IA,3)  =  STATED (IA+1, 3) 

IF(IA.LT.INFO(l,3) )  THEN 
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r 


100 

c 

c 


c 

999 

910 

900 

902 

904 


STATEY (IA,3) 

ELSE 


STATEY {IA+1, 3)  +  STATED (1,3) 


STATEY (lA, 3)  =  0. 

END  IF 
CONTINUE 
END  IF 

CHANGE  RECOGNITION  INFORMATION 
INFO (1,3)  =  INFO (1,3)  -1 
CHANGE  EXECUTION  STATES 

IF(INFO(l,4) .EQ.O)  THEN  (ROOM  ON  EXEC  SERVER 

STATEC(1,4)  =  JTEMPC 
STATEX(1,4)  =  JTEMPX 
STATE? (1,4)  =  JTEMPP 

STATED (1,4)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,4,3, JTEMPX) 

STATEY (1,4)  =  RTEMPY  +  STATED (1,4) 

END  IF 

IF( (INFO(l,4) .GT.O) .AND. (INFO(l,4) .LT. INFO (2, 4) ) )  THEN 

!MORE  ROOM  IN  EXEC  BUFFER 


STATEC (INFO (1,4) +1,4) 
STATEX (INFO (1,4) +1,4) 
STATEP (INFO (1,4) +1,4) 
STATED (INFO (1,4) +1,4) 


STATEY (INFO (1,4) +1,4) 

END  IF 

IF(INFO(l,4) .GE.INFO(2,4) )  THEN 

KILL (JTEMPX, 4)  =  KILL (JTEMPX, 4)  4  1 
IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILLED  CHARACTER: 
C 


JTEMPC 

JTEMPX 

JTEMPP 

SGAMMA (ALPHA, BETA, DATLAT, 

SEED, 2, 4, 3, JTEMPX) 
RTEMPY  +  STATED (1,4) 


!NO  MORE  ROOM  IN  EXEC  BUFFER 


EXEC  BUFF  FULL  ', 

JTEMPX, KILL (JTEMPX, 4 ) 


END  IF 
END  IF 

CHANGE  EXECUTION  INFORMATION 
IF (INFO (1, 4) .LT. INFO (2, 4))  THEN 
INFO  (1,4)  =  INFO  (1,4)  +  1 
END  IF 
CONTINUE 

IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  INPUT  AFTER  SREC  CHANGES  ' 

DO  910  lA  =  1,2 

WRITE (2, 10)  (INFO(IA, J) ,J=1,4) 

CONTINUE 
DO  900  lA  =  1,2 

WRITE (2, 10)  (STATEC (lA, J) ,J=1, 4) , (STATED (lA, J) , J=l, 4) 

CONTINUE 

DO  902  lA  *  1,2 

WRITE (2, 10)  (STATEX(IA,J),J=1,4),  ( STATEY ( lA, J) , J=1 , 4 ) 

CONTINUE 

DO  904  lA  =  1,2 

WRITE (2, 10)  (STATEP (IA,J) ,J=1,4) 

CONTINUE 
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END  IF 
RETURN 
END 
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SUBRUOTINE  SEXE 


SUBROUTINE  SEXE (STATEC, STATED, INFO, KROW, KCOL, STIMBUF, NELEK, 

C  KEEP, KILL, ALPHA, BETA, STATEX, STATEY, TIME, STATEP , 

C  KEEPSP,TIMESP,NOUT,DATLAT) 

INTEGER  KROW, KCOL 
INTEGER  NELEM 
INTEGER  NOUT 
INTEGER  INFO (2, 4) 

INTEGER  STATEC (30,4), STATEX (30,4) 

INTEGER  STATEP (30, 4) 

INTEGER  STIMBUF (820) 

INTEGER  KILL (41, 8) 

INTEGER  KEEP (41, 46) 

INTEGER  KEEPSP(5) 

REAL  ALPHA (2, 4) , BETA (2, 4) 

REAL  DATLAT(41,45) 

REAL  STATED (30, 4) , STATEY (30, 4) 

REAL  TIME  (41, 45) 

REAL  TIMESP (5) 

0  FORMAT(1X,4I4,10F10.4) 

MAKE  SURE  THAT  SEXE  IS  CALLED  WHEN  IT  WAS  POSSIBLE  FOR  PERCEPTION  TO 
HAVE  HAD  MINIMUM  DURATION 
IF(INFO(l,4) ,EQ.0)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  NOT  POSSIBLE  FOR  MINIMUM  TO  HAVE  BEEN  EXECUTION  ' 
END  IF 

SAVE  INFORMATION  ON  CHARACTERS  WHICH  MAKE  IT  THROUGH  RESPONSE 
AND  INCREMENT  RUNING  RESPONSE  TIMES 
KEEP (STATEX (1, 4) , STATEC (1,4) )  =  KEEP (STATEX (1, 4) , STATEC (1, 4) )  +1 
TIME (STATEX (1, 4) , STATEC (1,4) )  =  TIME (STATEX (1, 4) , STATEC (1, 4) )  + 

C  STATEY (1,4) 

SAVE  SERIAL  POSITION  INFORMATION  FOR  CHARACTERS  WHICH  ARE  CORRECTLY 
ENTERED 

IF (STATEX (1,4) . EQ. STATEC (1, 4) )  THEN 

KEEPSP  (STATEP (1,4))  =  KEEPS? (STATEP (1,4))  +  1 

TIMESP (STATEP (1,4))  =  TIMESP (STATEP  (1, 4) )  +  STATEY (1,4) 

END  IF 

C  CHANGE  EXECUTION  STATES 

IF(INFO(l,4) .EQ.l)  THEN  INOTHING  IN  EXECUTION  BUFFER 

STATEC (1,4)  =  0 
STATED (1,4)  =  0. 

STATEX (1,4)  =  0 
STATEP (1,4)  =  0 
STATEY (1,4)  =  0. 

END  IF 

IF(INFO(l,4) .GT.l)  THEN  ! MORE  IN 

EXECUTION  BUFFER 
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STATEC(1,4)  =  STATEC(2,4) 

STATED (1,4)  =  SGAMMA (ALPHA, BETA, DATLAT, SEED, 

C  1,4,4, JTEMPX) 

STATEX(1,4)  =  STATEX(2,4) 

STATEP(1,4)  =  STATEP(2,4) 

STATEY(1,4)  =  STATEY(2,4)  +  STATED (1,4) 

DO  100  lA  =  2, INFO (1,4) 

STATEC(IA,4)  =  STATEC (IA+1, 4) 

STATED (lA, 4)  =  STATED (IA+1, 4) 

STATEX(IA,4)  =  STATEX (IA+1, 4) 

STATEP(IA,4)  =  STATEP(IA+1,4) 

IF(IA.LT.INFO(l,4) )  THEN 

STATEYdA,  4)  =  STATEY  (IA+1,  4)  +  STATED  (1,4) 

ELSE 

STATEY (I A, 4)  =  0. 

END  IF 

100  CONTINUE 

END  IF 

C  CHANGE  EXECUTION  INFORMATION 
INFO(l,4)  =  INFO(l,4)  -1 
IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  INPUT  AFTER  SEXE  CHANGES  ' 

DO  910  lA  =  1,2 

WRITE  (2, 10)  (INFOdA,  J) ,  J=l,4) 

910  CONTINUE 

DO  900  lA  =  1,2 

WRITE ( 2 , 1 0 )  (STATEC (I A, J) , J=1 , 4 ) , ( STATED (lA, J) , J=1 , 4 ) 
900  CONTINUE 

DO  912  lA  =  1,2 

WRITE ( 2 , 1 0 )  ( STATEX (I A, J) , J=1 , 4 ) , ( STATEY (I A, J) , J=1 , 4 ) 
912  CONTINUE 

DO  914  lA  =  1,2 

WRITE (2, 10)  (STATEP(IA,J),J=1,4) 

914  CONTINUE 

END  IF 
RETURN 
END 
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SUBROUTINE  SPERBUF 


C 

C 

C 

C 

C 

SUBROUTINE  SPERBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, STATEX, 

C  STATEX, STATE? , NOUT) 

INTEGER  INFO (2, 4) 

INTEGER  KROW, KCOL 
INTEGER  NOUT 
INTEGER  KILL (41, 8) 

INTEGER  STATEC(30,4) , STATEX (3 0,4) 

INTEGER  STATEP(30,4) 

REAL  STATED (30, 4) , STATEX (30, 4) 

C 

C  FORMAT  STATEMENTS 

10  FORMAT(1X,4I4,1X,10F10.4) 

C 

C  PLACE  DECAXED  STIMULUS  IN  KILL 

KILL (STATEX (KROW, KCOL) ,6)  =  KILL (STATEX (KROW, KCOL) , 6)  +  1 
IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILL  DECAX  IN  PERCEPTION  STATEX (KROW, KCOL) 

END  IF 
C 

C  BEGIN  BUFFER  REORGANIZATION 

DO  100  lA  =  KROW, INFO (1,2) 

STATEC (lA, 2)  =  STATEC (IA+1, 2) 

STATED (lA, 2)  =  STATED (IA+1, 2) 

STATEX (lA, 2)  =  STATEX (IA+1, 2) 

STATEP(IA,2)  -  STATE? (IA+1, 2) 

STATEX (lA, 2)  =  STATEX (IA+1, 2) 

100  CONTINUE 

C 

C  DECREMENT  PERCEPTION  INFO 

INFO(l,2)  =  INFO(l,2)  -  1 
C 

C  OUTPUT 

IF(NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  OUTPUT  IN  STATE  AFTER  PERCEPTION  BUFFER  DECAX  ' 
DO  210  lA  =  1,2 

210  WRITE(2,10)  (STATEC(IA,J),J=1,4), (STATED(IA,J),J=1,4) 

DO  212  lA  =  1,2 

212  WRITE(2,10)  (STATEX(IA,J),J=1,4),  (STATEX (IA,J),J=1,  4) 

DO  214  lA  =  1,2 

214  WRITE (2, 10)  (STATEP (lA, J) , J=l, 4) 

END  IF 
RETURN 
END 
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SUBROUTINE  SRECBUF 


C 

c 

c 

c 

c 

SUBROUTINE  SRECBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, STATEX, 

C  STATEY,STATEP,NOUT) 

INTEGER  INFO (2, 4) 

INTEGER  KROW, KCOL 
INTEGER  NOUT 
INTEGER  KILL (41, 8) 

INTEGER  STATEC (30,4), STATEX (30,4) 

INTEGER  STATEP(30,4) 

REAL  STATED(30,4) ,STATEY(30,4) 

C 

C  FORMAT  STATEMENTS 

10  FORMAT(1X,4I4,1X,10F10.4) 

C 

C  PLACE  DECAYED  STIMULUS  IN  KILL 

JA  =  STATEX (KROW, KCOL) 

KILL(JA,7)  =  KILL(JA,7)  +  1 
IF(NOUT.EQ.l'  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILL  DECAY  IN  RECOG  ', STATEX (KROW, KCOL) 

END  IF 
C 

C  BEGINNING  OF  REORDERING  OF  RECOGNITION  BUFFER 

DO  100  lA  =  KROW, INFO (1,3) 

STATEC (lA, 3)  =  STATEC (IA+1, 3) 

STATED (lA, 3)  =  STATED (IA+1, 3) 

STATEX (lA, 3)  =  STATEX (IA+1, 3) 

STATEP(IA,3)  =  STATEP (IA+1,3) 

STATEY(IA,3)  =  STATEY (IA+1, 3) 

100  CONTINUE 

C 

C  DECREMENT  RECOGNITION  INFO 

INFO (1,3)  =  INFO (1,3)  -  1 
C 

C  OUTPUT 

IF(NOUT.EQ,l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  OUTPUT  IN  STATE  AFTER  RECOGNITION  BUFFER  DECAY  ' 
DO  210  lA  =  1,2 

210  WRITE(2,10)  (STATEC(IA,J),J=1,4), (STATED(IA,J) ,J=1,4) 

DO  212  lA  =  1,2 

212  WRITE (2, 10)  (STATEX (lA, J) , J=l, 4) , (STATEY (lA, J) , J=l, 4) 

DO  214  lA  =  1,2 

214  WRITE (2, 10)  (STATEP (lA, J) , J=l, 4) 

END  IF 
RETURN 
END 
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SUBROUTINE  SEXEBUF 


C 

c 

c 

c 

c 


c 

c 

10 

c 

c 


c 

c 


100 

c 

c 

c 

c 


210 

212 

214 


SUBROUTINE  SEXEBUF (STATEC, STATED, INFO, KROW, KCOL, KILL, STATEX, 

C  STATEY, STATE? , NOUT) 

INTEGER  INFO (2, 4) 

INTEGER  KROW, KCOL 
INTEGER  NOUT 
INTEGER  KILL (41, 8) 

INTEGER  STATEC (30, 4) , STATEX (30, 4) 

INTEGER  STATE? (30, 4) 

REAL  STATED  (30,  4)  ,  STATEY  (30, 4) 

FORMAT  STATEMENTS 
FORMAT (IX, 414, lOFlO . 4) 

?LACE  DECAYED  STIMULUS  IN  KILL 

KILL (STATEX (KROW, KCOL) ,8)  =  KILL (STATEX (KROW, KCOL) , 8)  +  1 
IF (NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  KILL  DECAY  IN  EXECUTION  ', STATEX (KROW, KCOL) 
END  IF 

BEGIN  BUFFER  REORGANIZATION 
DO  100  lA  =  KROW, INFO (1,4) 

STATEC (lA, 4)  =  STATEC (IA+1, 4) 

STATED (lA, 4)  =  STATED ( IA+1 , 4 ) 

STATEX (lA, 4)  =  STATEX (IA+1, 4) 

STATE?  (lA, 4)  =  STATE? (IA+1, 4) 

STATEY (lA, 4)  =  STATEY (IA+1, 4) 

CONTINUE 

DECREMENT  EXECUTION  INFO 
INFO (1,4)  =  INFO (1,4)  -  1 

OUTPUT 

IF (NOUT.EQ.l)  THEN 
WRITE (2,*) 

WRITE (2,*)  '  OUTPUT  IN  STATE  AFTER  EXECUTION  BUFFER  DECAY 
DO  210  lA  =  1,2 

WRITE ( 2 , 1 0 )  ( STATEC ( I A, J) , J=1 , 4 ) , ( STATED ( I A, J) , J=1 , 4 ) 

DO  212  lA  =  1,2 

WRITE(2,10)  ( STATEX (IA,J),J=1, 4), (STATEY (lA, J) , J=l, 4) 

DO  214  lA  =  1,2 

WRITE (2, 10)  (STATE?(IA,J),J=1,4) 

END  IF 
RETURN 
END 


o  o  o  o  o 


SUBROUTINE  SPERERR 


100 

110 


SUBROUTINE  SPERERR (JTEMPC, CONMAT, CONFUS,  SEED) 
INTEGER  JTEMPC 
INTEGER* 4  SEED 
INTEGER  CONFUS (41, 41) 

REAL  CONMAT (41, 45) 

RA  =  RAN (SEED) 

DO  100  lA  =  1,41 

IF (RA.LE. CONMAT (JTEMPC, lA) )  THEN 

CONFUS (JTEMPC, lA)  =  CONFUS (JTEMPC, lA)  +  1 
JTEMPC  =  I A 
GOTO  110 
END  IF 
CONTINUE 
CONTINUE 
RETURN 
END 


C  -  41 


oooo  ooooo  noooo 


SUBROUT' -JE  SGPM 


SUBROUTINE  SGPM (STIMDUR, INTERC,  INTERG, NGPM) 
INTEGER  NGPM 

REAL  INTERC (15) , INTERG (15) 

REAL  STIMDUR  ('' 3) 

STIMDUR (42,  =  INTERC (NGPM- 5) 

STIMDUR (43)  =  INTERG (NGPM- 5) 

RETURN 

END 


REAL  FTn^CTION  SGAMMA (ALP HAl, BETA1,DATLAT, SEED, JROW, JCOL, 

C  JPRO, JCHR) 

INTEGER* 4  SEED 
REAL  ALPHA, BETA 
REAL  ALPHAl (2,4) ,BETA1 (2,4) 

REAL  A,B,D,Q,U1,U2,THETA,W,V,Y,Z 
REAL  DATLAT(41,45) 

ALPHA  =  ALPHAl (JROW,JCOL) 

BETA  =  BETAl (JROW, JCOL) 

ASSIGNING  PERCEPTION  TIME:  PERCEPTION  TIME  *  OBSERVED  CORRECT  TIME  - 
RECOGNITION  TIME  -  EXECUTION  TIME.  BETA (1,2)  (VARIANCE  FACTOR) 
ASSIGNED  BEFOREHAND.  ONLY  A  MATTER  OF  ASSIGNING  ALPHA 
SO  THAT  ABOVE  EQUAITON  HOLDS 
IF( (JROW.EQ.l) .AND. (JCOL.EQ.2))  THEN 

RSUB  =  ALPHAl  (1,3)  *BETA1  (1,3)  -t  ALPHAl  (1,  4)  *BETA1  (1, 4) 

ALPHA  =  (DATLAT(JCHR,42)  -  RSUB) /BETA 
END  IF 

A  =  l./((2.*ALPHA  -  1.)**.5) 

B  =  ALPHA  -  LOG (4.) 

Q  =  ALPHA  +  l./A 
THETA  =4.5 
D  =  1.  +  LOG (THETA) 

10  Dl  =  RAN (SEED) 

IF(Ul.LE.O)  GO  TO  10 
U2  =  RAN (SEED) 

V  =  A*LOG(Ul/ (1.  -  Ul) ) 

Y  =  ALPHA*EXP(V) 

Z  =  U1*U1*U2 

W  =  B  +  Q*V  -  Y 

IF({W  +  D  -  THETA*Z) .GE.O)  THEN 
SGAMMA  =  Y*BETA 
RETURN 
END  IF 

IF (W.GE.LOG (Z) )  THEN 
SGAMMA  =  Y*BETA 


C  -  42 


RETURN 


ELSE 

GO  TO  10 
END  IF 
RETURN 
END 


C  -  43 


SUBROUTINE  SZERO 


C 

c 

c 

c 

c 

SUBROUTINE  SZERO (CONFUS, NCHR, INFO, NCNT, KEEP, KEEPSP, KILL, 

C  STATEC, STATEP , STATEX, STATED, STATEY, TIME, TIMESP , 

C  NPER) 

INTEGER  NCHR 
INTEGER  NCNT 
INTEGER  NPER 
INTEGER  CONFUS (41, 41) 

INTEGER  INFO (2, 4) 

INTEGER  KEEP (41,  46) 

INTEGER  KEEPSP (5) 

INTEGER  KILL (41, 8) 

INTEGER  STATEC (30,  4) 

INTEGER  STATEP (30, 4) 

INTEGER  STATEX (30, 4) 

REAL  STATED (30, 4) 

REAL  STATEY (30, 4) 

REAL  TIME(41,45) 

REAL  TIMESP  (5) 

C  —  NPER 
NPER  =  0 
C  — >  INFO 

INFO(l,l)  =  1 
INFO (2,1)  =  NCNT 
C  — —  CONFUS 

DO  100  lA  =  1,NCHR 

DO  100  IB  =  1,NCHR 

CONFUS (lA, IB)  =  0 
100  CONTINUE 
C  —  KEEP 

DO  105  lA  =  1,NCHR 

DO  105  IB  =  1,45 

KEEP  (lA,  IB)  =  0 
105  CONTINUE 

C  -  KEEPSP 

DO  110  lA  =  1,5 

KEEPSP (lA)  =  0 
110  CONTINUE 
C  KILL 

DO  120  lA  =  1,NCHR 

DO  120  IB  *  1,8 

KILL (lA, IB)  -  0 
120  CONTINUE 

C  -  STATEC, STATEP, STATEX 

DO  130  lA  -  1,30 

DO  130  IB  -  1,4 

STATEC (lA, IB)  -  0 
STATEP (lA, IB)  -  0 
STATEX (lA, IB)  -  0 


C  -  44 


130  CONTINUE 

C  IN  ORDER  TO  START  SIMULATION,  WANT  FIRST  CHARACTER  IN  ARRIVAL  TO 
C  BE  A  BLANK,  42  OR  43  IS  OK 

STATEC(1,1)  =  42 
STAIEX(1,1)  =  42 

C -  STATED,  STATEY 

DO  140  lA  =  1,30 
DO  140  IB  =  1,4 
STATED  QA,  IB)  =  0. 

STATEY (lA, IB)  =  0. 

140  CONTINUE 

C  IN  ORDER  TO  START  SIMULATION,  NEED  TO  ASSIGN  A  DURATION  TO  BLANK 
C  SO  THAT  SMINI  DOESN'T  THINK  EVERYTHING  IS  FINISHED 

STATED (1,1)  =  50. 

STATEY (1,1)  =  50. 

C -  TIME 

DO  150  lA  =  1,NCHR 

DO  150  IB  =  1,45 

TIME  (lA,  IB)  =  0. 

150  CONTINUE 

C -  TIMES? 

DO  160  lA  =  1,5 

TIMES? (I A)  =  0. 

160  CONTINUE 
RETURN 
END 


i 


C  -  45 


